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A 3 S 7 RACT 

Methods for obtaining iouiratiou rate contours, fractionation 
caters, atom-concentration is toss it;.- ratios, and particle site distribution 
fron the Ml Her fallout sod el are given. Included are complete computer 
profrass for estimating sublimation pressures, ionisation rate contours, 
and particle site pa rase t e rr. 

Graphical integration of activity iron fallout over the watersheds 
serving Houston, Texas asd iTev 'fork City is described. Maximum levels of 
selected contaminants at the water intake for these cities ur.der nost adverse 
wind conditions were calculated and reported In yc/ul as follows: Mew York 
City; Sr~=9, £-3 x 10' 3 t Sr-90, 5-9 x 10* 5 ; Ku-106, 7-1 x 10‘ 4 ; 1-131, 

7-6 x 10' 2 ; Cs-137, x 10" 5 ; Ea-li-O, 5-3 x lO* 2 . Houston; Sr- 89 , 1.9 x 10‘ 2 j 
Sr- 90 , 1.9 X 10' 4 ; Hu-10c, 2.5 X 10* 3 ; 1-131, 2-5 x 10* 1 ; Cs-137, 1.0 x 10' 4 ; 
Ba-lk 1 , 1.9 X 10* 1 . 

It Is concluded that the contribution of Induced radioactivity 
to the contf—l.r.stlon of water supplies would not be significant. 

Many commercially available Instruments of the survey type can 
be used to detect adequately the presence of substantial concentrations 
of radioisotopes In water under emergency conditions. For more accurate 
measurement of such concentrations down to safe levels more sensitive 
equipment or methods of preconcentration of the samples are needed. Such 
methods or instrument* should extend the range of sensitivity by a factor 
of the order of 100 and should be reaiily adaptable for field use. Outlines 
of procedures for the analysis of radioactivity in water are given. Efficiency 
of water decontamination by various methods Is discussed. 

Boiy burdens of individual radioisotopes nay be determined from 
a consideration of the rate of ingestion, the effective decay rate, and the 


.VUVW.v 





affinity for the Isotope of the critical organ. It Is heller**} that 
uptake of radioactive Isotopes nay he reduced or prevented by selective 
blocking of critical ore&ns vlth stable Isotopes. 
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r. rrrrpoDucnon 


The purpose of this study was to evaluate and s-.immarize available 
Information on the problem of water contamination bv radioactive fallout 
in the event of nuclear war. Consideration has been given to the current 
theories of the formation and distribution of fallout and the level of 
fallout that might result from a possible nuclear attack. 

(l ? ) 

Of the various attack models that have beer, recently described' ' 

(d' 

the or.e prepared by Technical Operations, Inc. ' hn3 been selected aa 


most appropriate for the present study. This model appears to have seen 


well thought out and carefully developed, giving due consideration to the 


relative Importance of various military, industrial, governmental, and 
power resource targets. This model assumes all detonations to be surface 
bursts, although perhaps a more realistic situation would include many 
air bursto, which would be more effective against an unhardened military 
or industrial target. As the type of burst has a considerable Influence 
on the fallout produced. It must be carefully evaluated In a study at 
water contamination. 

(M 

A thorough analycls of the fallout model of Miller has been 
made. Several Important functions derived from this model have been 
utilized In this study and are discuGscd in detail in Section II, "Summary 
of Analysis of the Fallout Model". 

Information from the Miller model has ulco been used in the 
study of possible levels of contamination by selected icotopco In the water 
cupply systems of Houston, Texas and Hew ’fork City. With appropriate con¬ 
sideration to the assumptions made, these results may he applied to the 
watersheds of municipalities other than those shown here. 
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This 3tudy Is corvee reed with the effect of Individual nuclear 
blasts, and does not consider the effects of overlapping fallout patterns. 

The complex problem of overlap has been treated by Technical Operations, 
Inc.^ and others, but additional studies on this subject appear desirable. 

Transport of fallout by surface eater has been studied. Evslua- 
tion of the effects of watershed characteristics on the concentration of 
radioactivity in surface water would be valuable. 

Decontamination procedures, both standard and emergency, have 
beer, described previously. ^ See Appendix A) These procedures are evaluated 
according to their efficacies for the removal of various radioisotopes and 
nlfo according to the feasibility of their use after a nuclear attack. 

For purposes of this study, it has been assumed that waterworks facilities 
have suffered only minimal damage, and that personnel will be available 
for operation. 

A study of published methods of radiochemical analysis and 
existing Instrumentation has been made. Due to the question of availability 
or personnel and facilities it appear* that this information will be of more 
value in recognizing long-term hazards than in determining potability la 
the immediate post-attack period. 

‘Ric biological uptake and resultant body burdens in man have been 
evaluated in relation to tho levels of water contamination considered likely. 
Existing deta haa been put into a workable form so *hat the body burden 
of any of the selected biologically important radioisotopes may be determined 
at any specific time after the initiation of intake. 

"The internal radiation hazard resulting from nuclear attack is 
of importance only to those persons who survive the catastrophe for a 
period of several years. As a primary onnihilator, it is of no consequence. 
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as the external dose concomitant with a lethal Internal dose la more deadly 
by several orders of magnitude. In general, radiation from sub-lethal 
concentrations of internally deposited radionuclides act3 slowly, and 
its results do not becore apparent for years. Therefore, the internal 
radiation hazard is not among the forces debilitating the population and 
reducing its Immediate ability to cope with a radically altered environ¬ 
ment. Rather it has long-term effects which culminate in producing diseases 
such as cancer, shortening of the life spaa, and possibly affecting the 
gem plasm of mankind to a small degree. 


W1 VVA'-VVVSiVV ..Ilklt.'Y-V -V 
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II. SUMMARY OF ANALYSIS OP TOE O.C.D. FALLOUT MODEL 

A. Ba3lc Assumptions of the O.C.D. Model 

During the initial phase of this project, the Postattack Research 
Division, Office of Civil Defense, requested that all studies be based on the 
fallout model as developid in "Fallout and Radiological Countermeasures"^. 

A thorough and quantitative study of the model vas conducted because the 
evaluations of vater contamination depend directly on previously established 
information concerning the distribution and properties of local fallout. The 
results of the detailed study are now presented In summary form. 

The model was based on both established physico-chemical theories 
and actual test data, brought together by the scaling method- The rationality 
and simplicity of tho basic assumptions provide this model with many advan¬ 
tages. The basic assumptions include: 

(1) The model applies generally to a "ground surface" burst, but 
it con be extended to other types of detonations by suitable transformations. 

(2) The required Input data consist of the weapon yield and the 
wind velocity. 

(3) The radioactive decay schemes, chemical properties and rela¬ 
tive abundances of the fission products are based on the data developed by 
Bolies & Ballou^, Kntcoff^, and Miller and Loeb^. 

(4) The Increase of the fireball volume with time Is beGed on 
the adiabatic expansion of an Ideal gas with the thermodynamic equations 
modified to Include a term for the change in free energy with altitude and 
with the external preooure proportional to exp (-mgZ/kt) . 


* For explanation of symbols, cce Glossary. 








(j) The radioective cloud is assumed to have the shape of an 
oblate spheroid; hence the Ionization rate contours or fallout patterns 
under a constant wind velocity are generally cigar-shaped.. 

(6) All values for the dimensions and properties of the fireball 
and fallout patterns are scaled from test data by the scaling method to he 
a function of weapon yield and wind velocity. 

(7) The condensation process of fallout is divided into two 
time periods: (a) the first period is characterized by the presence of 
gas and liquid phases and cuds when the bulk carrier or substrate material 
of the particles solidifies, nnd (b) the second period i3 characterized 

by the existence of ga3 aad solid phases. 'Hie temperature of l673°K Is 
used to divide these two periods of condensation. The radioactive material 
vhlch condenses In the first period will be fused inelde the fallout parti¬ 
cles and generally becomes Insoluble in water. The material which condenses 
during the second period will be adsorbed on the outside of the fallout 
particle and becomes readily soluble in water. With these basic assump¬ 
tions, the model will yield the intensity, the activity and atom-concentra¬ 
tion for a specific dowr.wiud area and other, related Information. 

B. Formation orvl Ceoactry of Stem and Cloud 

After a surface burst, the fireball forms the shape of a sphere, 
and ita radii, both horizontal and vertical, expend exponentially with 
altitude as it rises. Since the expansion coefficients are different for 
the two radii, the rircball grows to t.he form of an oblate spheroid, vith 
a circular tor-view and an clllpticul side-view. When the fireball reaches 
its flnul stabilized height, It is commonly known as the "cloud”. A stem 
io shaped by the trajectory of a continuous!;/ expanding fireball and 
assumes the fora of an inverted exponential horn. 
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The Initial spherical radius of the fireball, R , the final 

s 

horizontal semi-axis, a, the final vertical semi-axis, b, and the final 
height of the center, h, of the cloud ■ ere related to the weapon yield, 
V, through empirical data as shown In the following sealing functions: 


2.09 x 10 2 W 0,333 ft. 

V 

• 1 to 10 5 KT 

(1) 

2.45 X 10 3 V 0 - 1 * 31 ft. 

V 

» 1 to 10 5 KT 

(2) 

1.40 x 10 3 w 0 * 300 ft. 

V 

- 1 to 10 5 KT 

(3) 

0.66 x io u v 0 " 1 * 1 * 5 ft. 

V 

- 1 to 28 KT 

<«0 

1.68 x 10 U V 0,164 ft. 

V 

- 28 to 10 5 KT 

(4a) 


At a given altitude, Z, the horizontal and vertical serai-axle 
of the fireball can be obtained by 


*z - V 


\ z 


\ - V 


V 


(5) 

( 6 ) 


where a , to , x , X. can be solved from the known information of R , a, b, h. 

O O ft 0 0 

The geometry of both stem end cloud is shown In Figure 1. 

C. Particle Size Paraj-eter and Ionization Rate Contours 

As mentioned previously under the basic assumptions, the fission 
products arc either fused Inside or adsorbed on the outside of the dust 
particles and are Inhaled into the fireball. The ultimate ionization rate • 
contours, or fallout patterns, on the ground are related to the size of the 
falling particles. The particle size parameter, a, is defined as the ratio 
of wind velocity, V , to the average falling velocity of the particle, 
since It la more concerned with the falling rate than the actual physical 


size of a particle 
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For the 9tea, It is assumed that the particles at a given altitude in 

the rising fireball are concentrated on a horizontal plane through the center 

of the fireball. Therefore, all particles which have equal values of a fall 

out of the stem starting at the same altitude. Any particle starts to fail 

when Its instantaneous falling rate, V^, equals the rate of rise of the our* 
d2 _ 

rounding air m03S, tt* There i3 an empirical relation bctueen V , V , and 

UU Ct l 

the flrebalL altitude Z: 

V 2 /V f - p ♦ qZ (7) 

vhere p • 0-95. q ■ 1.02 x 10*^, where Z ■ 5»C00 to 50,000- ft 

with particle diameter d ■ 200 to 1,200 microns 
p - 0.58, q - 1-7 1 * x 10* 5 , when Z - 50,000 to 110,000 It 
with particle diameter d - 300 to 1,000 microns 


There io also an empirical approximation for Z: 


•M 

Z 0 (1 - e ^ 


for t » 20 to •‘BO seconda (B) 


where Z Q lo a yield dependent multiplier 
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>2 io 0.011 sec . 


Itous, 


Z ' g - >Z <*o - 2 > 


(9) 


For values of Z 


see Figure 5- 
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Therefore, the particle site parameter, a, for the particle group 
falling free the altitude Z la 


\ Vjp ♦ q2) 

V f ' V 2 o - _Zl 


' ( 10 ) 


Moreover, it can be determined from geonetry of the stem that 


- X * a^ ( ♦ for a 


min* 


for a ) 
max 


( 11 ) 


Therefore, Q_and a cay be obtained for the particles falling 

QLELX mill 

at a distance X from ground zero. In addition, the rising and falling time 
for particles of a given value of a can. bo determined from combining Equa¬ 
tions (8) and (10), as shovn below 


* r (sec) 


1_ 

*2 


In ( 


■Vo * V*o 


j* 


t f (oec) 


°(°V 0 * V } 

V v («^ ♦ V y q7 


( 12 ) 

(13) 


The outer dimensions of the cloud are defined by 



(14) 


Consequently, and a n for the downwind distance X con be 

derived, as shovn by the equation 


Q 

m 


hX 


t y? 


LL 


Va 2 ) 


« 2 ** \ t . 

(a - V ) [h ~ b (1 - 




u 


h 2 - b 2 (l - y 


/a 2 ) 


( 15 ) 


This calculation usually includes the l8o sec. overage delay time due 
to fireball toroidal circulation, which should be subtracted in order 
to obtain the actual rising time. 



- 11 - 


(where the sign Is ♦ for a B&x > 611(1 * for a ni i n ) 

Since the cloud is assumed to he stabilized, all particles with 
the same values of a fall along the same slope. Hence, the tines of arrival 
and cessation for particles of a given a can be determined from the geometry 
of the cloud 


t a • « (h ♦ z a ) / V y 


(16) 


t c - « (h ♦ z c ) / \ 


(17) 


vhere z and z are the smallest and largest intercepts of the slope to 

a c 

the outer dimension of the cloud, respectively. 


QtX - Oh) b 2 ♦ ab J (a 2 * a 2 !? 2 ) (1 - y 2 /b 2 ) - fx - oil) 7 

« 2 . A 2 


Z m 


(19) 
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(where the sign is + for i , and - for z ). 

c a 

Assuming that the ionization rate, or fallout intensity, 1^, 
vnrLeii exponentially with downwind distcnee, end its dlstrltortion is closely 
related to a and W, also confirmed by cr.plrleal test data, a complete set 
of scaling functions which defines the voriatioc of intensity versus down¬ 
wind distance and the bi-elliptical outline of the ionization rate contours 
is tabulated In Chapter 3 of "Fallout and Radiological Countermeasures"^^. 
Those scaling functions arc not repeated in this report, but a corresponding 
computer program was established, as presented in Figure 5. A typical 
ionization rate contour for a 10MT weapon yield and a 15 mph wind velocity 
is shown in Figure 6. Since the only inputs are weapon yield and wind 
velocity, this computer program is relatively simple to understand and 
to use by others to obtain the information on ionization rate contours. 

It should be pointed out that these contours have been corrected 
to a standard reference time, H ♦ 1 hour, by a decay correction factor 


I (t: actual time) 
d(t, l) » j x ( 1: « ♦ 1 hour) 


(19) 


This decay correction factor can be obtained directly from a 
typical close-in fallout decay curve ns shown in Figure 2 , 

Knowing the ionization rate, the exposure dose, D^, in the cloud 
fallout area, which is of greater Interest than the stem fallout a*-ea, can 
te estimated by 


D K */ t ¥ t > dt 


( 20 ) 



(u os/noissu a3d bH a 5( _oi > 3iva NonvziNOt 












D. Condensation of Fallout or.d Practlonatlon tfumbera 


During a nuclear detonation, a tremendous amount of heat Is created 
vy nuclear fission and fusion. This heat melts and also vaporises every¬ 
thing inhaled Into the air mass fron the ground and forma the familiar 
fireball. When the fireball rises, Its volume begins to expand, and lta 
temperature io assumed to drop according to the Ideal gos lav. The material 
swept Inside the fireball starts to liquldlfy and finally solidify to 
particles, in or on which the fission products begin to condense. The 
condensation process Is generally divided Into two time periods as mentioned 
In the basic assumptions. The dividing temperature depends largely on the 
composition of the fireball. For a model surface burst on a soil consisting 
of silicate minerals this temperature Is assumed to be l6?3°K, as cost ele¬ 
ments will have started to solidify e.t this temperature. Since the two 
periods are charactarlzed by the coexistence of either gas and liquid, or 
gas and solid, the vaporization or sublimation pressure plays an important 
role In the condensation process. The following formula for the evalua¬ 
tion of the pressure la presented In "Fallout and Radiological Counter- 
measures ' : 


log Pj - j ♦ B ♦ G (21) 

where Pj Is the vaporization or sublimation pressure for element J, 
gm/sq. cm., 

T lo temperature of the fireball, °K, 

A, B, C are the empirical constants for vaporization or sublimation 

reactions of elements contained in the fireball which were 

(4) 

tabulated In Vdller'a manuscript . 

. The established computer program for estimating Pj Is shown In 
Figure 3» 





Figure 3. Cooputer Program for Estimating Sublimation Pressure 
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Figure 3 (Continued). Computer Program for Estimating Sublimation Pressure 
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Figure 3 (Continued). Computer Program for Estimating Sublimation Pressure 
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Wgure 3 (Continued). Omputer Program for Katl mating Subllmticn Pres*ur» 
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figure 3 (Osutlnued). Cbaputer Program for Estimating Sublimation Pro*sure 
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The elements In the firebaiL are radioactive ond hence the number 
of moles of each species is constantly changing with time. If y,,(t) 
denotes the amount of radionuclides of element J ar.d mass number A present 
at time t after fission, then the total amount of element J present at tine t 


'fj(t) « atoris or ®oles 


Correspondingly, the total amount for the entire chain yield of 
mass number A Is 


Y JA (t) atoms or moles 


in which 16 constant, except for masu chains containing neutron emitters. 

For the first period of condencation, the material balance for 
element J between gas and liquid phase la 

ty° ■ Vj* (,) • V;a°< 1 > * Vja<‘> (2 ‘ 

where nj A °(t) la the number of moles of element j with mass number A, 
which is mixed with moles from other mass chains to 
fora n moles of vapor, and iu abbreviated to n®, and 
r. .(t) is the number of moico or eicment J with msec number A 

JA 

dissolved in the n(/) moles of liquid carrier, which 
Is the particle in the liquid phase, prior to solidi¬ 
fication, and is abbreviated to n.. 

J 

For a p.-rftct gas, the ratio between the mole fraction of clement J 
in the vupor pnasc and that in the ,’iquld phase is 
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N? n®/n n® KT/p V k 

Nj " nj/n(/J * * p^ (25) 


where 



n 


n(i) 


R 


T 


P J 

V 


is the cole fraction of element J in the vapor phase 

is the cole fraction of element J in the liquid phase 

is the number of moles of vapor 

is the total moles of liquid carrier 

Is the nolar gas constant 

is the absolute temperature 

is the total vaporization pressure 

is tiie molar volume, and 

is the Henry's Law constant 


From Equation (25), the following relation between n® and n can oe 

J J 


obtained: 


n.k. 


NO/ 


VI RT 


"j k j 


(26) 


where k° » k^/[n(/)/vj RT 


Substituting Equation (26) into Equation (24) yields 

■ <‘ * 9 Vj*<'> 


(27) 


The fractlor.atlon number of the first period of condensation, 

r (A,t), is defined to be the fraction of the moos chain condensed, or in 
o 

mathematical fora 


r Q (A, t) - ^ Ej n JA (0 



y .1A (t) 
l ♦ k? 


(28) 








vh* re 



where t, in seconds, is related to the fireball temperature, T, by 

T ■ U.66 x 10 3 w’°*° l ° exp(>0.546 «*°* 373 t)°K (33) 


and the data of yj A (t) ve rc calculated and tabulated as Hj(A, t) by holies 
and Ballou^ . ;ter:cc, the fractionation number, r o (A, t), can be evaluate i,r 
the elements contained in the fireball. 

During the second period of condensation, if an exceed of solid 
surface area iu present, the reiml condensed at r.ny time after colldlflea- 
tior. of the carrier, io given by 
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a i “ a i ’ °j 


( 5«0 


where o.' Is the amount of element J condensed on the surface of 

the solid particles 

n? is the amount of element J that has not condensed in the 

V 

liquid carrier, and 

n^j is the amount of element J In the vapor phase. 


For a perfect gas 


n i ‘ iar 


( 35 ) 


where 


13 the sublimation pressure and may be evaluated according 
to Equation ( 2 l). 


The material balance for element J is 


7 j(t) • Dj(t) ♦ nj(t) ♦ n”(t) 


(36) 


Therefore, employing Equations ( 27 ) and ( 35 ) 

. £ 

RT 




( 37 ) 


or 


VjA 


1 ♦ K 


Z A y JA ’ K? E A P «A 


( 38 ) 


In thla equation the partial pressure of each nuclide is propor¬ 
tional to Its abundance at uny given time, ar.d is given by 


■ < x 


(39) 
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‘ (I ’ r o )y JA* 

v O - o 
■J * Z k y ik 


The f rac t i on*. 11 on number of the second period of condensation 
r^(A, t) Is defined to be the fraction of the mass chain condensed up to 
a given time during the s-.cond period, or expressed mathematically 


r; (A, t) 


h 


V r a i 
0.23 SVRT J P JA J 


( 40 ) 


(A t\ - \ „ ^ / A f \ m - - - 

‘o lA ’ } 1 O’’ z; O^jtWK? 


£ / 3A/y >°A 

E J y jA 


( 41 ) 


vhere 


is the fireball volume and is given by scaling the function 
V - 7.62 x 10 11 W exp (0.510 W*°' 373 t) cm J 


B le the ratio of fl66lon to total yield, and 

0.23 Is a conversion constant changing weapon yield, W, to 

r.urter of fission-moles per atom. 

Physically In Equation (-1), p 1 ? V/RT glveo the number of molts of clement J 

J 

In the vapor phase sr.d 0.23 Wl? gives the total n'lmter of moles of the 

V 

element produced nr.d In existence at a given time In the second period. 

The combined tern in Elation (4l) gives the fraction still remaining In 
the vapor -tats during the second period of condensation. The 1 - r Q (A, t) 
torn glv< j the fraction which has rot. condensed in the liquid carriers during 


their existence In the 


rot period of condensation. 
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The fractionation numbers yield the information about the fractions 
of radioactivity condensed in, or on the particles, at any tine after a burst. 
Therefore, the fractionation numbers not only make the conversion between 
Intensity and activity possible, but also lead the vay to the evaluation of 
soluble and Insoluble activity vhich is important to biological availability 
and radiological counte measures research. 

S. Conversion between Intensity and Activity 

The conversion between Intensity (measured as ionization rate) and 
activity density in fissions per unit area is presented in "Fallout and 
Radiological Countermeasures",^**) based on the simple assumption that 


I x (t) 


I x (t)(ideal) 

ideal) 


yt) 


o*> 


where the conversior factor, K^t), * s further broken down into the fol¬ 
lowing components 

K^t) - Dq x (r x (t) i fp (t) ♦ l 1 (t)l 


(>3) 


in which 


r x (0 




i t (t) 


is an instrument response factor, usually assigned a 
value of 0.75 

is the terrain shielding factor, usually assigned a 
value of 0.75 

l8 the gross fission product fraction number, sometiaea 
aloo designated as r Q (t) 

is the air ionization rate per fission at 3 ft above 
an Infinite, ideal plane for a uniform distribution or 
the normal fission product mixture, and 
Is the cone unit for neutron Induced activities. 
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Details about the neutron Induced ionization rate, i^(t), are almost entirely 
classified, but generally assigned^ as 1,9$ of i^ p (t). The other components 
In equation (^ 3 ) can be readily derived. 

• Detailed information about the activities of the U-235 fission 

( 7 ) 

products has been provided by Bolles and Ballou . These data can be 

( 8 ) 

correlated to the data supplied by Katcoff' and applied to other types 
of fission elements, such as U- 238 , Pu-239 by taking proportions of the 
corresponding mass number. This conversion has been systematized and 
tabulated in "Fallout and Radiological Countermeasures". Moreover, Hiller ^ 14 ^ 
translated the activity data from Bolles and Ballou into Ionization rate by 
a disintegration multiplier, m. Both tables of correlation and of disinte¬ 
gration multipliers are reproduced from Miller's manuscript as Tables I 
and II. Therefore, the ideal or normal ionization rate of fission product 
mixture may be computed according to the following formula: 

i fp (t) - BnA t (kk) 

4 

where A fc is the activity of the nuclide at time, t, per 10 fission. 

In accordance with the earlier discussion, the fission products 
do condense and fractionate with time instead of condensing completely at 
the very beginning as in the ideal situation. Therefore, the real ioniza¬ 
tion rate, i f *(t), should be multiplied individually for each radionuclide 
by its fractionation number at any specified time. Mathematically, this 
may be expressed by 

i f *(t) - Bn A t (r Q (A, t) ♦ r^(A, t)J (Uj) 

The (A, t) term vanishes in the first period of condensation. 

As time Increases, the sum of the two fractionation numbers approaches unity. 


i 
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TABUS I 

(4) 

Cumulative Mass-Chain Yields of Fission Products 
(Values are la per cent of fissions} 


Mass 

U 2 ^ 

U 2 ^ 

Pu 2 ^ 

number 

Thermal 

Neutrons* 

Fission 

Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

Neutrons 

Pi 33 Ion 
Neutrons 

72 

1.6xlO* 5 

4.6xl0’ U 

5.0xl0" 6 

m 

•4* 

1.2x10 

• 

73 

-4 

l.LxlO 

0.0012 

3.7xl0” 5 

« 

2.2xl0* U 

- 

74 

(3-2xlO*V 

0.0034 

1.1x10 

0.001 

4.1xlO~ 4 

0.0011 

75 

(8.8x10^) 

0.0062 

8.3x10*^ 

0.0040 

J4- 

7.xl0 

0.0023 

76 

( 0 . 0029 ) 

0.012 

0.0012 

0.0078 

0.0014 

0.0051 

77 

O.OO 83 

0.023 

0.0038* 

0.014 

0.0026 

0.011 

78 

C .021 

0.048 

0.0095 

0.026 

0.0049 

0.025 

79 

(0.041) 

0.096 

0.019 

0.053 

0.0090 

0.043 

80 

(0.077) 

0.19 

0,045 

0.096 

0.016 

0.075 

8 l 

0.14 

0.21 

0.088 

0.18 

0.030 

0.14 

82 

(0.29) 

0.50 

0.20 

0.35 

0.056 

0.23 

83 

0.544 

0.80 

0.40* 

0.66 

0.10 

0.37 

84 

1.00 

1.3 

0.85 

1.02 

0.17 

0.60 

85 

1.30 

1.85 

0.80 

1.45 

0.28 

0.92 

86 

2.02 

2.5 

b+ 

WJ 

CD 

+ 

1.9 

0.45 

1.15 

87 

(2.94) 

3-3 

l.S» 

ro 

0.73 

1-5 

88 

(3.92) 

4.2 

2.45 

2.7 

1.2 

1-9 

89 

4.79 

5.1 

2 . 9 * 

3.17 

1.9* 

£ . ^ 


Continued 



TA3LE I (continued) 


C-JMulatlve Mass-Chain Yields of Fission Products 
(Values are in per cent of fissions) 


y235 -j23S 



Thermal 

lieutrons* 

Flsslor. 

2,’eutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

Nc-utrcas 

Fission 
Neutrons 

20 

5.77 

5.8 

3.2* 

3.7 

2.4 

3-0 

21 

5-e\ 

3.85 

3.6 

4.3 

3-0 

3-7 

9£ 

6.03 

6.0 

4.1 

4.8 

3.7 

4.4 

93 

6.45 

6.4 

4 >85 

5.2 

4.6 

5.0 

94 

6.40 

6.4 

5.3 

5-45 

5.5 

5A 

93 

6.27 

6.3 

5-7* 

5-6 

5.9* 

5.6 

9^ 

6.33 

6.3 

5.8 

5-7 

5-7 

5.3 

97 

6.0? 

6.1 

5.7 

5.64 

5-6* 

5.2* 

9* 

5.78 

5.8 

5-7 

5.6 

5.4 

5.4 

92 

6.06 

6.1*-* 

6 . 3 * 

6 . 2 *» 

5*9* 

5-9* 

100 

6.30 

6.7 

6.1 

6.4 

6.0 

6.4 

101 

5.0 

5-3 

5.5 

6.5 

6.0 

5-9 

102 

4.1 

2.9 

5.6 

5-9 

5.9 

5.3 

103 

3.0 

1.7 

6.6 

5.0 

5 . 8 * 

4.6 

104 

1.8 

0.95 

5.4 

3.2 

5.0 

3.5 

102 

0.90 

0.54 

3.9 

2.2 

3.9* 

3.2 

106 

0. j3 

O.jO 

2,7* 

1.5 

5 . 0 * 

3.6 

107 

0.19 

0.17 

i.35 ' 

1.0 

4.0 

3.1 

103 

(0.035) 

O .095 

0.67 

0.70 

3.0 

2.6 

109 

(0.0;9) 

0.0>3»** 

0 . 32 * 

0.43 

1.5* 

1.9* 

110 

( 0 . 00 ) 

0.030 

0.15 

0.33 

0.65 

0.81 

111 

(0.015) 

0.0..2*** 

0.0/3* 

0.23'** 

0.27* 

O .34 

li: 

(0.013) 

0,0_0*»» 

0.046* 

0.19 

0.10* 

0.14* 

113 

(0.0L?) 

0.018 

0.043 

0.17 

0.055 

O.OpO 
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TABLE i (continued) 

Cumulative Mass-Chain Yields of Fission Products 




(Values are 

in per cent of fissions) 


Mass 

Number 

U 2 * 5 

U 2 ^ 

Pu 2 *> 

Thermal^ 

Neutrons 

Fission 

Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thi-rrrifll 

Neutrons 

Ficslon 

Ne-tron; 

114 

(O.Oll) 

0.017 

0.041 

0.16 

0.046 

0.075 

U5 

0.0104 

0.017"** 

0.040* 

0.15*»* 

0.041* 

0.069* 

116 

(o.oio) b 

0.017 b 

0.039 

0.14 

0.039 

O.O 65 

117 

( 0 . 010 ) 

0.017 

0.039 

Q.l4 b 

O.O 38 

O.C 65 

118 

( 0 . 010 ) 

0.017 

0.4p b 

0.14 

O.O^ 

0.064 b 

119 

(o.ou) 

0.017 

0.041 

0.14 

0.039 

0.064 

120 

(o.ou) 

0.018 

0.042 

0.15 

0.041 

0.065 

121 

( 0 . 012 ) 

0.020 

0.044 

0.16 

0.044* 

0.066 

122 

(0.013) 

0.022 

0.046 

0.17 

0.047 

0.069 

123 

(0.015) 

0.030 

0.050 

0.19 

0.052 

0.076 

124 

(0.017) 

0.053 

0.055 

O .23 

O.O 58 

0.082 

125 

0.021 

0.095 

0.072 

0.33 

0.072* 

0.14 

126 

(0.053) 

0.17 

0.175 

0.48 

0.175 

0.35 

127 

(0.H5) 

0.30 

0.39 

0.70 

0.39* 

0.80 

128 

0.37 

0.54 

0.77 

1.0 

0.77 

1.9 

129 

0,90 

0.95 

1.45 

1.5 

1.45 

2.5 

130 

2.0 

1.7 

2.5 

2.2 

2.5 

3-2 

131 

(2.63) 

2.9 

3.2* 

3.2 

3.8* 

3-8 

132 

(4.31) 

4.3 

4.7* 

4.4 

5-0 

4.6 

133 

(6.43) 

6.1 

5.5* 

5.4 

5-27* 

4.9 

134 

( 7 . 80 ) 

7.3 

6 . 6 * 

6.5 

5-69" 

5.2 

135 

(6.40) 

6.3 

6 . 0 * 

5.9 

5.53* 

5.1 

136 

(6,36) 

6.4 

5-9* 

5.3 

5 . 06 " 

5-3 

137 

(6.05) 

6.0 

6.2 

5.85 

5-24' 

6.4* 


Continued 
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IA3LE I {continued) 

Cumulative Mass-Chain Yields of Fission Products 




(Values are 

In per ednt 

of flpelons) 



Mass 

1 'usher 

,235 

1 

238 

Pu 239 

'Ihemal 

Neutrons 

Fission 

t’eutro.-.L= 

Fission 

Neutrons 

8 -Mev 

Neutrons 

Chennai 

Neutrons 

Fission 

Neutrons 

133 

5-7** 

5-7 

6.4 

5.9 

5-5 

5.4 

139 

(6.34) 

6.4 

6.5 

6.0 

5-7* 

5-2 

l4o 

6.44 

6.4 

5-7* 

5-6 

5 . 68 * 

5.o« 

14 i 

(6.30) 

6.3 

5-7 

5-5 

5*2* 

4-7 

142 

(5 85) 

5-9 

5*7 

5.4 

6 . 69 * 

4.9 

1^*3 

(5-87) 

5.8 

5-5 

4.97 

;.4* 

5.0 

ill* 

5.67 

5.1** 

4.9* 

4.3“ 

5-29* 

4.6 

145 

3-95 

4 .2 

3.7 

3-7 

4.24* 

4.4 

1 ! *6 

3-07 

3.3 

3.1 

3.17 

3-53* 

3-7 

l4? 

2.38 

2.5*# 

2.6“ 

2 . 7 ** 

2 . 92 * 

3.0 

l'i 8 

I .70 

1.85 

2.0 

2.27 

2 . 28 * 

2.36 

H9 

1.13 

1.3** 

1.45 

1.9“ 

1*75 

■ 1.86 

150 

0.67 

0.80 

1.05 

1.45 

1.38* 

1.48 

151 

0,45 

0.50 

0.74 

1.02 

1.08 

1.16 

132 

0.265 

0.31 

0.50 

0.66 

0.83» 

O .92 

153 

0,15 

0 . 19 ** 

0.32 

0.41“ 

0 . 5 ? 

0.60 

151 

0.077 

0 . 0-56 

0.19 

0.25 

O. 32 * 

0.37 

155 

0.033 

0.048 

o.ll 

0.15 

0.20 

O .23 

15' 

0.014 

0.023“ 

0.0 66* 

0.092-*' 

0.12* 

0.14 

l? ( 

0.0078 

0.012 

O.034 

O.C57 

0.064 

0.075 

I r »- 

0.002 

C.9062 

0.016 

0.032 

0.034 

0.043 

131 

0.00107 

0 . 00 }-" 

O.OC 90 ** 

0.017“ 

0 . 020 ““ 0.025 


Con timed 
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TABLE Z (continued) 

Cumulative Mass-Chain Yields of Fission Products 
(Values are In per cent of fissions) 


Mass 

u 235 

^38 

Pu J 

39 

Number 

Thermal 

Neutrons* 

Fission 

Neutrons 

Fission 

Neutrons 

8 -Mev 

Neutrons 

Thermal 

Neutrons 

Fission 

Neutrons 

160 

3-5x10-4 

0.0012 

0.0036 

0.0085 

0.0092 

0.011 

l6l 

7.6x10-5 

4.6x10-4** 

9-4x10-4 

0.0044** 

0 . 0038 *^»* 

O.OOJi 


■•Seymour Katcoff, Fission-Product Yields From U, Th and Pu , Nucleonics , 

Vol. 16, No. 4, p 78-85 U95B) 

r. Bunney, E. M. Scadden, J. 0. Abriam and N. E. Ballou, Radiochemical 
Studies of the Fast Neutron Fission of L'-35 and U^38 t Second Inter¬ 
national Conference on the Peaceful Uses of Atonic Energy, A/Conf. 15/F/643* 
USA, June 1958 

«**G. P. Ford, J. S. Gilmore, et si, Plsclon Yields , LADC- 3083 , 1958 

*»«L. R. Bunney, E. M. Scadden, J. 0. Abrlam, N. E. Ballou, PI scion Yields 
In Neutron Fission of Pu^^° , USNTOL-TR-268, 1958, Unci. 

a. Parentheses indicate estimated values or where Katcoff 8 value was altered 
in order to adjust the yields to a gross sum of 100 in each peak. 

b. Line Indicates division of two peaks that woo used for Individual peak 
sumo. 





Nuclldr 


Zn-72 

2n-74 

Ca-72 

Ca -73 

Ca-74 

Ce-75 

0-77 

Ce-78 

Aa-77 
As -78 
Aa-79 

Se.-8ljB 

Sa*-8i 

Se -83 

Br -83 

Br-84 

Kr-833! 

Kr.-85« 

Kr^-85 

Kr -87 

Kr -88 

Rb-88 

Rb -89 

Rb-91 

8r-89 

8r-90 

Sr-9l 

8r-92 

3r-93 

V -90 

Y,-91n 

Yp-91 

Y-92 

Y-93 

Y-9i* 

Zr-95 

2r-97 


•S. 

TABLE 11 


Air Ionization Disintegration ftaltlvlien tor 
-na_ Fisa ion Product and oTh'er RaJlon^n?^^) 


n 

10* 9 7/hr-ft 2 

v^To/SecT”’ 


ITuclide 


€.75 

4,22 

14.3 

2.11 

12.4 

0.216 

13.1 

2.74 

0.0735 

2.51 

0.0 

0.LL7 

0.0 

12.8 

0.0174 

8.12 

0.0959 

0.973 

0.0316 

7.03 

7.40 

4.30 
12.0 
6.6 7 

0.0 

0.0 

5.03 

6.93 

1.20 

0.008 

3-24 

0.0135 

6.56 

0.817 

6.67 

fc-32 

0.346 


JTb.-95« 

Nbi*95 
Kb,-97 b 
Nbp-97 
Nb-98 

Mo-99 

Mo-101 

Mo-102 

Tc-99b 
T c-101 
Tc-102 

Bu-103 

Pu-105 

Ru-106 

Rh.l03« 

Rh.-105m 

Rhi-105 

Bh -106 

Rh-107 

W -109 

Pd-Ul 

Pd-112 

Ag-109* 

Ag-lii 

Ag-112 

Afi-113 

Ag-115 

Cd -115 

Cdf-115 

Cd>U7m 

Cd-118 

Cd-120 

In-115 

In-117 

In-ll8 

In-119 

Sn-121 
Sn-123 
Sn-125 
Gn -126 
Sr.-127 


a 

10' 9 r/hr-ft 2 
(dis/sec) 


1.41 
4.56 

4.42 
4.02 

11.7 

0.772 

6.85 

0.0 


0.782 

2.03 

2.45 


2.99 

4.36 

0.0 


0.150 

0.283 

0.559 

1.8l 

5-08 


0.0 

5-00 

0.0225 

0.172 

0.175 

*♦.59 

O .158 

1.28 

O.I 87 
1.11 
8 .H 
0.0 
3.08 


1.12 

1.36 

13.5 

0.0371 


0.0 

0.0 

6.19 

0.0 

3.08 
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Table 11 (continued) 


Air Ior.l2atlon Disintegration Multipliers for 


Ruellde 

the Fission Product and 

Other R'idlonuclldea 

a 

10* 9 r/hr-ft 2 

a 

20’ 9 r/hr-ft 2 

rtuellde 


(dlt/see) 


(dls/aec) 

Sb-125 

2.78 

Ce-l4l 

0.470 

Sb-126 

16.1 

Ce.l43 

2-34 

Sb-127 

2.70 

Ce*l44 

0.188 

Sb-i28 

13.0 

Ce-145 

4.22 

Sb-129 

6.10 

Ce-146 

1.60 

Sb-131 

3.66 




Pr-l43 

0.0 

Te. -125** 

0.0525 

Pr-144 

0.163 

Te,-127 

0.181 

Pr-145 

0.0 

Tei-127 

0.0249 

Pr-146 

6.42 

Tcf-129a 

0.398 

Rd-l47 

O .968 

Te^-129 

1.52 

Te,-131 

8 .U2 

Nd-149 

2.31 

Tei-131 

2.08 

IM-151 

6.44 

Te -132 

1.52 

Ra-l47 


Te.-133a 

0.790 

0.0 

Tei-133 

9.12 

P»-l49 

7.25 

Te-13l* 

7.89 

Pb-150 

4.99 


Pb-151 

2.21 

1-131 

2.3S 

Pm-152 

3-79 

1-132 

13.0 

Pm-153 

4.76 

1-133 

3.59 


0.0544 

1-134 

9.63 

Sou 151 

1-135 

9.60 

Em-153 

0.455 


Sm-155 

1.80 

Xe -131a 

0.207 

Sm -156 

1.20 

Xe*-133« 

0.1*13 

Sm-158 

3.08 

XeJ-133 

0.304 


0.314 

Xe^-135. 

2.69 

Eu-155 

XOp-135 

1.55 

Bu -156 

5-55 

Xe -138 

7.09 

Eu-157 

3-78 


Eu-158 

7.21 

Ca-137 

0.0 


0.430 

Co-130 

U .2 

Od-159 

Co-139 

4.78 

Tb-l6l 



0.123 

Ba-137« 

3.65 



Ba-139 

0.888 

tl-237 

1.13 

Ba-l40 

1.10 

U-239 

0.338 

Ba-lUl 

4.22 

U-040 

0.0 

Ba-l42 

5 .Y 6 




:ip-239 

1.11 

U-lUO 

13.0 

Kp-240 

2.07 

Lo-l^l 

0.400 

Kn -56 

9.40 

La-l42 

11.6 

La-ll*3 

6.67 
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The gross Ionization fractionation neater, r^(l), is defined as 
the ratio between the ionization rates at H + 1 hour of condensed and normal 
fission products. It nay be expressed as 

r x (l) - 


*f^> 


Thus, all the components for the conversion factor, K^t), have teen cited 
and derived. The intensity-activity conversion, shown as Equation (h2), 
is very useful in fallout contamination studies and the determination of 
other properties of fallout. 

F. Biological Availability and Atom-Concentration Intensity Ratio 

The fraction of radionuclides condensed on the outside of fallout 
particles Is potentially available for biological uptake. By an analogy 
to the Intensity-activity conversion presented in the foregoing section, 
the ratio of the number of atoms condensed on the exterior of the 
particle to the standard Intensity at H ♦ 1 hour is given toy the following 
formula 


N a <A) 

rnr 


Y A r;(A, t) 




(**7) 


where 


*& A) 


Is the number of atoms (of the end member of the mass 
chain and condensed on the outside of the particle) 
which land per square foot of ground. 

The other terms in Equation (^7) are the some as defined previously, 
ihe ratio Is colled the atom-concentration intensity ratio. 

Ac shown In Equation (12) the fireball rising time, t, is a function 
of the partlcle-slzo part meter, a. Because of tine dependence this ratio is 
also dependent on Q. However, the Or can be approximated and simplified by an 


average value o^, which is 
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X 

h 


Since ci^ does not vary with crocs-wlnd dis*ft r,ce > atom-concentration 


intensity ratio holds constant &s a first approxinatlth* f° r ®H valaco or" cross- 

Moreover, ns shall 


wind distance at a given value of downwind di6tnr.ce X 
be seen in the following section, a given value of ci yepresento a group of 
particles with a snail range in diameters. Therefore; f hc ratio N^/l ia essen¬ 
tially associated with a given particle size. 


The values of N a (A)/l(t) for six biologically important icotopeo at 


different downvind distances for a typical 10 MT burst have been calculated 
as shown in Table VI. 

This ratio can be applied to water contaxnllth t l° n studies to yield 
the activity of soluble radionuclides in water cupplie 8 > W t q(^> O 1° 
replaced by r Q (A, t), the ratio W Q fl will provide the Information for the actlvltyj 
of insoluble radionuclides, since those radionuclides fased into the silicate 
particles are assumed to be not readily soluble in wat®r* 

0. Ionization Rate Contour Ratios and Particle Size £>< strlbutlon 

Several ionization rate cootour ratios havs l ecu defined to deter¬ 
mine special properties of fallout. They aro: 

(a) The mass contour ratio, M r (t), which Id tho ratio of the 
fallout mass per unit area at any downwind distance 10 the ideal ionization 
rate at that location, or mathematically: 


M r (t) 


where 


Is the mass of fallout per unit area ft! tu V downwind 
distance, X 

I^t) is the ideal ionization-rate at 3 ft wLove an extended open 
area covered with fallout. 


(**9) 
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The mess of fallout, m, empirically determined to be a function 
of both a and V, la 



a - f(t») W °‘ 9L7 Rg 



(50) 

where 

f(a) - 5.33 X 10 * 12 u’ l ‘ 25 mg/fioslon. 

a 

■ 0.1 to 0.9 



- 5.64 x 10* 12 mg/floslor.. 

a 

- 0.9 to 2.0 



* 7.14 x 10 * 13 H^/flesion 

a 

>2.0 



According to Equation {’ 1^)1 lb * 8 found that 
yi) * Area - D(l)q x U a (l ' * M 01 W 


Therefore, the mass cohU™** ratio associated with H ♦ 1 hour, M r (l), 
for a ground surface burst will to# 


M r (l) 


f(a) W ^’ 083 _ _ mg/sq ft 

D^l) qj^Tr^Cl) l fp u) * l^ljj r/hr at 1 hr 


(52) 


(b) The fisulon-produOl contour ratio, FP r (t), which is the ratio 
of the number of atoms, or moles, of fission products per unit area at any 
do.wlnd distance to the Ideal iohl^Rllon rate at that location, or matho•• 
natlcally; 

PP r (t) * *rjh lt) (53) 

where R 1<J the numbef >'f atoms, or moles, of fission products 

fp 

per unit areft 1 

From the dlscuoelon of Ruction 6., this contour ratio, corrected 


to B ♦ 1 hour, for a ground surfsl lu burst may be expressed as 


-37' 


P? r (l) 


1-16 x 10 ~~* r u (r o (A) * r o {A)1 -A 
" Diir^ir Q (i) i.-jrijvi^ijr 


taoles fp/sq ft / c i •. 

r/hr at I hr' (5 ^ 


(c) CTe fraction of device contour ratio, FD^ft), which Is the 
ratio of the fraction of veapon device per unit orea at any downwind distance 
to the Ideal ionization rate at that location. This ratio when corrected to 
H 1 hour may be expressed as 

.o), 

_^1 x 10 

■*X c * t: u-~' *i 


™r {l) “ Pd T q, cv"wi) T ; -r i)% ' 1(0 1 r/hr/3 « ft at 1 br - 


(53) 


These conto-r ration predict various characteristics of local 
fallout, as they yield information about the mass, the number of atoms of 
fission products per ’-tit area, the fraction of weapon yield per or.lt Area 
and other properties of fallout at a specific location. This information 
is valuable for the design of fallout shelters and other radiological 
countermeasures. More-over, the use of these contour ratios ur.i their 
scaling functions maker possible the extrapolation of a limited amount 
of experimental data to cover a number of operational cases since the con¬ 
tour ratios are not constant but are point functions whose values depend 
oa many variables. Furthermore, the average density of the fission products 
can be estimated from the ratio of FP r (l) to M r (l). 

Another advantage or the C.C.D. fallout model ic that particle 
site groups for any giver, downwind location can be estimated by a graphical 
method, provided that the particle fall-rate is known. Since the fall-rate 
for a given else parameter a does denote a group of particles with a range 
of diameters, with the fell-rate for a given particle size parameter known, 
the range of diameter; for a group of ps:*iclcs can be determined, depending 
on the thickness of the cloud and the altitude from which the group falls. 


A family of curves which shows the size of ideal, spherical particles falling 
through a standard atmosphere from various altitudes as i function of fall- 
rate was prepared by D. E. CLark^^ and appears as a sketch in Pigure 8* 

When estimating the particle size range at a downwind distance, it 
should be first decided whether the source of the particle group originates 
from the stem or the cloud. The origin of the particle group can be generally 
Judged by comparing the downwind distance with the downwind range of the 
stem fallout. In case the distance falls within this range, the particle 
group predominant.ly originates from the stem; otherwise it is essentially 
derived from the cloud. Within the stem range, the maximum and minimum a 
vnlucs may be calcul-ied by Equation (11), which is r. laborious procedure. 

A computer program, cs shown in Figure 7, has been established to solve this 
equation. The height of fall for the stem fallout is derived from Equation (M) 
which is 


Z 


{Z 0 QX Z * V> 

°?8 * V * 


(56) 


For particle groups falling inside the cloud range at nny down¬ 
wind location (X, V), the crossvind distance Y should first be converted 

to the cloud coordinate, y, according to 

* 

y - a (57) 


where Y is the raxlri'in half-width of the croscwir.d distance on a 

m 

1 r/hr contour, and 
a is the cloud radius. 

Following this coordinate transfer, use equations (l$) and (l3) to 
Ktemine Q and the two extreme heights of fall, z and ; . For the actual 

ED 14 C 





distance of fall, the height of the cloud, h, should be added. The fall-rate 
is evaluated from Equation (10), which i6 the defining equation of a and 
which yields 



Therefore, with weapon yield and downwind distance known, and with 
the help of the fall-rate curves, the approximate particle size distribution 
at (X, Y) can be estimated. This information on particle size distribution 
is useful for studies of transport phenomena in fallout investigations. 

H. Conclusion - The Reality and Applications of the O.C.D. Model 

The O.C.D. fallout model mainly uses the technique of the scaling 
method, hence most of its formulae are derived from realistic experimental 
data. Therefore it serves as a very good bridge between the ideal theoretical 
hypothesis and practical empirical data. Moreover, this model has been 
applied to several real test shots, and was found to yield very satisfactory 
results. On the other hand, the fundamental thermodynamic theory on which 
this model is based is quite simple and convincing, and all the computations 
Involved are relatively easy. Another advantage of this model is its vide 
application, since it gives not only information on fallout patterns tut 
also that of activity, hence it is very valuable to radiological research. 
Although there arc still a few defects, such as constant wind velocity, 
the sole dependence on weapon yield, and the cigar-chape contour assump¬ 
tions, etc,, this model is generally considered very realistic and applicable 
to the problems of fallout research. Its main applications will be radio¬ 
activity prediction, countermeasure design and biological protection. 




Ja¬ 

in. COMPUTATIONAL METHODS OF THE FALLOW MODEL 

A. location Rate Contour* 

Although the physical characteristics and the mathematical approxi¬ 
mation for the rising stem and the cloud have been discussed In detail, 
little has been said about how the fallout contours are determined. Even 
though the available data from actual land surface shots Is small ar.d limited 
primarily to low weapon yields, some general conclusions have been drawn 
from the teat data. 

Because most of the observed properties of the fallout patterns 
are In terms of Intensity, the fallout contours are usually given In units 
of roectgens/hour (where the time was corrected to a standard time of 
H ♦ 1 hour). The general shape of the fallout pattern can be approximated 
by the overlapping of ellipses, for both the stem and cloud fallout. 

The dimensions necessary for constructing the ellipses are the 
distances from ground zero along the downwind axis where the lateral dis¬ 
placement la zero, (Y « 0). First, to construct the intensity vs. distance 
curve (see Figure **) the downwind distances of interest are the highest 
intensity and a preselected low intensity, which varies with wind velocity, 
but for a 15 mph wind la 1 r/hr at 1 hr. 

For the stem these locations have been labelled as follows: 

X^ and are the upwind and downwind distances, respectively, from 

ground zero to the 1 r/hr at 1 hr Intensity. 

Xg and are the distances to the upwind and downwind shoulders, 

respectively, of the intensity ridge. 

For the cloud the following notation is used; 

Xj and Xj are the arrival location for the group of particles of the 
the selected low intensity from the rear and front of the 
cloud, respectively. 





DOWNMNO DISTANCE X (MU.ES) 








X, is the distance from ground zero to an intermediate point 

o 

between Xj and X^. 

X^ la the downwind. distance from ground zero to the peak of 

the Intensity contour. 

Xg is the downwind distance from ground zero where the contour 

Intervals has a maximum half-width. 

The values for these quantities are determined by equations 

appearing In the Ionization or Intensity rate computer program. Eecause 

of the limited characters available on the computer printer, the notation 

used In the analysis of the model had to be changed. The corresponding 

notation between the model and the computer program is given in Table III. 

The following discussion describes the method for calculating the 

X *3 and how the scaling functions can be found in the computer program in 

Figure 5 . First, yield-dependent scaling functions for the size parameters 

a ■ a, y o^, a^, , . . , have been determined from data that has an 

effective, or average wind epeed of 15 mph. Thece ocaling functions are 
V V 

then multiplied by to correct for the wind velocity V^. I 2 ytho 

ionization rate between Xg and X^,is calculated by two formulas. For yields 

W * 100KT, I 2 3 - Kg >3 U) V M° Where fche values for K 2,3( l ) *^2,3 

and n were generated by fitting curves through the values of a Table in 

Chapter Three of Carl F. Miller's manuscript "Fallout and Radiological 

Countermeasures". For W < 100 a curve was fitted through the first three 

values Of I, , given in another table from Miller and adjusted for wind 
** > j v 

velocity by multiplying by The resulting equations appear on 

lines 07^1 through 08 C 2 in the computer program appearing in Figure'5. 


In all the equations that follow the scaling function formulas 
and constants are given in Table IV. For any scaling function dependent 




figure 5» Computer Program for Ionization Rate Contours 
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figure 5 (Continued). Computer Prog ran for Ioniration Rate Contour* 
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Figure 5 (Continued)., Computer Program for Zonitation Rate Contour. 
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Figure 5 (Continued). Coraputer Progma for Ionltatlon Rate Contour# 


lit) 

UH 

**i«a 

»n 

loti 

ntt 

ttt> 

it** 

iiai 

ini 

uu 

t»I 

1214 

• 

122 # 

1231 

1214 

124* 

12*1 

USt 

12*2 

12 ** 

111 * 

1)11 

1)42 

1)47 

1)47 

1)77 

1)77 

1)77 

141 * 




* V J *4 4 < tW4 •* W 4 4 4I4 »' 

V# • 7419 . lit I 

Il*t • I«S* • II) < 

*04 J * (1.1.♦» I 

UU) • toutiutl I 

4*1 TitI46CI.SCF1I I 

«W7*vt ceitro* j « ti.t.*7 > j . ro* j • (i,i,«» 4 , 

ro* 4 • 11.1«*» 4 JIJ1 . 4 * V4 • f4442## » 4 

FORMAT 4CFHlll.4ll).«J.,* •♦*41)*4#44,*| *,#41!•#»VA»X)».Of• *1*. 

IU.I<4 «7R».14l.»H **.41,4. 4 44 I 

COMMfMT COMruTtTION 0r VALUf* mmICm *tu #i :044141*1 U* *HH WDlXH 

VI* • FF«t#*«> _ .... . * 

T# • Ttt . lO#«l*l»/tOSU#2>7.. * 

111 • III) • «17» .... . * 

• • I0HTIH.1I • X47.X071 4 

117 ■ IOC I H7I / H*» I 4 

«lNl(X<l».1.4.Xt}t.12) 4 "17.*121 4 

i#Mlll(l>.l2).XI4t.ll*> | MS*.#(4t S 

SIMlIV|9I«I<$)*X|47*1141 • "J4.444I I 

*#M| t XI4), |(4I.*I7>»|(7I * #4l,*471 4 

UNtt#.4,121,7*.1.9 I #7*,i»7$» 4 

*fMtt«.I,1(71.0.1.4 « NT*.*7*1 * 

4*|T*mCCr»l * FORMAT ICF4I# |,.*I.,I1 - Xf X) • *«». 
#*.,X* » 1# • X7 • X# • «l*.»».»T.»tlH<>|1.»7.Cl0U0*il». 
*C#»##1*.«F0CI*.W<I 4 

COMMfUT UC NOV F1MQ TM| CNQFOlNT* OF IMF SEVi*AL f LI U*tl AND 

TM| FOCI OF fMf UR4I4D CL500 PLLIRIH I 

mao<okaici < IRRUT fAtcdoe.FQA j*ii.i.to:i i icon » 

FOR 4 • 11.1,IOC) 4 









Yi&'ir* 5 (Continued). Computer Program for Ionlta'lon Rate Contour* 
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figure 5 (Continued). Computer Program for Iofiliation Itete Contour* 
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Figure 5 (Continued). Computer Program for Ioniiation Rate Contour* 
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Figure 5 (Continued), Computer PrograHl tor, Xonlx&tlon Rate Contour# 
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Table III 




LIST CP 

SYMBOLS 


Natation of 

notation due 

dotation of 

Notation due 

computer 


primarily to 
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primarily to 

pro gran 
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prograa 

Killer 

A 


a 

ELLC 

Distance from center to 



0 . 0 
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contour In upwind cloud 
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(used to differentiate 
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studied 

LA3 

In & 

0 







Table III (Continued) 
LIST C. SYMBOLS 


Rotation of 
co^iter 
Tr-.o-raq 

LAO 

LI { ) 

“7 

IM 

M { J 


PHI 

<K 

RS 

SA 

8 VW 

V 

vw 

VUB 

WSI 

W 

X( ) 

XCl 

XC3 

xcu 


notation due 
primarily to 
_ Miller 

la a 

o 

In 1 ° 

1 

In V 

Array containing 
constants for 
scaling function* 

Uccd in equations 
involving X* 

V *6' f 7 

°^,3 or a i * 5 

R . 

SA(I) 

8 (^ v ) 

V 15 


Rotation of 
computer 
program 

XC 6 

XJPRC 

XPRS 

XX 

xxc 

xxrr 

X4U 

X7 

X8 

X07 

X9 

X 98 

yicc 


Rotation duo 
primarily to 
Miller 

7, 9 ~ 

X* 

X* 

X 

X 

lajfA 


In (V r 9) 

X98 2 . X 87 2 
X ? in feet 
Xq. 4 h feet 

*8 ' *7 
Xj in feet 

X 9 * X fl ln 

Maximum cloud width at 
chosen contour 


a 

YIC3 

Maximum 



chosen c 

T 

V 

YS 

Y 

V 


s 


YY 

Y 

X l’ X 2’ X 3» \> x 5 * 

YOS 

Y° 

V X 7' X 8’ X 9 

Y8 


l, 2 •* 



Y015 

V«) 

2. 

3, «* •* 

23 

5, 7 ** 


s 


20 

2 

0 
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Table III (continued) 


LIST OP SYMBOLS 
Continued 


Notation of 
computer 
program 


Notation due 
primarily to 
Miller 


Notation of 
computer 
program 


Z99WDlDVELOCrmiPH v y 299 YIEIDKX 


Notation due 
primarily- to 
Miller 


V 


* The slope and intercept, respectively, for a semi-logarithmic line 
through the two points whose coordinates are given. 

** The downwind distance (for the chosen intensity) to the fallout pattern 
between the points indicated. 
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table iv 

Data Table for Ionization Rate Contour Program 


for 

rt 

Hui 

m 

(a/b ) 2 

i 

0.486 

0.262 

In a/r 

a 

2 

1.070 

0.098 

In A 

3 

3.389 

0.431 

R 

s 

4 

2.319 

0.333 

H 

5 

3.820 

0.445 


6 

4.226 

0.164 

°^3 

7 

.O .509 

0.076 

In a 

0 

6 

2.880 

0.348 

X. 

9 

3.308 

0.496 

i 

10 

3-564 

0.319 

°5 

11 

-0.054 

0.095 

X, 

12 

3.850 

0.481 

6 

13 

4.255 

0.200 

X. 

14 

3.862 

0.586 

7 

15 

4.268 

C.305 

x„ 

16 

4.005 

0.596 

8 

17 

4.410 

0.315 

x„ 

18 

5.190 

0.319 

9 

19 

5.202 

0.3H 

< 

20 

3.223 

0.400 


Tor 

V 

A* 


Slim a 1 

a 

23 

1.371 

-0.124 

9 

24 

O. 96 O 

0.146 

*4 

25 

0.270 

0.089 


26 

-O .176 

0.022 

a 6 

27 

0.030 

0.036 

°7 

28 

0.043 

0.141 

“8 

29 

O .185 

0.151 

kJT 

30 

-3-286 

-0.298 

ya 

31 

-2.889 

-0.572 

vlX 

32 

-3.185 

-0.406 

5*a 

33 



K,X 

34 

-1.134 

-0.074 

6 a 

35 

-1.225 

-0.022 

K.JI 

36 

-O .989 

-0.037 

ra 

37 

-1.079 

-0.020 

Vc 

38 

- 2 c 166 

-0.552 

i/H 

39 

-0.431 

-0.014 

40 

-0.837 

0.267 


41 

1 

fO 

vn 

O 

-0.404 

*12 

42 

- 2,600 

-.337 


v 21 3-644 0.467 

5 22 4.049 0.186 
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variable, for example H (height of cloud), the program determines which 
formula should be used, The equation for H appears on line 06l3 of the 
contour program. 

The following chain is then started 

3 - r? (5 + F, V) line 0613 

FF(N, W) - EXP(lfT(X, V)) line 0454 


UFO?, W) - 2.3055851 L(h) * HI * M(tl) line 0 U 67 


The program then finds for W* $ 28 KT that F * 0 since 
V greater 28 is false and F Is a Boolean variable. Therefore for W • 5, 
one finds in the Table IV that L(5) • 3*820 nnd M( 5 ) « 0.445. The 
complete equation of H for W $ 28 KT la 

H • exp((2.302585l)(3*82O) ♦ log^l# • 0.445) , 

log^H - (2.3025851)(3*^0) ♦ log^w • O.W»5 , 

or 

log^H • 3*820 ♦ 0.445 logyjW. 


The variation of and X^ with wind speed depends on the wind modified 
parameter or (lines 0713» 0717)* The equations for X. end X, are then 


*2 • 


2*303 log o r 2 


* 1,2 


VWtf ’ 86 ; 36) 

(1 ♦ 9.273 X 10*\v y ) 


(59) 
(59a) 

(59b) 

(60) 


(61) 
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or 


a (a 2 - 1900) 
«• *+ o 

~~a * O'. 02 0‘ 


( 6 la) 


for 1= npH * lnd velocity. 

<lh« equation- for the distances X^ through X ? (lir.e 0995) are 


given by 


The 


For i - 6, 7, 8 and 9 

x t - 6 .C 0 x loV 3 ’** 1 * 5 a t , tf - 1 to 28 KT 

X : - 1.68 X 10 U y°* l6U a t , « - 28 to 10 5 KT 

atove equations also hold for ^ when 0 ^ le greater than a/h. When 


q 5 Is lc-ee thai or equhl to a/h, then X § lo determined by 


( 62 ) 

( 62 a) 


X 5 - 6»C-d X loV*^ -1.1*0 x loV '* 0 <n.OC#-** + *? » ^ - 1 to 28 KT ( 63 ) 


or 


lit/ , N loV-^oi , I.to x loV * 300 / 3 . 06 W 0 - 262 + ; V-28 to 10 5 KT ( 63 a) 

? 5 

v w 

jiL*t the Integrities 1 ^ I 4 , and are set equal to 33 mpS • 

The eqUA*. i0»-> for I;, «d 1 ? are on lines 0951 through 0961. Iq Is 
determined fn** a sesi-lornrl thole line between (X ? , I ? ) and (X ? , 1^ 
with the V* 4 '»e o: Xg as previously calculated. 

Al’U r x^ through X^ and their corresponding intensities have been 
,, cun then tc drawn of the Intensity profiles or elevations. 

x OUDG ^ il F • * 

, „ rh , 1(0 u^ir.-.t distance X for a 10 CT weapon is given in 
A plot or h-'t'io' V ~ 

Figure U. t'or any set o:' intensities the intersection (downwind distance X) 
of the concur cUlpse* v -‘ h V - 0 can be read off the graph. 
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In the computer program semtlog equations were found for the 
Intensity vs. Distance curve so that any pre-selected intensity contour 
intervals were Included la the computer printout. 

Next, to find the lateral locations or vidth of the contour 
patterns for the stem, the half width Yg for the stem fallout is needed. 

Yg is defined as the lateral distance from the center-lir.e of the sten 
pattern to the lr/hr at 1 hr (for 15 mph wind) contour (lines 1216-1220). 

A seallogarithmic line from (l 2 y Y ■ 0) to (1, Yg) is then calculated 
on line 1315- Using this equation the lateral displacement of any intensity 
contour of the stem can he found. 

The maximum pattern half width for the cloud, Yg, Is the maximum 
distance from Y - 0 to the lr/hr at 1 hr contour when the wind velocity 
is 15 mph. The calculation of Yg with its variation due to wind speed 
occurs on lines 105 *» to 1092 . 

To find the lateral displacements of the contour locations 
between ( X^, 0) and (Xg, Yg) a serailogarlthic line from (ij, Y • 0 ) to 
(1, D - 7 Yg 2 ♦ (Xg - X ? ) 2 ) is calculated. In lines 1236 through 13**7 
the X ♦ Y and foci are then found for the ellipses. Using the information 
lu the computer printout on ionization rate contour map can he drawn. 

(figure 6). 

B. Fractionation Humbert 
1. First period of condensation 

For a ground surface burst, the ratio of the total moles of 
liquid carrier to the molar volume n(x)/V is computed by Equations (32) and 
(33) with y - 10 MT or 10,000 XT, and T - l673°K. They yield 
t ■ 53 see 

(n(0/V]KT - 3-25 x 10 * 2 atm. 








\0 

















For & typical element Sr- 89 , the procedure of finding its fractionation 
number of flrat period of condensation r Q ( 89 ) is given below: 

Tecay chain of Sr- 89 : 

88.2 < Xr—Kb — w Sr 

Se-*-3r -' 



Here a neutron-emitter Is involved. I at the Police and Bnllou 
date have already included this effect, hence no alteration needs to be put 
on the number of atoms N(A, t) in their work. Moreover, the valuca calculated 
by them according to Clecdenin theory are generally considered to be more 
reliable it.d also used in the following calculations. 


Se 

Br 

Kr 

Kt 

Sr 

V’o k 1011,106 

1.04 


0.603 

4.6l x 10* 13 

k j ■ Uuij/vTKf 3,108 * i0 

32.00 


18.62 

1.4i8 x 10* 11 

N>, t) 

0.21 

309 

94 

3 

C<'A, t) 

0.21 

309.21 

403.21 

406.21 

MtJil 

1 + *J 

6.363 x 10* 3 

0 

4.791 

3 

r "LbJl 
u 1 * V 

1 K o 

6.363 x 10" 3 

(2) 6363 

4.797 

7-797 

r (39) 

0 

0.0)03 

0 

0.0H9 

0.01919 


Vito: data un.i<;tcd, usually r.'ellgible 


flic r \O'?) for ir, Kr, Pb un the by-product.,, 
o 

2, Stcoil period of coadvnrntlou 

The frvctior.stlon number of second period of condencatlon '.A, t) 
>p:nds very ra.tch or. the time and the e It rents involved in the decay chain. 
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When the time is the sublimation pressures of moot 

P.V 


elements are ctill quite significant, ar.d the values of are usually 
greater than those of 0-23 EW®. In the second period of condensation 

V 

4 portion of each such element remains In the vapor phase. Therefore, 
In computing, the ratio of these two values should be set equal to unity. 
As time increases, the sublimation pressures of all elenenta except the 
rare gases and As and Se, decrease sharply and eooo become negligible. 
Therefore, the second period fractionation number for a mass chain not 
Involving the rare gases and Aa and De will be greatly simplified} 


r; (A, t) 


1 - r Q (A» t). 


«*) 


Four typical examples are given belovi 

(a) 10 KT ground surface burst, 60 sec after burst, Kr (83)' 

Secay chainj 

.Se, 


Ce-*>As 




Br-*»Kr 


This Involves a split chain, but the Bolles and Ballou data also have 
taken this effect into account. 



Ce 

As 


*2 

Br 

Xr 

r 0 (83) 

— 

0 

0 

0 

0.0066 

0.0066 

l * r o 

— 

1 

1 

1 

0.993^ 

0.993k 

W(83-6o) 

— 

0,13 

21.80 

16.80 

12.30 

0 


— 

0.13 

21.93 

16.93 

51.03 

51.03 

K°.(l-r o )H 

— 

0.13 

21.80 

16.80 

12.20 

0 
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n° - 

0.13 

21.93 

16.93 

50.93 

50.93 

Pj (2)8216 

(9)3578 

34.53 

34.53 

0.2086 


* 10 7Al - 

7 AX 

7.41 

7.4i 

7.41 

7.41 

ZN" _ 

m 

1 

1 

1 

0.9980 

0.9980 

r;(83) - 

0 

0 

0 

0* 

0 

* Rotes r^(A, t) can not be a negative value. 

henoe the 

smallest Is 

zero. 


(b) 10 MT ground surface burst, 1.47 nln after burtt,Sr('i)(90). 

Decay chains 



er-—Kr- 

Br 

Kr 

—sr- 

Rb 

Sr 

V 

r o ( 9 °) 

— 

0 

0.036 

0.1378 

0.1378 

X * r o 

— 

1 

0.964 

0.8622 

0.8622 

H( 90 , 88 ) 

— 

67 

338 

106 

0 

a» 

— 

67 

405 

511 

511 

iT-d-r^l! 

— 

67 

325-8 

91-4 

0 

in" 

— 

67 

392-8 

392.8 

392.6 


( 4)3050 

00 

(6)4338 

0 

0 

y iA 0 

0.23BWKf X 10 

1.92 

1.92 

1.92 

1.92 

1.92 

ZN" 

231 

— 

1 

0.9698 

0.7686 

0.7686 

r;( 90 ) 

— 

0 

0 

O .0936 

0.0936 



(e.) 10 MT ground surface bu 


r o (T?) 

1 -r 

o 

:<( ?5,68) 
,N°-(l.r o )N 

* 108 

v * 

J 

Vp K° 

w 

r o<75) 


rst, 1.47 min after burst,Ge(75), 
t^cay chair,. 

C u ■" m in *• Ga- ~ ■■ Ce 


Cu 

Zn 

Ga 

Cc 

— 

l 

1 

0.6716 


0 

0 

0.3284 

— 

( 3)12 

0.0459 

0.0536 

— 

(3)12 

0.0460 

0.0996 

— 

0 

0 

0.0176 


0 

0 

(11)2208 

i.9e 

1-* 

1-92 

1.92 




(3)4263 

— 

0 

0 

0.0018 

— 

0 

0 

0.3266 


Koto: In computing Y y the individual mass number dependent conversion 
iuctor Q V A) for U 235 to u 23 ® should be multlpUed to each radionuclide. 

“J " O* a(A) can be calculated from Table I 


(d ‘) 10 W ground 6Jrfae * 6.77 »in after buret, S r (89).Decay chain 

Wt-s fhewn in Section (a.) 


Br 

Kr 

Kb 

Sr 

0.0303 

0 

0.0119 

0.0192 

O .9697 

1 

0.9881 

0.9808 








■‘.^4.^1 fx.'fi,' 1 .', 
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N(89,4o6) 

— 

=** 

84 

267 

36 

a 

— 

— 

84 

351 

ko 7 

U"«U-r 0 )X 

— 

- 

84 



O* 



84 

84 

6a 




os 

0 

0 

cr ci 



1 

0.2393 

O.'OO* 

r (89) 

— 

=. 

0 

0.71*86 

0.76-04 


3 . Cross Fractionation Number 

The gross fractional loti number, r^, or r Q , is defined by Ejuaticn C'?). 
The tvo Ionization rates, l f (0 ** nrt l fp ln tht3 formula are defined ly 
Equations ( 45 ) and ( 46 ) rocpectlVe^y* The computations are straightforward, 
though tedious, once the values of the fractionation numbers, disintegration 
multipliers and activities an? fclitwn. The variations in the ionization rates 
of normal fission products with tln e > or c Imply the decay of normal fls- 
products, for U-235, y.238 and Fu s 239 hive been calculated and are tab. _ ■ 

in Table V, From this tabic, tit*.* normal ionization rate of 1J-238 at H ♦ 1 
hour can be Interpolated and ol/lniftad as 

i f p{i) • 6.973 x 10' 11 (r/hr at 1 hr) / (fission/cq. ft.) (65) 

For a 10 KT surface burst, the condensed ionization rates for four different 
fireball rising times have been *‘t»»'p-*tcd and are listed with the values of 
gross fractionation numbers comp j* * d i rom them: 

t 1 *{l) r (t' 

’M* X' ' 

2,'|08 x 10‘ 13 0.388 


53 sec 
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60 sec 

3.255 x 10‘ 13 

0.467 

1.47 nla 

3.873 x 10‘ 13 

0.554 

2.15 ms 

5.0P9 x 10' 13 

0.721 

{since the fireball rlrirr 

time, t, 1 3 related to 

the particle 


meter. or, ty Equation *2x- r , r^(t) is also a function of 0 (, eooetlmea 
denoted as r„(t). 

•j'i Act 1 v ity •.Concent rvtf an tr. tensity Ratio 

Vlth Equalling. tiA (48), the activity concentration Intensity 
ju.tlo for any Isotope at v,y downwind location can be readily computed. The 
values oaf the ratio for ’Ax biologically Important isotopes in a 10 MT 
surface burst are c tle-JUt. &rA shown in Table VI, 
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TABLE 7 


Decay of floras! Fission Producte from U-235., U-233 and P*->»239 
( r/hr at 3 ft above an Infinite plnnc Tor 10* fissions r*~r sqT ft. ) 


Age _ _ U.235 U-2 38 _ Pu-239 


Year* 

Day* 

Hours 

Thermal 

Fission 

(8 Mev) 

Thermal ] 

Fission 



0.763 


( 8 ) 

(*) 
9977 v ' 

(8)9960 

( 8 ) 

9329 

( 8)8907 

(8)8750 



1.12 


( 8 ) 

6648 

( 8 ) 

6632 

(8 

6172 

(8)5966 

(8)5761 



1.64 


( 8 ) 

4149 

( 8 ) 

4153 

( 8 ) 

3827 

0)3592 

(8)3537 



2.40 


(8 

2453 

(8 

2484 

(8 

2256 

(3)2071 

( 8)2065 



3.52 


(8 

1410 

( 8 ) 

1450 

( 8 ) 

1303 

(3)1166 

( 8)1196 



5.16 


(9)8079 

(9) 

6418 

(9) 

7582 

(9)6642 

(9)7098 



7.56 


(9) 

4786 

(9)5014 

(9) 

4587 

(9)3986 

(9)4398 



11.1 


(9 

2964 

(9)3094 

(9 

2897 

(9)2555 

( 9)2821 



16.2 


(9) 

1804 

(9) 

11869 

(9 

1792 

(9)1626 

(9)1761 



23.8 

(10 

9716 

(9 

11094 

(9 

'1073 

(9)1010 

(S'1063 


1.45 

34.8 

(10 

6305 (10 

16428 (10 

16393 (10)6235 ( 10)6360 


2-L3 

51.1 

(10) 

'3730 (10)3786 (10) 

13817 ( 10 )j %9 (10)3811 


3.12 

74.9 

( 10 ) 

>2276 ( 

’10)2319 ( 

I 10 ! 

12365 (10)2470 ( 10)2362 


U.57 


(10 

11483 ( 

10)1524 

10 

>1556 ( 

10)161*5 (10)1546 


6.70 


(11 

19986 ( 

: 10)1031 ( 

10 

>1039 ( 

10)1099 (lo)l021 


9.82 


(11) 

16774 ( 

:il)6972 ( 

in 

16899 ( 

11)7244 (11)6655 


i4.4 


(11) 

14490 ( 

: 11)4587 < 

;n 

>4462 ( 

11)4650 ( 

11)4226 


21.1 


(11 

12910 ( 

In; 

>2940 ( 

[nj 

>2845 (11)2953 ( 

11)2660 


30.9 


( 

;n; 

I 1813 ( 11)1807 < 

:u: 

>1762 (11)1837 ( 

11)1645 


45-3 


1 

11 

(1061 (11)1039 ( 

I 11 

>1034 ( 11)1092 ( 

12)9777 


66.4 



12 

16055 (12)5807 ( 

12 

>5910 ( 12)6360 (12)5728 


97.3 



:i2 

»3676 (L2' 

)3497 ( 

[12 

>3559 (12)3896 ( 

12 ) 3543 


143 



;i2 

>2170 1 

12)2090 < 

12 

>2079 

12 2320 

12)2200 


208 



. 12 ! 

>1194 1 

[12, 

)ll64 | 

k! 

>1133 

[ 12)1287 ( 

12 )ll80 


301 



[ 13 ; 

>4874 l 

[ 13)4790 1 

[13[ 

>4733 

13)5707 ( 

13)5170 

1.2 

438 



[13, 

)1399 1 

13)1373 1 

13, 

1525 

13)2135 ( 

,13)1864 

• 78 

650 



14 

3884 l 

14)3758 1 

14 

5517 

14)9083 ( 

,14)7690 

2.60 




14 

2031 1 

14 

)1975 1 

14 

3160 

14)4964 f 

,14)4352 

3.00 




14 

1444 1 

(14)1432 1 

14 

2213 

14)2692 

,14)2594 

5-58 




(l4 

>1154 i 

(14 

)1158 1 

(14, 

)lC03 

[l4)l442 ( 

[l4)l6u 

8.18 




(14 

)1026 

(14 

)1021 1 

(14 

)1291 

(15)9971 (14)1225 

12.0 




15 

9432 

(15 

)9293 1 

(l4 

)l094 

15 3452 ( 

[14)1057 

17.6 




15 

)8310 

(15 

)8211 i 

(15 

)9l64 

(15)7377 ( 

,15)9160 

25.7 




(15 

)7l83 

(15 

)6987 

(15 

)7431 < 

(15)6219 1 

[15)7668 


Numbers In parentheses indicate the number of zcroo between t..e decimal point 
and the first significant figure. 
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vo 

3 

3 
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OJ 

UN 

UN 
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UN 

OJ 

A 

OJ 
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CM 

Ov 

^4 

OJ 
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UN 

CM 

r4 

-4 

OJ 

Si 

ro 

CM 

OV 

«n 
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D* Particle Site Distribution 

ttiere la almost no data available on particle size distribution, 

nor has any rigorous theoretical approach been made. However, using the scaling 

functions of the previously discussed simplified fallout model, estimates 

cm be made of the particle size groups that fall at any downwind location. 

Given the fall velocity and the initial altitude, the size of 

spherical particles falling through a standard atmosphere can be found 

in curves prepared by D. E. Clark.At any downwind distance X the 

first information needed is the maximum and mimimum a size for the particle 

groups arriving at that location. Because the formulae for a and a . 

max min 

depend on whether the a group descends from the stem or the cloud, it was 

decided that the X-coordinate intersection of the logarithmic line connecting 

(X^, Ij) to (X U . 1^) (point 3 to point U in Figure b) with the logarithmic 

line from (X^, 1^) to (Xg, Ig) - extended if necessary - would be the 

division mark. For any downwind distance lose than this intersection X 

coordinate, the source of the particle group ie assumed to be the stem. 

Then and for the stem may be solved by graphical methods or 

by the following computer approach. Essentially a Hevton-Raphsoa Iterative 

procedure is performed on the following function to find separately 

a and or . 
max mi a 


f(a) 


». (a • aft 

2.303(0 * ~) 
*Z 


a(z o a - V) 

log (X - -jy-r-1 * log a„ 


a ♦ 


^T 


( 66 ) 


where 


k & is the exponential constant for the fireball major axi3 exj/tnslon, 
Z Q is a yield dependent multiplier. 


t 
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kj Is ao empirical Inverse time constant, 0.011 sec*\ 

V y Is the viad velocity, usually assigned to be 10 mph, 

P Is the constant In the empirical, formula of Y^/V^, 0.95 whea 
Z • 5,C00 to 50,000 ft., d ■ 200 to 1,200 microns, 

q Is the constant In the empirical formula of V^/v^., 1-02 * 10*^ 
when 2 • 5,CC0 to 50,000 ft., d - 200 to 1,200 microns, and 
® 0 is the uajor cesl-oxis of the fireball at ground zero. 

The iterative procedure le 


0 L ■ Initial estimate 

O 

(67) 

“W * S ' 

( 68 ) 




If the sequence of a’a approach a limit, then this at is such that f(ot) i 0 } 

or a la approximately equal to o or a . depending on the sign of the 

iwu nu 

number enclosed in the absolute value symbol of Equation ( 56 ). 

When the downwind distance )» greater than or cquil to the X 
intersection coordinate, it is assumed that the fallout originates from 
the cloud. Por particle groups falling within the cloud area at any (X, Y) 
location, the maximum and minimum a values can be calculated byi 

a . —- ♦ <<? - v a ) h 2 . b 2 (i ■ v 3 /» 3 ) 

* h 2 - b 2 (1 . y 2 /. 2 ) 

where Y is translated to y by Equation (57). 

v v 

For any Ot value, the fail.rate is given by • JJ" • The height of fall 
for the stem fallout ia 

<y h • v> 


Z m 


(70) 





For the parr.lcle groups that originate from the cloud the maximum and 
minimus height of fall is given by 



Using this information calculated in the particle size parameter program. 
Figure 7, and a set of curves prepared by D. E. Clark, see Figure 8, the 
approximate particle size distribution at (X, jf) can be estimated. 

For a 10 MT detonation, and for various downwind distances the 

maximum and nini=r_n a; the falling velocity, V f : the initial altitude Z 

and particle diaojter u are given in lab Lb VII. 

(*»y 

Recent vork by Clirk and Cobbin' ' is closely related to the 
subject of this report. These authors used the same fallout model aa 
described herein for the calculation of particle size and radiation 
intensity at varying distances from ground zero. The results obtained 
are in good agreement with those herewith reported for weapons of similar 
yield. 



t 
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Figure ?• Computer Program fnr Parsunetera 
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Table VII 

particle Site Distribution at Various Downwind Distances for 5 WTT Veapon Yield 


_gn>_ u») 







- m - 

*T 

—nwj — 


f^»• d*•'1 

9 

— r~ 

- ran - 

•tlJI.I 


2^ r * _ e*J T 030 

iJ A 

_w 

i*ii 

7(0)Ufl 

10»’#*0«J 

JWMaOJtece 

. .1)5. 

10 

t»|H 

ItlMMT 


*m»a?2;co8 

m 

rf ■ 

-faiff- 

uur-'.i — 

—f;***?r*" 


. 

• rf ■ 

-1*1*’ * 

* fc—- 


”‘*»0»SJ|i?63 - 

— 



»«*C*.*T 

«•»*»«*! 

**4}>.?47000 


20 

*»i* 

■ i>»UJ 

Ud.fclftJJO 

4* >: ».«?'00f 



** 

-aril 

T* •••TJ. 

5TU7'T.^ 

T1**775?^ 

in 

71 

- HI* 



^'6T4*.9?30e “ 

—nr 

>0 

**ll 

nWtH 



__LE- 


*0 

*-h 


♦e** % rs;ioQ f±Q 


If 

. Amif ~ 

.»fi*irr. 

-ra~ 

« 


*71— 

ma 

- T»7ltni- 

—r.nim 

-rurrmiMo 

w 

*8 

un 

*•117*4 



- w 

♦ 0 

ii*l* 

W*7»Ut 


T, I**«4I'DO® 





*•*>)>** 

*»!)..»Mc0t 


4f 



1» 

41 

All# •• 

t#V41U 


m*vMH 

—jsr* 

II 

ArM 

#*f •*«*! 

*«4 4 .?>7 

*Wu’.r#CJ 


'1 

l«|l 




—«•_ 

» 


1«%* >!♦* 

"Jo'iCO’* 

»:♦**•* ip## 

to* 







)«• 

i*m 

U>2!' ,£( . 

|»«*Vi’__ 


_ja_ 

lie 

#■1* 

*«*i|4#l 

4IM)1« 

r.3c<iO;(*« 

-U1-. 

..ill.. 

1*0 

ion 

l*»|4 

.uufr.'A_ 

i#« 

**>•*•4*0008 

_12-. 

_Ji_ 





* * 

•lllr.HOCI — 

M 

/Cf ‘ 

—i«r* 


If# 

i*n 

ilt«l ill* 

«*•v*l. 

•oi •••***m 

.J2.. 

m 


(iirtitn 

l.;o» >*4 


«» 


• w~ 

-irn- 



• cr2T.vp.wf 


w 

— w( , -... 


..lino 

• nni.ijiKi’' 

r 


U 

5 


£ 


j 






























INITIAL ALTITUDE 


















■75' 


IV. IKPXED RADIOACTIVITY EC SOI IB 

Elements can be made radioactive by the action of thermal neutrons 

which excite the stable atoa making it unstable and radioactive. Ac thermal 

neutrons are liberated in a nuclear explosion, there will be some induced 

radioactivity In the soil around the detonation site. A small portion of 

this coil will be taken up into the fireball of the bomb (for a ground burst 

or near ground burst) and become part of the total fallout. 

A. V. KlsmonV^ in his report concerning the potential radionuclides 

produced In weapons detonations, lists the activities of the various nuclides 

26 

that are induced In the soil. A total of 10 neutrons rcr KT was assumed to 

27 

be liberated, which agrees within a factor of 15 with that of 1.5 x 10 1 

(l4) 

neutrons per KT calculated from Classtone . The composition of the coll 
V 0 . 1 t van used for those calculations Is given in Table VIII. 

Figure 9 gives the decay of the Induced radioactivity In coil. The 
data van obtained by taklog the contributions of the various nuclides given 
in Table IX and sunning then for the various times after detonation up to one 
year. 

Figure 10 , "Integrated Neutron Flux as a Function of the Slant 

Range In Air of 0.9 Sca.Devel Density for a 1-Klloton Explosion" chows tho 

neutron flux verouo distance for a 1 KT bomb. These curves, as given by 

Classtonewere Integrated over ar area vlth a 2,500 yard radius ond 

were used to correlate Klement'o values. 

Senftle and Charpioo^^ give a detailed discussion concerning 

Induced radioactivity. A method for calculating the Induced radioactivity 

Involving tables and formulas is alto given. This method was used to pre-paro 

' 13) 

the data In Table X In order to correlate the data given by Klcnent' 


and ehovn In Tabic XX. 











,er.*nt 

% by weight 

Blcmont 

£ by velgl 

2 

1.00 

SI 

28. kO 

C 

0.00 

P 

0.10 

0 

50.16 

K 

2.50 

Ha 

1.60 

C* 

5.70 

M« 

0.80 

Tl 

0.20 

A1 

6.80 

Mo 

0.04 



fo 

1.90 


D^ntlty: 1.18 ea/cu.en. 


Va^*V*1>£S££v 






Figure 9. Taui Activity of Soil 





TABLE IX 


Radioactivity In 3olla (Mc/MT) 


Huelid* 



0 

1 hr 

1 day 

1 month 

1 year 

IW* 

1620 

1430 

542 

a 

a 

m/ 7 

504 

7.05 

a 



Al 28 

83,670 

15.9 

a 



Si 31 

991 

773 

1.78 

a 


P 32 

1.21 

1.21 

1.15 

0.276 

a 

J4 O 






K 

a 






204.5 

194 

50.0 

a 


Ca** 

0.300 

0.300 

0.300 

O .263 

0.057 

c U9 

1100 

135 

a 



Ti 51 

156 

0.0015 

a 



Mn 56 

21 22 

1655 

3.62 

a 


P *55 

0.110 

0.110 

0.110 

0.104 

0.087 

Pe 59 

0.014 

0.014 

0.014 

0.0095 

0.00006 

a • less 

than one curie per megaton 





The formula for the radioactivity induced iu any element as given 
by Senftie and Champlon^"^ isi 

A. - (o flfh)(l - e' 0,693 t/T )(e -0.693 »/*) (72) 


where 


activity, In disintegrations per cccond, after the nuclide 
has been removed from the flux for a period 0 
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SLANT RANGE FROM EXPLOSION (YARDS) 

Figure 10. Integrated Neutron flux as a Function of the Slant Kango la Air 
of 0.9 Sea-Level Density for a 1-Kllcton Explosion 
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As a special case, -where t is toova to he less than 15S& of T, 
a further simplified formula was used. 

\ • a : ' 2 =| 2 ) {.- 0 - 693 0/,r ) ( 75 ) 

Values used for A r obtained by Senftle and Champion^^ are listed 

below; 


rtl. 

at 

Ha 

- 1-39 x 13 i0 

_ M» 

Ca 

- 1.73 x 10® 

Ml 27 

- 1.3 x lY 

ca 49 

. 2.54 x 10 7 

Al 2 ® 

- 4.79 x 'Jp 

ti 5L 

- 9-01 x 10 7 

SL' 1 * 

- 6 .74 x 10 7 


- I.39 x 10 U 

F 32 

- 3-69 x iP 

Fe 55 

- 4 .56 x 10® 

K* 

. 1.32 x Up 

Fe^^ 

- 2.26 x 10 7 


In calculating the data of Table X , a neutron flux of 0.91 x 10 1 ® 
r.uutrons/sq ee per megaton via assumed. The total weight of the coil involved 
was taken us 5.92 x 10^ grssu. IMs value was calculated using a 2500 yd 
rudlu3, a l*ft depth, an averige density of 1.18 gni/cu cm, and with a soil 
composition as that given !c Table VIH. The irradiation time (t) who taken 
03 or.e minute. 

The data presentee it Table X agree rcaeonubty well with the 
data presented by Klement la Table IX. This agreement is within an expected 
factor of 15 ss the neutron flux used vac a factor of 15 higher than that 
assumed by Klecvnt. 

Knnleviile^^ sta'.*s that the radioactivity of the fission products 
from u '-DcL'.f.r explosion overrun lows the InduceJ radioactivity by a fuctor of 
approximatvlo - 10®. Ibren tho^a this tends to make the contribution of the 
induced radioactivity almost r. j .-liglble there are conditions under which the 



induced radioactivity night play a detectable part. The induced nuclides 
that form soluble salts will be readily dissolved In both surface and 
ground water and therefore will add to the contamination of the water 

supply. 

It has been generally recognized that the activating neutrons 
will penetrate the roil before the blast wave arrives. Some of the 
activated coil will then be token up into the cloud due to the cratering 
effect and thereby becomes part of the total fallout. According to 
Dr. C. F. Miller ^ the induced activity in the coil will be about 
0.019 (l«9 pei cent) of the ionization rate due to normal fission (l fp (t)] 
1 ht3 fraction (fill be a small percentage of the total ionization rate. 
TMo estimate of the induced activity is applicable for n 10 Of, risslon 
yield, assuming 0.8 neutron capture per ficslon from U- 238 , rcgnrdleoa 
of weapon size. 

It has also become established that the induced activity will 
predominate for a total fusion bomb. For a JO£ fusion and jjOjt fission 
bonto, the induced activity will be relatively small, except for a short 
period Immediately after detonation. However, activated 6 odlum and 
mn^tcalua may contribute materially at a later time^). 

It la therefore concluded that the contribution of Induced 
radioactivity to the contamination of voter supplies will not be signi¬ 


ficant. 
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V. TRANSPORT OF FALLOUT PARTICLES FT SURFACE WATER 

Transport of the insoluble portion of fallout may be divided 
into four general phases: transport in air* transport in overland flow, 
transport in stream flow and transport In reservoirs. Each of the four 
phases may be treated as a separate topic. Of the latter three, only 
transport by stream flow nas been considered in any detail) overland 
transport and reservoir transport are mentioned only briefly. Transport 
in air is discussed in Section II of this report. 

There are three major limitations on any stream-flow transport 
analysis. These are: (l) the complex nature of sediment behavior in 
streams; (f?) the paucity of data concerning fallout, especially data 
concerning particle size distribution or the iclatlon of radioactivity 
to particle size; and (3) the obvious impossibility of obtaining utrcoa- 
flow data for every reach of stream. These limitations precluded an 
analysis of the quantity of sediment resulting from the fallout of a 
given burst, although work on the quantity of fallout la now being done. 
It could not be determined that a certain weight of sediment would fall 
on a water surface or what the particle size distribution would be. 

Computations arc available, however, from Glass tone ^ and 
the OCD fallout model which gives certain relationships between particle 
sizes, downwind ulctancc carried, and total activity. 

Glasstone chows that 99 per cent of the total activity of 
fallout is carried by particles of ^00 nieivno or less in diameter. 

This relation of particle size to total activity was derived from 
Figure 9-18? of CLaestone's boo)/^^ and is chcrvn as Figure 11 of this 
report. The method contained in this section is based on this relation¬ 
ship. 
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Particles of bC0 nieron diameter are on the order of magnitude 
of normal stream sediment* This means that the particles carrying 
essentially all the activity will remain in suspension, or at least 
will net settle out immediately. It follows that most of the fallout 
particles which land on the water surface or find their vay into the 
stream from overland transport will be carried by the stream and will 
present a hazard at 6ome distance downstream. 

A relationship of particle size versus downwind distance 
curried was computed In Section II of this report. (Table VIl). Data 
from that table has been plotted in this section as Figure 12. 

Hue following general assumptions were made: 

(1) Uptake by biological organisms la negligible. 

(2) The stream chemistry is not such as te dissolve ordinarily 
insoluble particles. 

(3) Effect of flocculation on eettllng rates is negligible. 

(1) Weapon size (5MT), wind speed (l5mph), arl fallout pattern 

are used to correspond with other cectlons of this report. 

(5)' Radioactivity or particles over 500 microns may be ignored. 

Further assumptions are made and discussed in the text. Some 
of these asc unptlons depend on the nature of n specific stream, some 
require further investigation, ond rome which ore considered negligible 
in an approximate analysis may assume a proportionally greater importance 
In a more detailed computation. 

An estimate must be made of stream width, volume, and surface 
area. Ibis cuiy be done by equating the cross-sectional area of the 
.itream (at & gaging station) to a rectangle of equivalent area whose length 


Cumulative pekentaCC W tOTM. *ctiviiv 


8 ? 



r**ticu w*mcii*. kichoh* 

Figure 11, Total Activity Percentage in Relation to Particle Diameter 
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figure 12. Particle Diameter versus Dovuvind DlGtaocc (MaxIbub and Minimum Curves) 
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is equal to the strean width. nils ia shown in Figure 13* length 

of the stream between any two sections ic considered to be a straight 
line equal to the distance along the 6 treaa channel between the two 
sections. Each section of the stream nay then be considered to be a 
truncated pyramid lying on its side with its bases the equivalent rec¬ 
tangles. (Or for the section containing the source, a complete pyramid, 

since the source is considered to te a point.) Concentrations nay be 

( 18 ) 

found using the volume of the pyramids. 

Intensity (roentgen/hour) contours for a given burst are super¬ 
imposed on the watershed. The Intensities may then be converted to activities 
(atorac/sq. ft.) for insoluble (and soluble) portions of the various lcotopes 
under consideration by use of Table VI. Details of this method are 
discussed on page 95 through page 98 of this report. 

If fallout landing on the water surface is acoumed to mix uniformly 
within the reach of channel enclosed by a given contour, an activity con¬ 
centration <aay be computed using the idealized channel shape. This concen¬ 
tration will move downstre nn as a unit, and the time of arrival and time of 
lapse of this concentration may be found at a point downstream by concideiirg 
the particles to move at mean stream velocity. 

By plotting the various concentrations from unlto under given 
intensity contours aa ordinate against time as abcioca, a cerlea of horizontal 
bars lo obtained. These may be averaged graphically to obtain a time- 
conccntrotion curve for the otrervn. !tote that concentration is exprcc.-.c-d 
as activity rather than particle concentration. 

By use of Figure 12 a maximum particle size at a given point 
downwind from the blast may be found. The maximum size particle landing 
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oa the stream reach under consideration is compared to the maximum particle 
size which the stream will carry In suspension. If the maximum size particle 
1» carried then it" la to he assumed that smaller particles will also he 
carried, and the method remains valid. 

If the stream will not carry the maximum size particle, then 
some percentage of the activity must be regarded as going into bed load, 
which travels at a much 6lcwer rate. The only way now available to find 
the proportion of activity which would go into bed load is by assuming 
that if particles of a given size are found to settle, that all particles 
of this size or larger will settle out, and all smaller particles remain 
la suspension. The relationship of particle size to activity given by 
Glasstone^ 1 *^ may then be used to find what percentage of the total activity 
would settle out with the larger particles. 

lfce ability of a stream to carry particles of a given size or 
smaller is best defined by what may be termed the effective carrying 
velocity of the stream, i.e. that velocity below which settling of particles 
of a given olze is most likely to occur. The effective carrying velocity 
is obviously not related to the idealized stream channel used to compute 
concentrations, but is rather a function of both channel characteristics 
(actual channel width and depth) and flow, as well as particle size and 
sediment distribution. Stream flow data, no matter how accurate or com¬ 
prehensive, cannot cover every mile of channel. Conditions which will 
produce the greatest possibility of settling must therefore be generaliza¬ 
tions of what is known of the nature of an individual stream. 

Particle Gize distribution in air is, unfortunately, not known, 
except for maximum and minimum distance which a particle of given size will 
travel. There is no way at present to compute the number of particles 





per unit area landing on the watershed. Sediment concentrations cannot, 
therefore, be computed. ThlB is the reason for computing concentration 
in terns' of activity. Additional fallout data, particularly on particle 
distribution, and additional study or. sedirr.cnt concentrations are needed 
to support the assumptions that the fallout particles mix .ir..u'or«ly, 
that particles will go out of curpension according to size, and that 
particleo irove downstream at mean stream velocity. 

In order to approximate actual stream conditions, close scrutiny 
of m'-.ps ar.d existing stream-flow records vill be required. Large-scale 
topographical mans should, yield much information on channel conditions 
and stro'-m-bed slopes. The U. S. ecological Surveys and Corp3 of Engineers 
publish or have on file not only gaging-etatlon records but much special 
measurement data. 

So analysis has been made for transport resulting from overland 
flow. Fallout landing on the land portion of the watershed is conoidcrad 
less of an Immediate problem than that lending on the water surface, 
and will probably present no Immediate problem unless It Is raining 
at time of fallout arrival, (in which case, incidentally, the fallout 
nvvy be intensified.) Subsequent rain may produce a second activity 
peak st a water intake by washing particles Into the stream. It Is of 
interest to note that at Covcetn Hydrologic Laboratory It war. found that 
much of the runoff which was formerly considered to flow overland actually 
gcc-u under the surface, or alternates between sub-surface and above-surface 
fj -• In work with radioactive tracers at Ooweeta, great difficulty was 
experienced lr. finding a tracer which would follow the water at the came 
ijved r.a the water. Most of the Isotopes used at Coveeta tended to otay 
near U.e point where they were applied. These studies were carried out 
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on steep slopes on which there was a higher percentage of surface runoff 
than on the same soil at a lesser slope; however, the Isotopes used were 
not ah sorbed within partlclea.^^^ 

Partielts which go into bed load nay present a long term hazard 
because of the slower and more Irregular movement of bed load. A bed 
load analysis will be required If It Is found that large particles, In 
spite of their possession of only a small percentage of the total activity, 
are found to actually carry a dangerous amount of activity. The smaller 
the particle size that Is found to settle, the more significant becomes 
the bed load analysis. Note that such a bed-load analysis complements 
assumption 5 on page 86. 

Transport In reservoirs Is more likely to produce settling 
because reservoirs are more likely to approach quiescent conditions. 
Calculations will probably be simplified because the conditions of reser¬ 
voir transport more closely approaches the condition of settling in a 
quiescent basin than do the other phases of the transport problem. 

This method was presented not only because It is the one most 
likely at present to produce usable results, but because It provides an 

outline of the problem of stream transport. Assumptions, questions as 
to their validity, and areas requiring further study have been pointed 
out. Computations have not been Included because it Is wished to ascer¬ 
tain the validity of some of the assumptions, and because sufficient 
study has not been made on specific streams. This treatment therefore 
constitutes only a good start; however, further study In the same manner 
should produce rcoults of known validity. 
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71. EVAIUATICS C? IffSTTCIPAI. WATER StTCTIX CCCuAMMAITIGB 
A. Introduction 

An approach to the problem of water contamination by radioactive 
fallout was taken which would give a no re reliable set of concentration 
values than those calculated previously. The values presented In Quarterly 
Technical Report Bo. - -were maximized, in that they represented a case in 
which the highest c once at rat ion. that could reasonably be expected at H + 1 
hour froa any kind of attack or envlronaent, except for possibly a very high 
megatonnage attack. PoLlovizg is presented a more detailed and realistic 
approach to the problen. The values herewith reported represent specific 
caces in which envirosneattl and other factors have been considered, and 
It is believed they are at least of the came order of magnitude as vould 
be expected following a nuclear attack. 

Only tha cities of Houston, Texas and Beu York were chosen for 
a full scale evaluation. Since the values calculated for these two water¬ 
sheds were in good correlation, it was felt that the results could probably 
be applied without serious error to the other cities previously considered. 
Excellent watershed data was available on both cities, thus facilitating 
the calculations cons lie rally. 

A physical Integration waa performed over the considered watersheds 
and reservoirs to give reaJLirtic specific Isotope concentrations at H + 1 
hour for variously directed lj aph winds. As before, tha case of fallout 
contamination from that falling in the reservoir alone was considered as 
well as the case of rmoff ctatanination free the entire watershed. 

B. Calculations 

Since large scale naps of both the Houston and New fork watersheds 
were available, transparent or*-lays were superimposed over the watershed 


See Appendix A 
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ureas tr.d appropriate grids scaled ir.. Thete grids consisted of five mile 
i-tj'-i ire a for water :hed areac, ar.d two xlle c qua raj for the reservoir areas 
when large es.o>igh to apply. Figure l.h bfcovc the grid on the watershed 
serving hcu-itoa- The gro.ind -^ro of a f- M? weapon detonation via located 
so as to give a fnllo'-t pa.hV.--n which fell on the watershed for the 15 mph 
wind .-colei. Vr-.riow- wind directions wore also chot.-n so as to give some 
bvis for comparison and to account for passible seasonal variation. It 
is obvious fron the geography of How York and Kou-uton that the wind 
directions frer the east and north would give rise to a negligible aiourt 
of fi.ilort over shell- respective watersheds, so that these wind directions 
were twitted. 

The attacks that were evaluated were U3 follows; 

Case la Houston - Crovr.d zero 13 downtown Houston ur.d the wind direction 
is fro* the south. 

Case Th Houston - Ground zero la downtown Houston tr.d the vlr.d direction 
Is from the vest. 

Case I la hew 'fork - Ground zero is Central Park and the wind direction 
is from the south. 

Case lib Hew York - Ground zero is Binghamton, Kev York, ur.d the wind 
direction is from the vest. 

tv pirate evaliiat Ions were made for ench of these cases with respect 
to cor. feci nan Ion from tne reservoir and that from runoff of the entire water¬ 
shed. T:i the case of Houston, the concentration of activity was calculated 
for one ol the f-.ed str:ons i'rsr. direct fallout co:,t.urination for comparison 
purposes. 

For the reservoir study, the ^re-a of the respective reservoirs 
within each two mile grid v»> found by piacireter. In the Houston watershed 









Figure 14. Method of Graphical Integration Used for the San Jacinto 

Watershed Serving Houston, Texan 










there Is only one reservoir and in the New York watershed all six rcoervoir* 
empty into one so that the problem is somewhat simplified if complete mixing 
is assumed. The centroid of each selected grid area wsuj estimated, and lta 
coordinates measured with respect to ground zero. An Intensity was then 
assigned to each set of coordinate^ and it was assumed to be constant over 
that specific area. The intensity values were chosen according to whether 
the specified area was either covered by the upwind or downwind oif'ud. 

Those values related to the upwind cloud were arrived at by interpolation 
between previously estimated intensity contours, and those of the downwind 
cloud were arrived at by interpolation of computer calculated contours. 

To convert from the intensity over each area to activity in 
atoos/sq ft in each area, the value of N*(A)/j^j from the fallout model, 
for each isotope considered, was taken as its value corresponding to the 
downwind distance X. By multiplying the appropriate n*(a)/j^j value by 
the intensity aver each square, the concentration in atomc/oq ft is found 
for each area. 3y then multiplying etch atom concentration by its cor¬ 
responding area, the total number of atoms in each square of the grid is 
found. The total number of atoms of each isotope at K ♦ 1 hour in the 
reservoir v&3 then found by adding up the contributions from each square 
over the reservoir. By assuming that complete mixing has occurred, the 
concentration of activity in the reservoir in atoros/litcr may be obtained 
by dividing by the total volume of water in the reservoir. The concen¬ 
tration of each isotope considered at H ♦ 1 hour, for the Houston and 
How York reservoirs is shown in Tables >1 and XT along with the equivalent 
activity in nc/al. 
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Table XI 

Activity Concentrations for Direct CootaMnotioa 
of the Houston Reservoir 



Couth Wind 

West Wind 

Isotope 

Atoits/iiter uc/ir.l 

Atoms/liter uc/ml 

Sr-89 

^•6 x 1C 5 1.9 x 10* 5 

b.6 x 10 9 1.9 x IO' 5 

Sr -90 

10.9 x IO 9 2.3 x 10' 7 

10.9 x 10 9 2.3 x IO' 7 

Ru-106 

5-1 x 1 C 9 3.0 x 10 * 6 

5.1 x 10 9 3.0 x IO' 6 

1-131 

9.8 x IO 9 2.7 x 10* U 

9-8 x 10 9 2.7 x 10'* 

Cs -137 

5-3 * 10 9 1 . 1 * x 10* 7 

5.3 x 10 9 1.1* x IQ' 7 

Ba-li *0 

13.7 x 10 9 2.3 x 10'“ 

13.7 x IO 9 2.3 X 10 '* 


Table XII 

Activity Concentrations for Direct Contamination 
of the New York City Reservoirs 


■I 

South Wind I 

West Wind 

Atoms/liter uc/ml 1 

Atoms/liter uc/al 

Sr-89 

2.9 x 10 10 1.2 X IQ - * 

3.9 X IO 10 1.6 x 10 "* 

Sr -90 

5.6 x 10 10 1.2 x 10’ 6 

6.9 x IO 10 1,5 x IO* 6 

Ru -106 

2.5 x 10 L0 1.5 x io* 5 

3.0 x IO 10 1.8 x IO' 5 

I-13I 

5.6 x 1 C 10 1.5 x 10 ' 3 

6.8 x IO 10 1-9 x IO’ 3 

Cr-137 

3.7 x 10 10 7.3 x 10’ 7 

5.1 x IO 10 1.0 x IO' 6 

Ba-ll *0 

7.0 x 1 C 10 1.2 x iO* 3 

8.1* x IO 10 l.U x IO' 3 
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Tke activity concentrations for the Houston reservoir (Table XT) are 
the earne for south wind and vert wind beenuoe the come intensity contour falls 
across the center of the reservoir in Doth caoes. When the wind is from the 
south, the right hand portion of the contour crosses the reservoir. With wind 
from the west, the left hand ptrtion of the contour is involved, 

A separate evaluation of stream eontamlnution vac mode using a 
randomly selected creek feeding the Houston reservoir for the purpose of com¬ 
parison only. This selected fred stream (Cyprce 3 Creek) ran perpendicular to 
the downwind fallout axis and emptied into the Houston reservoir. The evalua¬ 
tion was carried out in the rare manner ao that done on the reservoir. The 
stream va3 assumed to approximate on isocelco triangle with the base repre¬ 
senting the width of the streoa south. Two mile lengths were laid off and the 
trapezoidal areas cnlculuted. Intensity values and activity conversion values 
were chosen as before to obtaim the total activity falling directly on this 
strewn, To present an estimation of the concentration of activity In the 
stream, complete mixing was as rosed. However, it is evident that the activity 
will really move In a elug type of formation, completely mixing, if at all, 
only when the radioactivity retches the reecrvoir. Although some mixing during 
stream flow will occur due to turbulence and diffusion, the presented case 
will never be actually attained, since countercurrent mixing is unlikely. 

From the presented values in Title XIII it may be seen that some build-up of 
activity may occur in areas iirte -1 lately surrounding the stream fecd-ln points. 
However, this build-up will r.ot represent a curious problem in the studied 
cace because of its distance from the municipal intake, and the relative time 
elapsed from the time oC dctor.v. Ion to the intake time. 

The second case studied was that of the contamination to be expected 
assuming runoff from the entire watershed. Tills stuuy was made for both the 
New York and Houston watersheds, the method being similar to that of the first 
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Table XIH 

Activity Concent rations for Cypress Creek (South Wind) 


Isotope 

Atoms/liter 

jir/ rr). 

Sr-89 

2.2 

X 

10 12 * 

9.4 x 10* 3 

Sr-90 

5-4 

X 

10 12 

1.1 x 10 

Ru-106 

2.6 

X 

10* 

1.6 x 10“ 3 

1-131 

4.4 

X 

10* 

1.2 x 10* L 

Cs-137 

2.5 

X 

10 12 

4.9 x 10" 5 

Ba-luO 

6.8 

X 

10 12 

1.2 x 10" 1 


case. Five mile grids were used, however. Instead of two mile grids. Wind 
conditions were selected Identical to those of Case I. To calculate the actual 
amount or radioactivity (soluble fraction only) that reaches the reservoir, 
the runoff coefficients supplied by the respective municipal voter works were 
assumed valid In that they represent a maximum.value to be expected. The 
calculated activities are therefore considered to be maximum values. The 
actual radioactive runoff coefficient will actually be lesB than the aqueous 
runoff coefficient depending on (l) the Instantcous moisture content of the 
soil, ( 2 ) the duration of time froa detonation to rainfall, ( 3 ) lon-exchange 
and absorption in and on the coll, and ( 4 ) plant uptake, Thv calculated 
values of activity concentrations of the selected lsotopco at H ♦ 1 hour, 
assuming complete dilution by the composite likes and streams In the watershed 
for both Houston and Hew York Is presented In Tables XIV and XV. The factors 
Involved obviously vary widely from one environment to another, so that general 
assumptions made from ony particular reservoir should ho applied with caution 


to other situations 
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Table XIV 

Activity Conccr.t rat loss for Contamination from Runoff 
for the Houston Watershed 



South 

Wind 

Vest 

Wind 

Isotope 

Atom/liter 

nc/-1 

Atons/liter 

pc/ml 

Sr -89 

4.5 x 10 12 

1.9 x 10‘ 2 

3.6 x 10 10 

1.6 x 10* 4 

Sr -90 

8.0 x 10 12 

1.9 x 10* 4 

3.5 x 10 10 

1.8 x 10* 6 

Ru-106 

4.1 x 10 12 

2.5 x 10* 3 

4.0 x 10 10 

2.4 x 10* 5 

1-131 

9.2 x 10 12 

2.5 x 10* 1 

8.4 x 10 10 

2.3 x 10* 3 

Cs -137 

12 

f 5*3 x 10^ 

1.0 x 10‘ 4 

4.3 x 10 10 

8.5 x lu* 7 

Ba-llO 

12 

11.3 « 

1.9 x 10‘ 1 

. 11,0 x id 10 

1.9 x 10* 3 




Table XV 


Activity Concent rations for Contamination from Runoff 


for the Key York City Watershed 



South Wind 

Wcct Wind 

Ioot.ope 

Atoma/llter Mc/ril 

Atomo/llter uc/ral 

Sr -89 

0,36 x 10 12 1.5 x 10* 3 

1.6 x 10 12 6.8 x 10" 3 

Sr -90 

0.65 x 10 12 1.4 x 10 5 

2.8 x 10 12 5.9 x 10* 5 

Ru -106 

0.29 x 10 12 1.7 x 10’ 4 

12 -4 

1.2 x 10 7.1 x 10 

1-131 

O.65 x 10 12 1.8 x 10 * 2 

2.8 x 10 12 7.6 x 10* 2 

Co -137 

0.48 x 10 12 9-5 x 10* 6 

2.2 X 10 12 4.4 x 10' 5 
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C. Cone las Iona: 

It la not the purpose here to draw any concLualona In regard to 
the biological hazard resulting from the concentrations of activity presented 
In Tables XI - XV. This aspect is being considered elsewhere In this report 
and also by other investigators. 

It may be prudent, however, to consider the relative value of 
the previous calculations. It must he remembered that only a 5 KT weapon 
was considered and that all concentrations arc for H ♦ 1 hour. Thus 
for any other size weapon, appropriate revisions in the fallout model will 
...ive to he made before it may be applied. As for a multi-bomb attack, 
one cannot but accept some method based on additivity. Complications 
arrive here, aa far aa forming a model la concerned, becauao of ground 
zero locations. Since H ♦ 1 hour values of concentration nro given, one 
may apply suitable decay curves to adapt them to any specifio otroam flow, 
reservoir flow-throu£>, or intake time, etc. for a particular watershed 
aystem. 

It Is not to be Implied here that the Houston and New York values 
nre representative of all municipal reservoirs. However, they do present 
a reasonable value to be expected In a rcaliotlc case, and also demonstrate 
the ease with which an evaluation may be made. By making ouch an Integrated 
evaluation for any specific watershed and applying local conditions, the 
relative radiological hazard may be estimated. 

The Houston and flow York values will be seen to vary from the 
potential valueo If unusual conditions prevailed such as would occur 
during a dr/ spell. However, the effects vlll be sonewhat counterbalancing 
in that high reservoir concentration due to lav volume will bo somewhat 
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offset by the lack of runoff contamination. It la doubted that unnatural 
conditions, unless very severe, will change the concentration of activity 
by no re than an order of nagnltwie. 
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vn. ATIALVSIS OF RADIOACTIVITY lit WATER 

A. Radiochemical Methods 
1 . Introduction 

In most radiochemical analyses the steps are concentration, 
separation, purification, and counting. Evaporation and Ion exchange 
are the evo principal methods of concentration, but may not be necessary 
If the level of activity Is aufflclcntly high. The individual radio¬ 
nuclides sre usually separated by means of precipitation, solvent 
extraction, or Ion exchange elution. Purification of the Individual 
radionuclides usually follows stan'lard chemical procedures vhlch depend 
mainly on the solubility characteristics of compounds of the radionuclides* 
Coryell and »<garman have described (urlflcatlon of many of the more 
Important radionuclides. Radloicotopta are counted using appropriate 
instruments, depending on the specific Isotope and the desired accuracy. 

"Carriers" are of prime importance In vet radiochemical procedures. 
The mechanisms end usage of radioisotopic carriers have been veil defined 
by Overman and Clark A carrier, vhlch is normally the stable isotope 

of the clement being determined or a stable element vlth very slmllaur 
properties, lo added to lncrcace the total concentration of the element, 
because the radioactivity vlll normally be present In micro (juantltles only. 
Hence, if the carrier is not added, part of the radioactive material will 
usually oe lost, even though employing the moat scrupulous analytical 
techniques. By weighing the final solid to be counted, and knowing the 
original nnount of carrier added, the percentage chemical and physical loaa 
r.ay be calculated. By assuming complete interchange of carrier and radio* 
nuclide, the came percentage of each vlll be lost, and the original concen¬ 
tration of activity may readily be calculated. The carrier itself Is 



normally added as a solution of a soluble salt of carrier element 
prior to any chemical procedure. Carrier addition la used extensively 
throughout the described analytical procedures, although scoe carrier-free 
determinations are now cooing Into use. 

There are several compilations available vr.tch describe detailed 

(20 22 23 24 } 

procedures for the quantitative determination of many radionuclides v ' * ' 

tost of these methods have been proven and are In general use. 

2 . Procedures for the Determination of Specific Eallcnuclldea 
a. Radtostrenttua 

Strontium -90 has been generally recognized as being a very bio¬ 
logically hazardous radionuclide as It accumulates In bone. IXte to its 
relatively dangerous nature, vide lute rent has been aroused to find a simple 
and rapid method of analysis for otrontlua-90. Sb.'irlques that have been 
In use for several years are Still finding widest use because of their 
accuracy, the time consideration being less important la peacetime. 

Strontium- 69 > perhaps somewhat less of a hazard than strontlun-90, 
because of Its shorter half-life, will also be present In fallout contaminated 
water. Separate analyses of strontlum -89 and streetium -90 are thus desirable. 
Strontlum-89 activity Is normally found by subtraction of strontlum-90 from 
total radlostrontlum concentration. 

(25) 

( 1 ) Radlostrontlum and Radlobarlua by Mtrate Separation* 

Strontium and barium carriers are added v.d the Croup II cations 
are precipitated as the carbonates. Partial separation from calcium Is 
accomplished by nitrate precipitation In fuming nitric acid. The remaining 
calcium nitrate Is then extracted with acetone. Rare earths ar.d other 
trlvaler.t cations are removed by t vo hydroxide ocawr-glngs. Barium Is 


finally separated fron strontium by precipitation barium chromate. Stron¬ 
tium Is then collected as the oxalate and counted. The barium Is converted 
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f/'oin the chromate to the chloride for counting. The lowest level of detection 

*8 

f.,r radiostrontium by this method is U x 10 uc/ml and for radlobarlum, 
l0 ' 7 nc/ml. A precision cf about 1 C# cay be obtained. 

The calculation of the total amount of radioactivity due to either 
jailostrontiua or radlobariuz, using an internal proportional counter. Is 
t „fc.ie by use of the f clloving equation: 


uuc _ nut cpta f 

nr * 


Vijare 


A ■■ efficiency factor, 

B * chemical yield, 

C « aelf-acsorptton factor, and 

D * volume In liters. 


<iiila formula also holds true for determination of the other radionuclide* 

( 26 ) 

A variation of the preceding technique vae introduced by Kool v 1 
A/i.er the separation of strontium by carbonate and nitrate precipitations, 
t, 6 plum reaioval is effected by barium chloride precipitation in a hydrochloric 
iId-ether system. Any lantha.num -140 vhich may be present fj-om barlura-l^O 
(3r;i:ay io removed by a ferric hydroxide scavenge. Sensitivity in the 10 * 10 uc/ral 
,.»ngo is obtainable by this procedure. 

(f>) Strontlur.-r* by Solvent Extraction of Yttrium- 9 i 27) 

Strontium carrier is added and carbonate precipitation is performed. 
j-, r . carbonate precipitate Is dissolved in hydrochloric acid and the. solution 
(avenged with hydroxide, strontium io again precipitated as the carbonate. 
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vhich 1 8 then allowed to stand, thus permitting yttrium -90 ingrowth, n* 

yttrlua la extracted with tributyl phosphate equilibrated vlth lb» nitric 

«td. Storing extraction, the yttrlua Is washed from the organic phase 

vlth 0 .^ N dried la a plnnchet, and counted In an anti -coincidence 

beta counter. 

A deeontiu'lr.afcion factor of approxlmately ^ la obtalne<i> ^ 

chemical yield is greater than 70 f for strontium; the recovery of yttrlua 
by the solvent extraction is 85 ^, 

Strontium -90 determinations, cuch as that shown above, are normally 
nade by measuring the activity of its beta decay daughter, yttrium- 90 . After 
a suitable period ^ elapsed, the lnerowth Is separated from the 

bulk of the strontium by some suitable means. Ohe yttrluo -90 is then purified 

and counted. Zhe strontium may then be determined by calculation using the 
following fomula:^ 2 ®^ 

Strontlum- 90 , uuc/i • - net cpm_ , 

(aT(b)(c)(D)Tb>[F;(2.22) 

• efficiency factor, 

• per cent extraction factor, 

• per cent Ingrowth factor, 

• chemical yield, 

• sample volume (liters), and 

• decoy factor; calculated from 



where: 

A 

fl 

C 

0 

E 

F 


vhero: 


activity remaining after a tine interval t, 
activity of cample at core original time, 

_ CJ21 _-_- , and 

half-life of yttrlun-SW (c2.“* hrs) 


t • time (hours) fron separation to counting. 

3y knowing the total otrontium-90 activity ar.d the total activity dee 
to radiostrontium, the strontlum-89 concentration may be obtained from the 
difference. 

(OQt 

( 3 ) Strontium -90 by Direct Precipitation of Vttrlum- 90 v 

Strontium -90 activity may atso be easily measured by a variation of 
the preceding method. After a suitable time has elapsed Tor yttrlum -90 in¬ 
growth, yttrium carrier is added and subsequently precipitated as the hydroxide 
by the addition of codlura hydroxide. The radloyttrium la then purified by 
standard radiochemical techniques, which will depend on the concentration and 
species of the other rare earth radionuclides present. Resolution of the 
decay curve will be necessary if yttrium -91 is present. 

(U) Radlor.trontlum by Ion-Exchange Methods 

Kahn and Reynolds^ 0 ^ utilized ion-exchange resins for the concentra¬ 
tion of radiostrontlum. Strontium carrier is added to a 10 -liter water sample 
and the oolution is pasced upflov through a cation exchange resin in hydrogen 
form. The strontium is then eluted with l^K nitric acid. The radiostrontlum 
Is then separated from other alkaline earth3 and purified by previously 
mentioned methods,( 2 S)( 27 )( 29 ) Jepondlng on vhether total radiostrontium con- 

centratlon i 3 desired or strontium-90 alone. 

By means of concentration of activity with lor.-cxchonge resins, great 
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sensitivity can be obtained, possibly down to the 10 uc/ml range. 
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A similar iou-ex change method vas described by Kahn, Eastwood, and 
Using appropriate radiochemical techniques, other long-lived radio¬ 
nuclides could also be determined- In every case, the recovery of activity 
was greater than $9jf. 

Bryant, Sattiznhn rnd Warren^ ^ utilized an ion-exchtngc procedure 
to separate yttrium from strontium after a ten day ingrowth period. The 
yttri'# 1 was then selectively el^ced from a cation exchange resin, followed by 
a rca'isorptlon on another resin and finally counted. The radiochemical re¬ 
covery was greater than 97 $* Gravimetric measurement is not required when 
using this procedure. 

( 5 ) Radlostrontlum by Evaporation Methods 
( 13 ) 

According to Libby' strontium -90 may readily be determined in low 
concentrations by concentrating the activity by evaporation. 

The sample (l liter) is evaporated to dryness, and dry strontium 
nitrnl? carrier is added. Following addition of calcium chloride, to the dry 
residual the sample is dissolved in phosphoric acid. Yttrium ia then removed 
by the addition of lanthanum carrier, followed by hydroxide scavenging. The 
precipitate is redlseolved and the milking repeated twice more. The third 
precipitate is saved, Ignited and counted. The analysis itself takes two 
hours td perform. It is necessary to utilize decoy curves to determine the 
yttrlufii» 9 C activity, due to the presence of yttrium -91 and other rare earths. 

h. Radloceslua 

Padiocesium i3 considered to be one of the more important of the bio¬ 
logically hazardous radionuclides. Cesium is chemically similar to the other 
alkali petals, such as sodium and potassium, which are commonly found in the 
internal organs. Thus ionic interchange within the body is easily accomplished. 
Co -137 “!>d Cs-ljU with half-lives of 33 year3 and 2.2 year3 respectively 
constitute the greatest radiocesium danger. 
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Selective alkali precipitation and ion-exchange elution have proven 
to he the test net hod a of separation for the analysis of radiocesium. Deter¬ 
mination of the various isotopes appears to be most feasible by gamma spectro- 
ftcanr- 

(1) Radloceslura by Fhosphomolybdate Precipitation ^ ^ 

After the addition of cesium carrier to the sample, phosphoric acid 
and ammonium molybdate are added and the cesium is precipitated os cesium 
ammonium phosphonolybdate. After dissolving the precipitate in sodium 
hydroxide, the cesium is reprecipitated as cesium eobaltinitrite. This 
precipitate is dried at 100 °C on an aluminum planchet and counted .in an in¬ 
ternal proportional counter. A 93 $ removal of cesium activity was attained 
by the phosphomolytdate precipitation. Concentrations as high as 10 pc/ml 
and as low as 10 ” 1 ® uc/ml verc measured with good accuracy. The too principal 
cesium radioisotopes measured by this method vere cesium -137 and cesium- 13 b. 

( 2 ) Radiocesium by Cobaltlnitrlte Precipitation 

After the addition of cesium carrier, the cesium and other Group I 
cations vere precipitated as the eobaltinitrite after the addition df sodium 
nitrite and cobaltous chloride. After vashlng, the eobaltinitrite is dis¬ 
solved in hydrochloric acid. Cesium siiicotungstatc is then precipitated by 
the addition of silicotungstlc acid. The silicotungstate is then dissolved 
in dilute base and the solution scavenged vith ferric hydroxide. The excess 
tungsten is removed as the insoluble trioxide. Cesium and sodium perchlorates 
are then precipitated by the addition of perchloric acid and absolute alcohol 
to the solution. The sodium perchlorate is then removed by washing the pre¬ 
cipitate vith absolute alcohol. The final cesium precipitate is then vashed, 
dried and counted. This method has been proven to be quite accurate and is 
nov in vide use. 
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Osmond, et ad,™ co-precipltated cesium cobaltinltrtte on potassium 
cobaltlnitrite and followed vlth a similar purification procedure. Thle 
method Las found little popular -J«. 

(3) Padloceslua V.~ rillcotungstnlc Precipitation ^^ 

This procedure is booed re a method originally described by Yienagata 
and Yamigata^"^. The cesium Is first precipitnted as the sllicotur.gi tnte. 
Dissolution of the precipitate is then followed by n ferric hydroxide 
scavenging. The cesium is then ^precipitated as the dlplcrylunlnnte. This 
salt Is dissolved In 4-!Bethyl-2-ps=.tanone nnd the cesium Is extracted by moons 
of 2 M hydrochloric acid. The cerium 13 finally precipitated ns the perchlorate, 
In which form It is dried, velghei, and counted. The chemical yield lo 0 Of>, and 
eight analyses may be performed — eight hours. 

(23) 

( 4 ) padloceslum by s Co-crystnlllrntlon Procedure 

A method has been describe! by which cesium 13 separated from the bulk 
of the alkali elements and the miasd fission products by eo-erystallUntion 
with amrsonlum aluminum sulfate. The ammonium salts are decomposed by heating, 
end the cesium is precipitated for counting from a dilute hydrochloric ncld 
solution as the chloroplntlnate. This method has been in use as standard 
procedure In son* laboratories. 

(5) Rndloceslum by I.-c.-Hxchangc Methods 

* IT 1 

Kahn, Eastwood, and Lacy' have developed n separation scheme lor 
the analysis of the more hazardous radionuclides by ion exchange. The npeclflc 
mdloelencnts studied were cesium, cerium, cobalt and strontium. The lower 
limit of detection waa decrease! a hundredfold by concentrating the radio¬ 
activity of n large sample with s cmtlon exchange resin. 

Cesium was selectively el-ted by means of CH hydrochloric acid 1 rom 
the cation exchnnge resins stuuiei- Standard purification procedures were 










then U6ed to prepare the sample for counting. The cesium vas separated from 
the other alkali cation contaminants by precipitating It as the olllcotu-.g- 
state, dissolving the precipitate In sodium hydroxide, and then repreclpltattng 
It as the perchlorate for counting. The activity recovered for all radio¬ 
nuclides vas greater than 99*2^. 

( 37 ) 

Tsubota and Kitano found that an ammonium formate - formic acid 
buffer of pH 3-2 selectively eluted the alkali metals from a cation exchange 
resin. This buffer vas employed particularly for the determination of radio- 
cesium, being somevhat superior to hydrochloric acid or citrate buffer as 
an eluant. 


C. Radioiodine 

Radiolodlne may be either beta or gamma counted, although gamma 
scintillation Is preferred for determination of the Individual rndlolsotopea 
of Iodine. 

(l) Radlolodlr.c by Chemical Methods 

(llendenln and Metcalf determined radiolodlne activity by on 

extraction jXiriflcatlon procedure. Carrier sodium iodide Is added to the vater 

sample, and Interchange io accomplished by oxidation to the lodato vlth sodium 

hypochlorite In basic colutlon, folloved by reduction to the iodide by codlura 

bloulflto In acid solution. Sodium nitrite Is then added to oxidize the Iodide 

to elemental Iodine vhich is then extracted into carbon tetrachloride. The 

Iodine Is further purified and concentrated by back-extraction into sodium 

bisulfite solution vhich la finally gamma counted at the lodlnc-131 photopeak. 

•8 

The sensitivity of the method Is approximately 10 * c/ml for a liter sample. 

The problem of Incomplete carrier Interchange Is overcome by this 
procedure vhich employs oxidation-reduction. 
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The same authors also describe & procedure in which the iodate Is 
reduced directly to the iodide by hydroxylaalne hydrochloride. The iodide 
obtained after the sodium blsunflte bach-extraction Is precipitated vitta 
sliver nitrate as silver iodide vhlch is then dried, veiled, and beta 
counted. 

A similar but faster nethod for radiolodlne determination va3 de- 
veloped by Levis . He developed a continuous extractor employing tvo 
centrifugal pumps, to promote mixing ar.d extraction. The iodine was extracted 
into carbon tetrachloride and backcxtraeted into bisulfite solution In one 
operation. TRie gaaaa emission of a pipetted sample vas then measured. The 
time of analysis vae cut from 2 hours to 30 minutes. Ar. accuracy of 
93.0 - d.of va a achieved. 

d. Total Radio Rare Earths Determination 

(1) Total Rare Earth Activity by Fluoride Precipitation 

Home and Martens developed a rapid method for the determina¬ 

tion of rare earth activity by which the beta and gamma activities arc deter¬ 
mined separately. The method la quite rapid, the gamma determination taking 
only a half-hour. 

Tvo rare earth fluoride precipitations serve to remove zirconium and 
niobium, the principle gamma emitters of fission material. Barium and stron¬ 
tium, vhich arc heavily eopreclpitatcd, are removed by hydroxide precipitation 
after the addition of holdback carrier. Tne cample may then be gamma counted. 
To estimate tho beta activity, the rare earths must be precipitated as the 
oxalates before counting. 

(2) Total Hare Earth Activity by an Alternate Mjthod 

( 1 * 1 ) 

ocldrldge and Hume developed a method similar to that above 
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fcr the det* — ir.atlon of total rare earth activity. It Is somewhat lengthier, 

Vat zcre arr.rate, particularly far beta counting. 

-'able c-rium Is added as a carrier for the entire group of rare 

earths. !2at rare earths are then precipitated aa fluorides, rediasolved in 

boric and -_i*.ric acids, and repreelpltatcd as the hydroxides. Zirconium la 

then recover by precipitation aa the lodate after reduction of all the rare 

eartha to '-it trlvalent 6tate. Following hydroxide scavenging for alkaline 

earth remrrsl, the rare earths arc precipitated aa the oxalates and beta- 

counted It an ir.tem&l proportional counter. 

<U2) 

Alctad and Pappas* ' employed a similar technique, using lanthanum 
as a carrier for all the rare eartha. Following separat'on and purification 
of the gro-t ;f rare earths, individual separation vaa attained by elution 
from a catiot exchange column with omonlura lactate. The ion-exchange separa¬ 
tion step vti performed without additional carrier, and the recovery of 
activity ftc each radionuclide was assumed to be the same as that for the 
lanthanum- Thr: advantage of the non-carrier separation Is that a "weightless** 
sample is rd'air.ed, thus eliminating cclf-aboorptlcn and self-scattering 
factors. 

e. Kir—llancouq Radionuclides 
(ij Fad lobar jura 

> termination of -adlobariun r.ay easily be made during the radlo- 
8tro,ntlum ar-vlyrls. toth t;ron»lum and barium form Inooluble carbonates, and 
subsequent ::-*ration may bo made by precipitation of barium as the chromate. 

This cethoi. by :->hn and Straub, was doccribcd earlier under radio strontium. 

(1*3) 

Hunter and r-.-oina followed the tame procedure except that the precipitate 
Is counted *j tie carbcnatc Instead of the chloride. 
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Mlnkklnen^ has developed a procedure for the analysis of radio* 
barium In vhlch barium chloride carrier la specifically precipitated ns the 
nonohydrate from a concentrated hydrochloric acid-ethyl ether mixture. This 
precipitate Is dissolved In water to vhlch ferrto hydroxide io subsequently 
added at a scavenger. The barium 16 then Isolated ag the chromate. The 
chenlcal yield of barium chromate Is about 70 $, Tne barium chroouite precipi¬ 
tate la then set aside for 13 ** hours to nllov the barium -140 and Its daughter 
lanthar.um-lUo to come to equilibrium. 

(2) RqJiocerlua 

( 143 

Borgus and Engelkemelr ' have developed method of analysis for 
radlocoi 1 uni (cerlura- 144 ), employing a series of precipitations. After the 
addition of cerium carrier and alkaline earth holdback carriers, cerium was 
separated along vlth the rare earths by fluorid« p.ocipitatlon. Alternate 
hydroxide und fluoride precipitations are then employed to remove zirconium 
and alkaline earths completely. The cerium Is then oxidized to the tetrnvnlent 
state vlth perchloric acid and precipitated as the lodate. Thorium activity la 
then removed as the lodate following the reduction of cerium vlth sulfur 
dioxide. After scavenging vlth thorium end further rcpreclpltutlon, the cerium 
lo finally precipitated as cerouo oxalate, vhlch Is subsequently Ignited to 
the oxide. After velghlng, the cerium activity io courted on on Internal 
proportional counter. 

Tho percentage of chemical recovery Is only about 7 / 0 $, but the radio* 
active purity of the final proclpltnto is very high. 

It lo likely that this process could be shortened bccauce of the 

probable absence of thorium from fallout contaminated vatcr. 

(40 

Ames employed a similar procedure for the analysts of total 
radiocerium. The concentration of cerium -144 Is obtained by counting la- 
mediately, using a 217 mg/cm aluminum absorber. This shields all the 




_ </j WLl'i 


.V. V- 







cerlua betas, counting only those t rom prascodymium-lUt, the daughter product 
of cerlua-lbU. The difference between the total radloceriua concentration and 
cerium-lW* concentration la assum'd due to cerlum-ll*l, a neutron induction 
product. 

i 1 * x 

Barnes' 7 describes a mel'iod by which radiocerium Is icolated by an 
extraction method and counted for 0 *riua-l 4 W. Carrier interchange is accom¬ 
plished by means of an oxidation-reduction cycle. Two dibutyl phosphate ex¬ 
tractions are then carried out to remove the heavy elements. The cerlua la 
then precipitated as the fluoride^ oxidized to the tetravalent state, and ex¬ 
tracted into hexonc. Finally it id converted to the dioxide and counted for 
cerium-lMr by means of u beta counter using a 217 mg/c» aluminum absorber, 
ttie chemical yield varies from “pO* 

Kahn and Reynolds^ ^ have developed another method for the determina¬ 
tion of radiocerium activity employing the use of ion-exchange reslno. Con¬ 
centration of rudlocerlut and added carrier was attained by passing a ono-llter 
water sample through either Dowex^t-O on the hydrogen cycle or 1R-220 on the 
sodium cycle. Cerium was oclcctively eluted from the reoln with 25 ml of 
3M nitric acid. Purification of tit" rndloccrium is obtained by successive 
precipitations as the fluoride, th<> hydroxide and finally as the oxalate. The 

cerium oxalate is weighed to dete»tnlr>* carrier loss and then counted. 

-8 

Sensitivity of detection was incrcdotd to approximately 10 pc/ml. 

(3) Radlqruthcnlua 

Glendenli/ 1 *^ has developed f>n analysis scheme for the determination 
of radloruthenlum which hoo been utilised us a standard method of procedure. 
Separation of the radloruthenlum vl<b carrier ruthenium is accomplished by 
oxidation with perchloric acid to volatile tetrooxide which is subsequent¬ 
ly distilled off. Sodium bismuthnl 0 *6 added prior to distillation to prevent 





volatilization of the halides by oxidizing then to the oxyaclds. The ruthenium 
is absorbed into a sodiun hydroxide solution and precipitated in the forn of its 
lever oxides toy reduction with ethanol. The ruthenlun oxides are then dis¬ 
solved in hydrochloric add, and ratheniua la precipitated in the metallic 
ctate by reduction with magnesium metal. The ruthenium metal la weighed for 
carrier loss and counted. The chemical yield is about 65$. 

Melnlck^®^ , using an almost identical chemical procedure, describes 
a process for the determination of ruthenium-106. By using a series of heavy 
aluminum absorbers with a beta proportional counter, a correction factor for 
the presence of ruthenium-103 (half-life - U2d) may be obtained by extrapola¬ 
tion. Extrapolation for the correction la necessary, in that while counting 

a 

the betas from ruthenium-106 through a 210 mg/cmf aluminum absorber, the gansna 
rays from ruthenium-103 are also counted. If the activity due to ruthenium-103 
alone la desired, gcuma scintillation may be used with a 2,000 mg/cm aluminum 
absorber. 

Merritt^) describes a method or radioruthenium analysis in which the 
ruthenium Is determined In the presence of strontium, cesium and cerium by 
fuoion and extraction. Following the addition of carrier ruthenium to the 
water sample, the sample is fused with potassium hydroxide, sodium nitrite, 
and sodium carbonate at 550°C for two hours. The melt is leached twice 
with water to dissolve the ruthenato and cesium, leaving strontium and curium 
In the residue. Strontium, cerlim and cesium are then analyzed by previously 
mentioned procedures. The ruthenium is extracted from the leach, after the 
addition of periodate, with carbon tetrachloride. The ruthenium Is then 
stripped from tho carbon tetrachloride with 6t! hydrochloric odd. Magnesium 
Is then used to reduce the ruthenium to Its elemental metallic state in which 
form it Is subsequently counted. This method has proven satisfactory for 
soli leaching analysis, etc., but could veil be used on water samples. 
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(*») Rad lo z i rcon lua-Rad ton icb lea 

Steinberg^®) has developed & method for the determination of zirconium- 
niobium activity which has been tested ar.d proven. Following the addition of 
both zirconium and n lob lira carriers, the two elements are completed by the ad¬ 
dition of oxalic acid, thorium cxalate scavenging is carried out by the ad¬ 
dition of potassium chlorate to the acidified solution. The niobium oxide 
is redlssolved and reprecipitated as the cxlde hydrate. The water is then 
driven off by ignition and the nloblun oxide is weighed and mounted for count¬ 
ing. the zirconium Is precipitated from the original solution as the phen- 
ylarsonatc by the addition of phenylarsonle acid. Thla precipitate is dlosolved 
In oxalic acid, forming the soluble oxalate, the zirconium la reprecipitated 
as the phenylarsonate, which Is subsequently ignited to zirconium dioxide, 
tills precipitate Is weighed and counted. 3oth precipitates (niobium and 
zirconium) are beta counted on an Internal proportional counter. 

Brady and EngelKcmcir^^ have developed a phosphate method for the 
determination of activity due to zirconium-niobium which Is somewhat lengthier 
than the preceding method. Separation of zirconium 1 b based upon the precipi¬ 
tation of zirconium phosphate. The phosphate Is then redlssolved as the 
fluoride and precipitated as the phenylarsonate which Is subsequently Ignited 
for counting. Ilioblua is precipitated la acid solution, dissolved a6 the 
fluoride complex and finally weighed at.l counted as the oxide. About 10 houra 


lc required for the entire analysis. The chemical recovery of zirconium la 
about 75& and that of niobium io about S0$. 


A method Is alro described by Stanley'"' ' for determination of 


zlrconlum-95. Carrier Interchange io effected by the formation of the fluoro- 

.2 

zlrconate complex, ZrF^ . Bare earth activities are removed by lanthanum 
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fluoride scavenging. The zirconium lo then separated by means of three barium 
fluorozirconate precipitations. The zirconium is finally precipitated vtth 
mantle lie add from hydr ochlo ric acid solution and ignited to the dioxide. In 
which fort It Is weighed and counted. A chemical yield of about 755& io 
obtained. It is essential to not begin the analysis until zirconium-97 
{tjy 2 ■ t? hours) has had time to decay to a negligible amount. After the 
final precipitate has been obtained, counting should begin Immediately so 
that no appreciable nloblum-95 has time to grow in. A beta proportional 
counter should be employed for the counting. 

3. 3wiuj«i/ of Radionuclide Analyses 

For the convenience of the reader, additional references, as well 
as those previously cited, are summarized In Table XVI, 








1. Radiostrontium 


25, 26, 27, 29, 
30, 31, 32, 33 


1*9, *>3, 54, 55, 
56, 57, 6l> 62 


2. Radiocesium 


23, 31, 34, 35, 
36, 37 


29 , *i9, 54, 55, 

50, r4, 65 


3. Radioiodiae 


28, 38, 39 


55* W 


1*. Total Radio Rare 
Earth a 


4o, 4l, 42 


55* ^ 


5. Radlobarlum 


24, 25, 43 


6. Radioccrium 


30, 44, 45, 46 


29* 4» ^9, 63 


Radioruthenium 


47, 48, 49 


8. Radiozirconium 
Radioniobium 


50, 51, 52 
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3. Instrumentation 

Quantitative analysis of radioactive elements in trace concen* 

♦ ration* la done most conveniently by meens of r&diat ion-detecting instru¬ 
ments. . Due to the very small amounts °f radionuclides present, sueh 
Instruments are usually superior in npei and sensitivity to conventional 
vet chemical analyses. Hovcver, wheil complete chemical analyses are 
required, extensive use is made of vet chemical procedures for the separa¬ 
tion of various isotopes, as descrited °n pages 105 - 121, Such separations 
are not required vhen only gross fission product activity is to be mea¬ 
sured or when only one radioisotope it* present. 

If the activity of a sample 18 f-ficiently high, this activity 
ray be determined directly using a simple detector and ratemeter. At 
very low concentrations in aqueous efeHipl* 6 or for beta emitters, for vhich 
self-absorption in the cample is a tiihjor factor, preconcentration of the 
cample and specially designed instruments tuch as a flow counter may be 
required. 

A survey of commercially aV§lls!l* portable radiation Instrument* 
hue been carried out and the information obtained is summarized in 
Tables XI ^ and XjJC. These are all illPtruser.ts of the ratemeter type 
capable of recording count rates dowll t® a few counts per minute. To 
establish their usefulness for the determination of fallout product 
concentrations in water it is necessary to correlate the sensitivity of 
the equipment with different types of detectors and reasonable sample 
volumes with the level of contamination concentration expected in emergencies. 

For rav, untreated uqueous polutiona, Figure 15 shows the relation 


of nuclide concentration to measured d®&« rate for various isotopes of 
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Interest for a representative eonsrorciHl instrument. It ha? toen shown ir. 
Tables and X' that fallout product cencshfrations of the order of 

10 5 • |C tic/r.I r.i/ ce expected In the water supplies of typical watersheds 


durirg the early path-attack 
'£* y 

Xshr»’ rhovs *.i* relative 
survey meters vita a lo.v-st- 


period. Table M’.. token from l*»cy and 
detectability of selected radioisotopes vith 
scale sensitivity down to 0.01 mr/hr assuming 


ideal counter gi-tmetry. 

If it 'J merely desired to determine the presence or absence of 
these feign-level fallout coacent rations., It is evident that available 
survey inatroett.* sooul'd be capable of supplying this information, pro¬ 
vided that they are Iftted with & detector head suitable either for immersion 
into the liquid rumple or for flow-through measurement.,?. 


‘fr/vever, in sort cases the concentrations rwwared ore veil 
above the ICKF it-stlmx: permisaible levels of concentration in voter- ' 1 or 
close tc the limit; of detectability or the instruments listed in Table XV a 11. 
For this reason it seems necessary to develop end make available procedures 
end facilities for the sore accurate meas'.urement of lav-level fission product 
concentrations to ensure the safety of available sources of drinking water* 
Thic eevelopeeat car. follow twe paths: 

1. f-eeccccctration of the cample to increase the number of 

radioactive shea? to a icvel that is readily detected by conventional 

•• 

survey instruments; or 

2. Increasic; the eensitivity cf the detectors by an increase 
in r.ample volume, cackgrcmd reduction by edded shielding and coincidence 
circuit Arrangers* tie, or more intimate contsct between cample and detector 
volumes, t? in projor:lonal flow counters or .’n liquid scintillation 
counting. 


Best Available Copy 



Table XVII 

Approximate Minimum t-f-tectable Concent;~atlor.s In Water (uc/nl) 


Tor Survey Meters vlth Pull Scale Sensitivity 


idlolrotope 

0.01 mr/hr 

0.1 nr/hr 

1.0 rcr/hr 

10 mr/hr 

_ 106 

Ru 

5 x 10* 6 

5 X 10* 5 

.1 

5 x 10 

5 x 10*- 

Pa 1 ** 0 

1 x 10* 7 

1 x 10*" 

5 x 10 4 

5 x 10* J 

P 32 

5 x 10* 6 

1 X 10* U 

1 X 10' 3 

1 x 10*'' 

,20i» 

-i, 

-3 

, -2 

-2 

T1 

5 X 10 

1 x 10 J 

1 x 10 

5 x 10 

Zr 95 

1 X 10 

1 x 10* 3 

5 x 10* 3 

5 x JO* 2 

I 131 

5 x 10* 4 

1 x IQ* 3 

.2 

1 x 10 

5 x 10* 2 

Ka 24 

1 x 10* 5 

.U 

1 x 10 

1 x 10* 3 

1 x lo* 2 

K 1 * 2 

5 x 10* 6 

5 x 10” 5 

5 x lo* 4 

5 X 10* 3 

Sc** 6 

1 x 10* 3 

5 x 10* 3 

5 X 10* 2 

1 x 10* 1 

Cr' 1 

5 x 10* 4 

1 x 10* 3 

1 x 10* 2 

5 x 10* 2 

CO 60 

1 x 1O* 3 

5 x 10* 3 

5 X 10* 2 

1 x 10* 1 

Sr 90 

1 x 10’ 5 

_k 

l x 10 

5 X 10* u 

5 x 10*^ 

Au 193 

1 x 10* 5 

l x 10 

5 x 10' 4 

5 x 10* 3 

c* 137 

1 x 10* 3 

5 x 10’ 3 

5 x 10* 2 

1 x 10* 1 


Values extrapolated from Mcy, W. J. anl Kahn, P., "Survey Meters 
end Electroscopes for Monitoring Radioactivity In Wntcr," J.A.W-W.A. , 
*•6, 59 (195 1 *) 
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Preconcert ration Methods 

Many methoda iucooucrui.tai.joa have disadvantages which limit 

their use in on er-ergeocy. Concentration by evaporation, for example, 
requires a great deal of care and skill to o'oialu reasonably accurate 
results. The use of Ion exchange resins to rerove the activity from the 
sample can give very good rc-pt'oduceability and In some applications can 

//fQ \ 

be alnost completely automated. Basons and Lauderdale' ' have developed 
a continuous water monitor which employs a long C.M. tube at the center of 
a column of cation exchange resin. By employing two columns, one may be 
used while the other Is being regenerated with strong acid. With thlG 
system sensitivities of around 10* 7 uc/ml have been obtained. 

The mixed resin bed 13 potentially capable of sensitivities In 
the 10* 7 - 10*® uc/ml range when used with law background counters. 

Although the measurement of radioactivity on resin presents problems of 
self absorption in the cample, reasonable efficiencies and rcproduccabllitles 
can be obtained if sufficient care io taken In sample preparation. 

One method which may be employed to overcome self-absorption 
in the resin lo to detect the beta particle# with a scintillation resin. 
Little work hao been done on scintillation resins, but when technique# 
are developed this should become a very valuable tool in contamination 
evaluation work. 

A highly accurate way in which the radioisotope content of water 
can bo determined Is by radiochemical analysis. By chemically separating 
the various elements of Interest and then measuring the uctivity of each 
constituent, the activity concentration as well as the Identity of the 
radio contaminants may be obtained. This type of analysis requires a 
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skilled radiochealst and such an individual nay he diff'l calt to flnd ln 
emergencies. Solvent extraction methods have also beof» l,fled > but at 
present a method involving & mixed-bed ion-exchange re? in 6cee » to offer 
the simplest practical approach^^. 

Detector Development 

Improvement ln detector sensitivity can be aelti ,iVcd by a suitable 
choice of detector characteristics, detector geometry* ? l,d electronic 
circuit parameters. In order to measure the contamlnal* evel vater 
supplies the radiation emitted from the water oample muef- fce determined 
as efficiently and as accurately as possible. Because uf the statistical 
nature of the radiation emission, higher accuracies cat) b c obtained only 
by increasing the total number of counts regiotered by t-l* c detector, 
while at the same time keeping background counts to a F ’ or th i« 

reason, at low levels of activity individual counting dlrP uito or scalers 
aro superior to integrating counting circuits ouch as fal-croetera. 

If the radioactivity is composed largely of K*l.« articles, a 
mnjor problem In measuring beta activity in voter is the short distance which 
a beta particle can travel in this dense medium. It is; therefore* desirable 
to bring no large a volume of voter os possible into thb Immediate proximity 
of the detector. 

One of the most convenient monitors for this Jj'ITP° oe 1° ^he flow 
through type which allows the water to come ln contact the detector for 

a brief period and then be flushed out. The moot common d'-tector u=ed ln 
flow-through monitors is the G.M. iu'oe^^K Various coiH lg>ra„lor.3 cf the 
G.M. monitor type have been used to check cooling water *.lllusnt at r hny 
of the water-cooled reactors in the country. 



Lately, the geometrical flexibility of plastic selntllle.tor3 

ha3 increased the volume of vater which can be brought close to the active 

volume of the detector. Thor? are presently on the market, water monitors 

which use large, spherical scintillators capable of measuring activities in 

the range of 10*^ pc/nl^^. Another method which displays much of the 

geometrical advantage of liquid scintillation but has none of the chemical 

( 7c) 

problems is the use of scintillation fibers . Although this method 

lh 

has been used for the most part to measure C and tritium in water. It 

should yield good results with fission products and extend direct monitoring 

-8 

capabilities below the 10" yc/ml range. 

Liquid scintillation techniques have yielded very accurate results, 
high efficiencies and good sensitivity in the determination of radioactivity 
in water. This method may prove quite valuable in the hands of an experienced 
radiochemist. 

The limiting factor In the sensitivity of a detection device 
Is the background or "zero activity" response of the device. Both electrical 
and physical approaches have boon pursued in the reduction of this background. 
Physical methods include the shielding of detectors and the reduction of 
detector sensitivity to the gorana radiation which can penetrate the shielding. 
Since the cross cectlon of plastic scintillators is quite high for beta 
radiation while being fairly low for gamma radiation, very thin scintillation 
discs have been used to produce good beta detection effioienceo with a 
background of less than one half count per minute for a one Inch diameter 
sample 

The electrical method of reducing background Involves a shielding 
detector ar.d anli-colncider.ce circuits. The detection device is thus fade 
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insensitive to most of the background radiation which enters the detector. 
Most of the anti-coincidence detectors presently on the market reduce the 
apparent background to a range of 0.5 to 2 counts per minute. 

The only reliable method presently explicable to civil defense 

use for censuring the amount of specific isotopes in water is gamma-energy 

analysis. Most gemma analyzers are large ar.d quite complicated to operute, 

although son/C that are portable ar.d easy to use arc r.ow appearing cn the 

( 74 ) 

market. ' Vith a small amount of training civil defense personnel should 
be able to make good use of these instruments. 

Summary 

Preparation for the post attack monitoring of drinking water 
supplies for radioactive contaminants, may now be divided into four areas 
of required action. These include the installation of continuous monitoring 
equipment at treatment plants, dissemination of information on the use of 
available instruments for water activity, insuring the availability of 
battery-per.-ered equipment designed for water or general liquid radioactivity 
assay, and continuing the research on equipment capable of meeting emergency 
requirements. 

A large number of water treatment plants would probably still be 
operational after a nuclear attack. The danger of contamination from 
fallout would require that the water be continuously monitored for contamina¬ 
tion even at levels well below the maximum permissible levels. In this way 
alternative water sources could be located or decontamination procedures 
could be brought into operation as soon as the possibility arises that 
existing wa*er supplies become unfit for consumption. It would be desirable 
for this permanent equipment to be able to determine the relative concentra¬ 
tions of the various radlocontaalnants. 




wv.v.v.v.--.v.v.v.-v^ 
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Many commercial instruments which detect nuclear radiations can 
he set up to taeasure the amount of radioactivity In voter to a degree of 
accuracy suitable for emergency ur£« This measurement would require a 
given procedure and a cct of conversion tables for each typo of Instrument 
or class of Instruments. 

Of the widely available Instruments, the most common type, the 
C.M, survey meter, can bo adapted for emergency uce by either dipping the 
probe Into the water or holding t-ho probe close to the water surface, to 
measure radlo-contamlnatlon levels down to 10 ^ uc/ral of some of the 
Isotopes. Some scintillation metcl'P can reduce this lower limit below 
10* 6 yc/ml with the absolute limit depending on background radiation. 

Most lor. chamber 3 ur/ey meters are not as sensitive as the G.M. or sclntll* 
latlon types, although some of the more sensitive ones might be used when 
no other equipment io available* 

Companies marketing nuclear detection equipment were asked to 
supply detailed specifications on their survey-type instruments. A 
summary of Information received by January 30, 19^3 la given in Table XV. 'I. 
The Instruments are listed In order of decreasing sensitivity (m 4 nlmua 
full scale readings). 

Pertinent specifications In the tablo Include the maximum full 
scnle range, battery life, type of ludlatlon detected, type of detector 
used, weight, the manufacturer's model number, and unit cost of the 
instrument. The ability of on Instrument to discriminate between various 
types of radiation Is also ir.cluJ"!, ao veil as the type of radiation to 
which a given Instrument Is sens!'. Ive, is indicated by an "X" In the 
appropriate column. If an "0" appears under the Indicated type of radiation 
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it simply means that It can be shielded out by a movable absorber. In the 
event an instrument may be powered by standard 'D" cells, the number of 
batteries required is indicated. Instnasento which arc specified to be 
moisture proof or immersion proof are indicated to be "weather proof" by 
an "X" in the appropriate column. It is believed that the nbove cited 
characteristics are the tost Important to be considered in selecting an 
instrument for surveying and measuring water supplies and for other general 
uses. The rrianufacturer of each instrument may be determined by referring 
to the "Commercial Source" column In the table and then to the company 
number listed in Table XIX. 


Table XIX 

Equipment Manufacturers 


1. Atomic Accessories, Inc. 

611 V. Herrick Road 
Valley Stream, 'lew York 

2. Baird Atomic, Inc. 

33 university Rood 
Cambridge 39, Haas. 

3 . Eborline Instrument Corp. 

P. 0. 2ox 279 

Santa Fc, New Mexico 

U. Vlctorecn Instrument Co. 

5606 Heugh Avenue 
Cleveland, Ohio 

5 . Lionel/Anton Electronic Laboratories 
1226 Flushing Avenue 

Brooklyn 37, ilev York 

6. Nucleonic CoiTpomtion of Am.erica 
19c Dcgraw Street 

Brooklyn jl, New York 

7 . Radiation Co>intcr Laboratories, Tnc. 
5121 W. Grove Street 

Skokie, Illinois 


8 , Radiation Equipment and 
Accessories Corp, 

£65 Merrick Road 
Lynb rook, Hew York 

5 . Technical Associates 

l'*0 W. Trovldcncla Avenue 
Burbank, California 

10. Trace rlab 

lC'Ol Trapclo Rond 
Walthr.m 5 1 *, Massachusetts 

11. Iluclcar-Chlcugo Corporation 
333 East Howard Avenue 

Deo Plaines, Illinois 

12. Franklin Systems 
273 b Hillsboro PoaJ 
West Pain Death, Florida 

13. Nuclear Corporation of America 
T.'.strvr.’ U* and Control Division 

Rlchwood Place 
Dcnvtile. How Jersey 
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Using the Information in Figure If. to convert the instrument 

sensitivities In Table XVIII into equivalent concentration value*, Tatle XV'.I 

has been compiled to chow the concentrations of typical Isotopes which can 

• f C (\ 

Just be detected with lnatrJM.-r.tc of the sensitivities stated' . It mur-t 
be bomt in mind that the dote calibration in mr/hr for mod Instruments 
is done with a cobalt or radl-ja ground source, ar.d that the dorc extrapoia* 
tlon to other Isotopes and tr.crc' rangoc will not be the same for different 
detectors. 

It is evident from the foregoing diccurslor. that theix; appear to 
be adequate commercial ln. tru.ner.*.a or. the- msrkc-t which can be adapted to 
give a reliable Indication of high-level contamination of water supplies. 

The at thy of icotopic contamination ar.J the accurate detcralnatlon of cafe 
concentrations of radloicotope-s in uater by field instruments under emergency 
conditions, however, doeu require .-.orec further research and development 
work. Such equipment vil). have to be made available in sufficient quantity 
for extent-lee water unaly.-es at. soon after on attack as possible-. 
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VIII. DECONTAMINATION OF HATER SUPPLIES 

A* dlceuosed la Quarterly Technical Report #3 (see Appendix A) 
rodloIogTcnUy contaminated waters have been subjected to decontamination 
procedures of both conventional and non-conventional nature. Of the six 
radioisotopes of biological importance, all have been found to be amenable 
to removal to some degree. The degree of treatment required has been 
predicted in previous reports, predicated by the Maximum Permissible 
Concentration for peacetime consumption of water or exposure over a con¬ 
siderable period of time, These requirements are very unrealistic when 
considered In the light of the occurence of a nuclear Incident. The 
removal efficiencies for the most part have been derived from laboratory 
data on small volumes of water treated under rigidly controlled conditions. 
Of practical Importance, only those data obtained by actual treatment plant 
operation are to be considered reliable on a major scale. Even in this 
respect, it must be pointed out that these data were derived from analysis 
of water supplies which had been contaminated by long-range weapon debrie. 

It has been discussed elsewhere in this report that the long-range fallout 
la more soluble than short-range or early fallout. It is for this reason 
that we may expect full-rcaie operations to be more efficient in radio¬ 
isotope removal of early fallout, than the presently reported efficiencies. 

It appears that water decontamination in the early phases of a 
nuclear incident aftermath are not beyond .the realm of practical application 
of the basic sanitary engineering principles of water treatment. The major 
problem involved in voter decontamination will be the long-term exposure of 
the surviving population to concentrations of the long-halflived isotopes 
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of biological Importance. The reduction of these contaminants to (l>-f'treablc 
levels vlll constitute the primary long-term problem to the water tiVntia'.nt 
field, but Is by no means an insurmountable task* 
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IX. BCC? 3UKDHI FROM C0KTAMINA7ED WATER 

In the analysis of tody burden from continued ingestion of water 
contaminated with radioactivity resulting from nuclear attach. It has been 
found that the developr.er.t of a rational formula is desirable. A formula 
has been developed in this report which gives values of body burden com¬ 
paring closely with those obtained by empirical means. 

Many of the empirical values used in this analysis were derived 
from animal studies, and cay or may not be applicable to man per se; 
however, the values obtained by formulas developed in this report are in 
good general agreement with data collected from experiments with man. 

Data from man has been used wherever possible; unfortunately experiments 
with man are not sufficiently comprehensive to servo as the entire basis 
of on analysis. 

The relation of body burden to many of the radioisotopes found 
in fallout is not easily understood for short periods of time, that is, 
when the study period is a very snail fraction of the physical half-life 
of the radioisotope. In this report, in addition to study of the immediate 
effects, several cases are analyzed for periods of time In excess of fifty 
years. (Figures 18 and 20) 

Hunwrous studies liave been made on the affinity of certain 
types of tissue for particular radioisotopes. These selectivity studies 
ore of special Interest to anyone who wishes to set up theoretical values 
which may not bo exceeded without permanent damage, that is, such values 
as "l-taxlnuB Permissible Concentration," "emergency limits," etc. Although 
it is not within the scope of this report to determine such values, it is 
of Interest to note that allowable body burden ia the total concentration 
of a 6lveu isotope allowed in the most critical organ. The moot critical 
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organ is that organ which receives its damaging dose first during the chain 
of biological processes. 

The phenomenon involved in the rate of turnover of various isotopes 

by a particular organ Is of Interest since this rate nay be equal to, 

(7 

greater than, or less that the rate for the entire body. 1 '' a sample plot 
of the critical organ versus physical half-life is shown in Figuro 16. 

Both biological and physical decay rates are considered in this 
report. In order to demonstrate the difference in body burden levels ob¬ 
tained by (l) considering only physical decay and (2) considering both 
physical and biological decay, l.e. effective decay, values for both 
cases are plotted in Figures lo through 21 as curves B and A respectively. 

values plotted ao curve A are parameters of the rate at which 
the radioactivity level In the body la reduced by a combination of bodily 
processes and physical decay. These values are therefore parameters of 
the effective decay cates of the various radioisotopes. The amount of any 
radioisotope which passes through the body must further be related to the 
amount retained by the critical organ. This relation is represented by an 
f value which gives the amount of a radioisotope taken up by the organ 
critical for that Isotope. 

In the lntcreot of generality ve have developed a method of compu* 
tntion (using A R values) rather than a Hot or specific figures for a limited 
number of cases. 

ilo attempt has been made to evaluate cell damage resulting from 
various levels of body burdens. It is obvious that the greater the body 
burden the more extensive is the radiation damage. Examples are cited in 
this report of sooe of the biological blocking processes which may bo 
applied to effectively reduce the body burden. These A n vbIucs will give 
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the highest burden revolting aiilJ fro ® drinking contaminated vater while 
eating a normal uncontcminated d* et and breathing uncontaninated air. 

With an. easily derived scaling factor, however, these same graphs may be 
applied to ingested solids. A /( values nay be obtained either from formulas 
to be derived later In this see'. Ion or from the graphs in Figures 16 
through 21. The graphs are easier to use and are included as an aid to 
computation. 

No limitations of physical death or early genetic death of the 
individual under study have be?fi placed on the data shown here, since 
both these Limitations require Sft analysis of damage. We cannot overlook 
the possibility that At eertaih levels of contamination, or consumption 
physical or genetic death will l* mit the continuity of our data, as will 
any departure from the size and nK a of a standard man. 

Establishment of a Criteria 

(a) The data to be presented should be in the most useful and 
general form applicable to the g'lbjcct matter. 

(b) Analysis of damage that rcti j r occur at the various levels of 

body burden is not to be attempted; however, comment regarding damage in 
general will be incorporated to dc P tb to the conclusions. 

(c) The development Q r a rat ional method for evaluating the body 

burden is to be outlined and e/pl alr * cd ' correlation of the values 

Obtained by the use of this im-tifOd and those obtained by empirical means 
will be discussed. 

(d) Assumptions oho/M fce QS f<fv end as reasonable as possible 
in order to maintain the lntegfMV of the data to be presented. 

(e) All values used !' or the various calculations shall be 
those derived from the most o-j*. norltatlve information available. 








Assumptions 

(a) The biological decay rate used Is that of a standard man. 

31r.ce the age arid overall metabolism of the Individual will affect the 
turnover and replacement rate of radioactive materials in the body, and 
therefore the biological decay rate, it is necessary to standardize these 
parameters for the present study. 

(b) The contamination of an individual’s water supply takes 
place all at one time, and no further contamination occurs. The analysis 
of a fluctuating contamination level would require that we assume ratea 
of fallout, wind direction and numerous other meteorological conditions. 

In order to retain the generality of our Investigation we have taken 

a case in which weather and attack data were not opplicatle. cur/es 
may be added to each other if new contamination occurs (possibly from 
fallout from e. different burst or fallout return from wind reversal), 

(c) The radioactive contamination of the water supply is reduced 
only by physical decay. This assumption is conservative in nature, since 
it disregards entirely the amount of activity which might settle out or 

be consumed by life forms found in the water supply, e.g. plankton, algae, 
fish. 

(d) Effective decay rate6, even if extrapolated from animal 
data, are applicable to man. In the absence of experimental data on man, 
the animal data 'used was the best available. It nwy be of in f erest to note 
that many contemporary investigators are not satisfied with many of the 
extrapolations of animal body burden data that have heretofore been accepted 
as correct, or at least useable, for the case of humans. 

(e) Tne daily intake of water is constant and takes place 
instantaneously. Virtually all of an Individual's water intake occurs 




in sixteen hours, and over half occurs in eight houro; however, the error* 
associated with this amplification tend to he compensating if the individual* 
drinking hah its are consistent* For reasons Of generality vc axe forced to 
consider that the man under study is a creatul 0 of habit. By drawing a 
smooth curve through the dolly points we have essentially integrated 
over the time periods between the daily values and have described a contln- 
uous function as if the individual had spread hi® intake over the eutire 
day. Tho combination of the two effects result 0 i n a close approximation 
of the actual build-up of body burden. 

Development of a Rational Method of Analysis i 1 *' body Burden for Continuous 
Consumption or Water Contaminated with RadioRi'Uve Element* 

Based on the Gelger-Nutall Law of radioactive tlecay we may say that the 
amount of activity present at any time 1* equal to 

At 

A • A e v 
o 

when only physical decay is considered. iS physical decoy constant 
for a particular radioactive element, and A q la amount of activity 
initially present. 

Furthermore, if only biological decay 1* considered we may say 

that 



where is the biological rate of elimination constant. 

If ve define "Effective Half Life” p« th<? time required for a 
radioactive element fixed in the tissues of nil animal body to be diminished 
by fifty per cent as a result of the combined action of radioactive decay 
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and biological elimination, then we nay say that the effective half-life 

»< 7 « 


B iological half-life x Radioactive half-life 
Biological haU-life + Radioactive half-life * eff 

An effective decay constant (k eff ) may be defined-/^ 



Let us assume initially that for any activity taken into the 
body the entire amount will be fixed in tissue. It is then apparent that 
immediately before initiation of consumption of contaminated water, th-' body 
burden iu t-quai to whatever radioactivity is present a3 a result of pre- . 
vlous activity. For the development of our argument we ahall take thi 3 
Initial body burden to be zero and correct it in later developments for 
the caoe of Initial body burden not equal zero. 

The following reasoning may then by applied: 

(1) Just before initiation of consumption activity in the body 
equals zero. 

(2) Immediately following the 3tart of consumption the activity 
in the body equals A^, that is, the activity per unit volume of the water 
(nc/ml, ntorou/gal, etc.) times the volume of water consumed. Since we 
have assumed that the amount of liquid is constant we moy drop the volume 
term and ^incorporate it later in the development. Furthermore, by 
considering A q equal to unity we will develop an equation which when 
multiplied by the actual Initial activity will give activity present at 
any given day n. Wn wish to develop a plot which is independent of both 
amount consumed and initial activity, and is dependent only on the radio- 
lrotopc ill question. 
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(3) Immediately prior to consur.ptloa on the ee» ,c,n< l 3a y the body 
burden would be 

A i “ V 


and Just after consumption the body burden is 

A eff (l) -*Jl> 

AJ . A* eff ♦ Ae p 

L O Q 


(U) Just before consumption or, the third day \H a body burden 


would be 


- AJ. ' k ' frl - (A o , A «"‘ « 


and immediately following ingestion, the expression for* t' ,v V burden may 
be written 

•1 1 -A l A l A 1 A 

Ag » (A q 6 P )e P + (A Q e eff ♦ A o e P )« 


( 5 ) Continuing this reasoning, a recursion fO)W>l a 1° developed 
of the fora 

A . A I ."*'" 1 ♦ « * V ‘ 

n o 


♦ e 


-u eff + (n*l5\.lj -n\_l 


•nil 

P " ♦ e ] * V P 


However, rather than deal with the (n + l) terms nccebsHf/ 1® evaluate 
this equation at n days ue simplify the e-^.tion os folld '"' 0 


h n ~ a 

A . A b (H- ' . r .) 
n o b - p ' 


where h • A t 
0 




and 


p - A o e 


*y i > 
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(6) This general expression for A q is the value plotted as 
curve A in Figures 16 through 21. 

(?) We may now incorporate into our formula the portions of total 
tody burden which have thus far been neglected. 

(a) Since we have considered A q to be unity, we may find 
the total body burden by multiplying A q by the true original activity, 

Aq tnte » the volume consumed daily, D. 

Therefore Body Burden, (ggS) - A q ^ (ggl) x D(M^) x 


A Of consumption * 

n body ' 


(b) The assumption that all the activity taken into the 
body is immediately incorporated into the intestine is highly improbable. 
We may however, consider the amount Incorporated into tissue ns n fraction 
of the ingested activity reaching an organ (f y ) and multiply the body 
burden of (a) above by f^ to obtain a more reasonable value. 

Thus Body Burden.. * f x A xDxA 

*1 w o true n 

Some values of t for Isotopes under consideration are given in Table 

(c) While reducing the bedy burden by the factor f we have 
neglected the activity, A^,, which is present in the intestinal tract 
on day n. Thio activity is an addition to the body burden, therefore: 

Body Burden, XI - x * 0 true x D x A r ♦ A Q , 


where by the reasoning and definitions previously expressed 




o true 


A (n-0.25) -1.25V 

P )(e ^ 






Since the tine required for ingested material to pass through the body 
is 31 hours or about 1.25 days, the A^j tem is strictly transitory. 

Tt nay or may not be negligible depending on the period in which it takes 
place. On the first day, for example, it would be the largest portion of 
the total, but by the time the fractional half life (n/T^) is equal to 
orie, the term may in some cases be considered negligible, (n equals 
time in days since beginning of consumption and equals radioactive 
physical half-life in days.) 

(d) The only tern still missing from the mathematical expression 
Is that for the activity present in the body at the lime drinking of 
contaminated water started (Aj). That which was already present in the 
tissue at time of initial ingestion will, after n days of effective decay, 
equal 




V 


A 

e 


rf° 


That present in the gaatro-inteetinal tract at the time of initial ingestion 
will give an additional activity equal to 

A J( • f y x Age offU 

where A^ is the pre-ingestion activity in the intestinal tract. 

(8) After combining the various components, the expression 

for body burden from any particular radioisotope at day n may be written: 

-\ (n- 1 . 25 ) -1.25V 

Body Burden ■ (f^)(A o tntrue® ^ 

(*! • Vo>* • * A ot * \ * \ 
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The full explanation of this development has been presented in 
order to allow the user of the quatlons to decide which, if any, terms 
may be neglected for any specific analysis. Should the individual In 
question be an average IT. S. citizen who has sought shelter during an attach, 
the term 

*i * f A ' " *T * *x 


will for practical purposes equal zero. As previously noted, the value 
of A QI will generally be negligible after a relatively short period of 
consumption. If the (A^ ♦ A^) and terms are dropped, the practical 
expression for body burden is (f y H* 0 true^ D ^ A n ^* 117,6 A CI tcrn » h< * reve *“> 
may be of major ir.jvrtance when the physical half-life io much greater 
than the effective half-life. 

It muot nlco be streoced that to evaluate the effect of Aqj, 

A t and on any individual that individual’s hlBtory of exposure to 
radiation must bo known. It is only by thorough Knowledge of the develop¬ 
ment of this rational method that one may make a reasonable calculation of 
body burden. 

Reduction of Body harden by Blocking 

It has Iron found that in certain coses uptake of radlolcotopos 
may be reduced by biological blocking. For example, by administering a 
soluble, non-toxic, inert compound of iodine one can prevent uptake of 
radioactive iodine Vy the tlyroll. Ten to one hundred milligrams of 
potassium iodide token at one time will prevent the uptako of a ay more 
iodine, radioactive or otherwise, by the thyroid for about one week. To 
continue this effect, it is believed that an additional '.00 milligrams evory 


two days would be required 


( 77 ) 


Blocking lovers the voluo of f in body 
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burden calculations. 

Id the case of Strontium 89 and Strontium 3C, which behave like 
calcium, it 1* believed that blocking nay be completely or partially 
accomplished by consuming a calcium compound far in excess of normal 
requirements. Ibis in, however, a conclusion made on the baaia of extrapo¬ 
lated data; the only case in which it has definitely ceen proved that blocking 
Is effective is that of iodine. 

Should the body be Initially deficient in a given stable isotope 
(possibly from deficiency of a common mineral in the pre-ingestion diet) 
the value of f for the radioactive isotope of that mineral or a chemically 
.timllar mineral may be lncreaoed because of the increased affinity of the 
body for that mineral. 

Blocking should not be considered a panacea to the problem of 
body burden, since it may result in the cure being worse than the disease. 

In the iodine example, for instance, the prescribed mounts of potassium 
iodide constitute a daily Input of iodine into the blood one thousand 
times normal. Although it is believed that for the majority of the U. 8. 
population this treatment would not be harmful for to to three weeks, 
the secretion rate of the thyroid is decreased. Extended treatment will 
produce other still more harmful side effects vhich bare not yet been 
fully evaluated. Blocking, however, may play an important role in com¬ 
puting body burden valuco and should not be overlooked when selecting 
f values. 

Discussion 

(l) In figures 16 through 21, curve A gives the value resulting 
fron the effective decay rate. Curve B, included for purposes of comparison, 
represents the rate at which physical decay alone take* place. It accns to 





point out that because of the body's tendency to eliminate certain radio¬ 
isotopes, or not to incorporate then into tissue at all, body burdens are 
significantly lowered. 

(2) T^ble XXTY gives a fev selected values of f for six 
radioisotopes and various organs. These are values which result from a 
normal uncontaminated diet. Values of f nay vary from this table according 
to the excess or deficiency of certain minerals in the body. Note that the 
concentration of specific Isotopes must be known before body burdens can 
be computed. 


Table XX 

Selected Values of t for Various Organs 


Isotope 

Organ 

f v 

Sr 09 

Bone 

2.5 X 10* : 

Sr 90 

Bone 

2.5 x 10* : 

Ru 106 

Kidney 

2.0 x 10* ! 

x 13l 

Thyroid 

2.0 x 10* 


Whole Body 

1.0 x 10° 

C.137 

Liver 

U.8 x 10* 

B. 1U0 

Bone 

7-0 x 10* : 


Motes: 

(1) The f values shown are for critical organs. The critical 

organ will determine the A since It Is the one used in computing X ... 

n cx 1 

(2) The f value for the entire tody is the ruunnation of f 

v w 

values for all the Individual organs concerned. 
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Figure l6. vereuo Tlisc for Iodine 1J1 
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FRACTIONAL HALF-LIFE 


Figure 18. versus Tine for Strontium 90 















Figure 19. A q versus Time for Ruthenlua 106 


IV kvr »> VMUJir,* 



































L53 












-154- 



« i 

FRACTIONAL HALF-LIFE 


Figure 21. A n versus Time for Barium l40 
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XI. GLOSSARY 


A 

A 

\j 

H 

A o 

a ; 

\ 

\ 


a 


enpirical constant for sublimation or vaporization reaction 


radioactivity after come period, of decay 

radioactivity already present in the Intestinal tract at time 
of ingestion of contaminated water 

radioactivity present In the intestinal tract at some time after 
Ingestion of contaminated water 

radioactivity already incorporated Into body tissue at time of 
Ingestion of contaminated water 


a parameter of the effective decay rate in the body of a given 
isotope 


initial radioactivity per unit volume of contaminated water 
taken Into the body 


(o f Nk) 

t-%1 A 

MT* 


■ 


radioactivity present In body tissue resulting from decay of 
that activity initially present in tissue (A^; after some time 
of effective decay 


activity at time t 

activity, In disintegrations per second, of the radionuclide 
In a target, after the nuclide has been removed from the flux 
for a period 9j 

radioactivity incorporated into tissue from that activity originally 
present In the Intestinal tract <V 

cloud horizontal semi-axis 


a fireball horizontal eemt-axis at ground zero 

o 

*2 fireball horizontal semi-axis at 2 

B empirical constant for sublimation or vaporization reaction 

B ratio of fission to total yield 

Body Burden radioactivity incorporated into body tissues as a result of 
Ingestion of radioactivity 


b 


cloud vertical semi-axis 






c 

D 

D 

D X 

d(t, 1) 
FD r (t) 
FP r (t) 
f 

f w 

e 

b 

yt) 

ljft) 


k 


k 

k 

a 

\ 


G10GSAXC 

(Coatlaitd) 

fireball vertical send-*/!* at ground zero 
fireball vertical semi~tala at 2 

empirical constant for zublination or vaporization reaction 
instrument response faetor 
daily water intake, volume 
exposure dose 

decay correct ion factor for H ♦ 1 hour 
fraction of device contour ratio 
fia6ion-produet contour ratio 

thermal flux in neutronr per square centimeter per second; 

fUgacity of the element in the liquid phase 

a coefficient relating invested radioactivity to Body Burden 

acceleration due to grertty 

cloud center height 

fallout intensity at time t 

air ionization rate per fission at 3 ft above an. infinite, 
ideal plane for a uniform distribution of the normal fission 
product mixture 

air ionization rate per fission at 3 ft above an infinite, 
ideal plane for a uniform distribution of the condensed 
fission produdt mixture 

air ionization r*ta per fission at 3 ft above an infinite, 
ideal plane for a uniform dlstrioution of the neutron Induced 
activities 

relative abundance of tie isotope from which the radionuclide 
is formed 

arbitrary constant 

fireball horizontal set!-axis expansion constant 
fireball vertical semi-axis expansion constant 
Henry's Law constant 




maaatam 
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GLOSS ART 
(Continued) 


M r (t) 

a 

m 

N 

V 


»J(A) 


n 

n(i) 

n. 


n 


i 


n ,5A (t) 




n 


HT 

inverse time constant 
mass contour ratio 

mass of fallout per unit area at any downwind distance, X 
mass 

disintegration multiplier 

the total number of atoms of the element In the target; 
number of atoms, or moles, of flesion products per unit area 
mole fraction, of element J la the liquid phase 
mole fraction of element J in the vapor phase 

number of a tome (of the end member of cass chain and condenccd era 
the outside Of tho particle) which land per square foot of ground 

number of moleo of vapor 

total moles of liquid carrier 

abbreviation of nj^(t) 

amount of element J condensed on the surface of the solid 
particles 

amount of element J In the vapor phase 
abbreviation of "Ja«> 

number of moleo of element J with mass number A dissolved 
In the n(i) moleo of liquid carrier, which is the particle 
In the liquid phase, prior to solidification 

number of moleo of element J with mass number A, which Is 
mixed with nolcn from other mass chains to form n moles of 
vapor and not condensed in the liquid carrier 

empirical constant In particle falling velocities relation 

sublimation or vaporization pressure 

partial pressure of the liquid 


<1 

<*X 

B 


r o (A, t) 
r;(A, t) 
r x (0 
T 
T 




t 

t 


V 

V. 


w 

x 


sublimation pressure 

<V9»; 


empirical constant in particle falling velocities relatioa 

terrain shielding factor 

molar gas constant 

initial fireball spherical radius 

fractionation number of the first period of condensation 

fractionation number of the second period of condensation 

gross fission product fractionation number 

absolute temperature 

half-life of the radionuclide formed 

time required for a radioactive element fixed in the ti66ue 
of an animal body to be diminished 50 per cent as a result 
of combined action of physical decay and biological elimination 

radioactive (physical) half-life 

time of irradiation 

time 


time of arrival for particle falling out of the cloud 
time of cessation for particle falling out of the cloud 
particle falling time in the fireball 
particle rising time in the fireball 
molar volume 

average particle falling velocity 
vind velocity 

instantaneous particle falling velocity 
weapon yield in Xlloton TNT units 
longitude axis in ground coordinates system 
longitude axis in cloud coordinates system 


x 





GLOSSARY 

(Continued) 




Jr M (t) 


latitude axla In ground coordinate system 

total amount of radionuclides for the entire chain yield of 
=«33 meter A 

total tis,y .int of radionuclides of element J present at time t 

total number of atoms of element J of all mass numbers 

maximum half-width of the crossvind distance on a 1 r/hr 
contour 

latitude axis in cloud coordinate system 

amount of radionuclides of element J and mass number A 
present at time t after fission 

number of atoms of element J of mass number A not condensed 
in the liquid carrier 

altitude axis of ground coordinate system 

yield dependent multiplier in fireball altitude equation 

altitude axis in cloud coordinate systea 

scaliest Intercept of particle falling slope with cloud 

largest Intercept of particle falling slope with cloud 

particle size parameter 

maximum or minimus 

X/h 

a constant expressing the biological rate of elimination of 
radioactivity 

a constant expressing the rate of attenuation of radioactivity 
in the tody as a recult of both physical decoy and Moiogical 
elimination 

disintegration constant} the fraction of the number of atoms 
of a radioactive nuclide which decay in unit time 

activation cross section in square centimeters for 00 m/sec 
neutrons 


0 


tine of decay 
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APPan>.TX a 

Ecrr-pts f’-o-a Qu-.rt r iy Te.obn-.caI P rp crt tl 

F• rroVi.nltJ~.ry Evaluation of Wotcr Supply Ccyitamln at lrr. for Selected Clttea 
1. Ir.t redder lor. 

Tr. order to have some taels fc-r evaluating tho relative hazard of 
nuclear fallout to the water supplies of different cities, attacks have been 
assumed on c-ach city thv-t would result In the ssrr.e intensity over each 
watershed invol/ed. 

From this, the centred ra tion of 'he voter due to the suluble portion 
oi the fallout has beer, calculated for K 1 hour. Simplifying assumptions 
were made in regard to an average intensity over the specific watersheds and 
c.j to the time factors involved. There ssctjnptior.; were made to provide a 
hasie fc-r the calculations of the activity concentrations shown in Table IV. 
There values are shown principally to d-. mono t rote tV.c approach being taken 
5nd to give s.-a iiua of the maximum contamination that could be expected 
immediately following a nuclear emergency. 

The approach to be taken in future study will be more specific. 
Factors that will be taken into account to give comeuhat more reasonable 
values, will be . .'l) a realistic sup.-rimposlt ion of intensity contours upon 

tho specific votersheJt with regard to eu.sor.ai wind directions and target 
sites, (i) consideration of decay factors to gi.c c.vprcted activity at times 
hater than ’1 + l hour, ; 3) consider* f lor. of vat-, r-hed characteristics such 
as runoff co^t l lclent, feed-stream •ciocltiee, ivser/oir draft, transport 
pV :.c-.w-n», nnd reservoir mixing propvities, 'L) integration of intensity 
contours to give h no-~e precise value for up. clfic isotopic activities or. 
the watermen e.r d (5) a conversion of ictetialty to activity for sites off the 
dc.ro 1 r.t ar,lu er.ptoying a known rcia*.;or.uhlp to tho predetermined values 
directly do-noir.d. 
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T&fcle IV 

Activity C.ar.ccnt rations lr. V»'jt<*r Selected Isotope? 
Folio'-.;i»t 3 MI s.r attacks 

Hailonnclide lorctr.tratIon in Vinter (uc/ml) 


**- i v.’ 

Sr -09 I :-90 R-i- 15 c I-l 

31 

0.-13- *a-l : -0 

T ir.vcr 

i 

2.7 x 10 *“' 3 .- X 10 *-’ 3.0 x 10 *- 3.0 X 

10* 1 

1.9 >- n** 3.2 x 10* 1 

1’ ->lo. 

ii 

f.l x 10* 1 t .5 >• 10* 5 5.6 v 10*' 5.7 


3.5 X 11*- 6.1 

"* * 1iih 

•» 

7.3 x 10 ' h 7.3 x 10 * 6 e .2 x 10 * 7 8.0 x 

X0* J 

4.9 >. 1C*- 5.7 X 10* 3 

Mi t.n. 


lot Applicable 



to:: 

- 

6.5 x 10* J 6.1 x 10* 5 5-0 x 10’* 5.0 x 

io* 2 

3.1 x 10*" *.6 x 10* ? 

le/s.- 

ii 

6.5 x 10 1 1.1 x 10 ■ 7.0 x 10* ' 3-0 


3.1 x 10*' 3-6 

' <?w ':ork 

i 

10.1 x 10*’’ 9- 3 x I0’ 6 i.l x 10*** 1.2 x 

10* £ 

7.2 x 10** 8.3 x 10* 3 

:j. y. 

ii 

3.2 x iO* 2 2-9 x 10* U 3.5 x 10* : 3.5 x 

10* 1 

2.e x i:*~ 2.6 x 10*-' 


^ C |i ^ ^ # 

. prir.gf ield I 1.2 x .10* 3 l.j x 10* ' 1.6 x 10* 1.6-x 10* 8.5 >. 10*'' 9-9 x 10 “ 

X « p i - *. m % 

vi 9.3 x i.0 8.5 x .0* i.l v 10*- 1.1 6.5 x 10** 7-6 x 10 * 
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By employing the aforementioned characteristics and variables. It la 
believed that activity concentrations will be reached that vould approach a 
realistic case. 

2. Calculation of Ac* Wtty Concentration 

The calculation of concentration of activity used In the preliminary 
study Is based upon n cons^ration of the area of the water surface, the volume 
of vater, and the total watershed area. In the case of Duluth, Minnesota, where 
the entire Municipal water e'ipply is obtained from Lake Superior, only the depth 
of the lake a. the Intake point was considered. 

Current Information on the characteristics of the watersheds was obtained 
in each case from officials of the municipal vater supply system Involved or 
their engineering consultant. Data used in this report lo tabulated In Table VI. 

TABLK VI 

Characteristics >» f Watersheds Serving Selected Cities 

Total Mutorched Total Water Total Water 


Denver, Colo. 

Area 

jf 

hi1 

Surface ( 

ft. 2 ) 

Volume 

(liters 

8.29 

X 

10 10 

U .36 r 

V? 

2.34 

X 

10 U 

Houston, Tex. 

7.84 

X 

10 l ° 

7.85 x 

10 8 

2.50 

X 

10 1L 

New York, H. Y. 

4.24 

X 

10 10 

1.37 x 

10 9 

1.90 

X 

10 12 

Springfield, Mass. 

1-30 

X 

10 9 

1.68 x 

10 7 

1.97 

X 

10 10 


Ac an example of »•!•« procoduro used, the derivations of the concentra¬ 
tions in the Denver, Coloi**' water supply prior to any municipal treatment are 
presented here in detail. 'iV other cities listed in Table TV vere analyzed In 
the same manner us Denver, »-n:ipt for Duluth which is presented separately. 
Danver Case I 

In case I, only radioactive fallout deposited directly on water 
surfaces was considered ai " ,u «’ urco *-he water contamination. Tho following 
assumptions were made prior 10 *•*«« calculation. 
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1 . 5 HT bomb attack. 

2. Vlnd speed of 15 mph. 

3 . ->n intensity of 5,500 r/hr over the entire watershed. 

4. Complete homogeneous solution of the soluble portion of the fallout In the 
voter within one hour. 

5. A ratio of ooluble radioactive ato ns per square foot to intensity was 
chosen as that value for each isotope at 51*5 miles directly downwind. 

C. That there will be no change in concentration of radioactivity either duo 
to reservoir draft or feed-in. 

7. That only Sr-89, Sr-90, Ru- 106 , 1-131, Cs-137, and Ba-llO vould be conoldcred. 

ft. That the volume and area of feed stroans will not affect the final concentra¬ 


tion appreciably. 

The coluble activity density, for each of the six isotop co being cod- 
flldcred, in atoms/ft 2 waa calculated as follows: 

Given: l(l) > 5,500 r/hr, 


• 6.17 X 10 atomo/ftyr/hr. 


fy.(90) iq , 2/ a. 

a 11.1 x 10 atomc/ft Jr/hr, 

JI.(106) 10 21 . 

" ^.87 x 10 atoms/ft Jr/hr, 

<( 137 ) 10 . 2/ , 

- -y^yy - a 10.7 X 10 QtOKs/ft Jr/hr, 

< (l37) 10 2 /, 

■ ■ a 9.19 X 10 atonc/ft j r/hr, and 

N. (l**0) i_o , 2 1 . 

— a 12.- x Kj atoos/ftyr/br, 

where the subscript 4 is the value of the particle site parameter, a, correcpondirg 


H u U37) 10 

' -y^yy - a 10.7 X 10 at OK. 


to the downward distance. 


By multiplying 5,500 r/hr by <(A)/i(l) for each isotope, the total 
number of soluble atoms of each isotope per square loot, K^(A), will be obtained. 
This value for each isotope was found to be ns follows: 

<(A) - KD 


<(89) 

m 

3.1*0 X 10 l4 Btoms/ft 2 

<(*> 

m 

6.10 x 10 l4 (4toms/ft 2 

<(106) 

• 

2.68 x 10 l4 ntoms/ft 2 

<( 131 ) 

■ 

5.e9 x 10 l4 ntoms/ft 2 

<( 137 ) 

• 

5.05 x 10 l4 ntoms/ft 2 

<( 140 ) 

■ 

6.82 x 10 l4 atoms/ft 2 


Multiplying the <(A) value for cocli of the isolopt-s by the total area 
of water surface will give the total number of atoms, N(A), of each isotope in 
soluble form, deposited in the water. 

<(A) x Area . HU) atoms 


N(89) 

■ 

3.1*0 x 10 l4 X <1.36 x 10 8 

■ 

11* .0 x 10 22 

N(90) 

• 

6.10 x 10 l4 x ( |.36 x 10 8 

■ 

26.6 x 10 22 

K( X06) 

■ 

2.68 x 10 l4 x »*.36 x 10® 

« 

11.7 x 10 22 

N(130 

9 

5.89 x 10 l4 x ( t.36 x 10® 

m 

25.7 x 10 22 

N(137) 

m 

5.05 x 10 l4 x ><.36 x 10® 

m 

22 

22.0 x 10 

N(lUO) 

m 

6.82 x 10 l4 x <*.36 x 10® 

m 

22 

29.7 x 10 


atone 
atoms 
atoms 
atoms 
atons 
a tons 


Thus Knowing the total number of soluble storm of each isotope; divlolon by 





tba total vater roliaae, VV In liter* will giro the concentration of each Isotope 
la atooa per liter. 


C, m ll( a) atoms 
V liter 


* 

Ca 

9 

14.8 x 

10 22 

• 

6.33 X 10 11 

9? 


2-3*- x 

10 U 



# 

9 

26.6 X 

10 22 

m 

11.4 x 10 11 

90 


2.34 X 

To 53 - 



* 

c ^ 


11.7 X 

10 22 

9 

4.99 x 10 U 

106 


2.34 X 

To 11 



* 


25.7 X 

10 22 


.. 11 

C 131 

■ 

2.34 X 


m 

11.0 x 10 

c* 


22.0 x 

io 22 

m 

9.40 x 10 11 

137 


2.34 x 10 U 



C *| 


29-7 x 

,«22 

10 

m 

12.7 x 10 11 

l40 


2.34 x IQ 3 * 




ter 


< J 


We may now convert the activity In atoms/liter to mlcrocuries/mllliiiter 
(pc/nl) tiding an appropriate conversion factor for each Isotope. 

First, to change from atomb/llter to disintegrations per minute per liter, 
d P°? > Me following relationship Is used. 

°« * C A X ( 22 ) 

Where: C n - Concentration of activity In d ym/t. 

C A “ Concentration of activity In atoms/i, and 

^ * Characterlstlc time constant defined by the e<juatioo 

\ m 

V2 

Where: • physical balf-llfe of tb# specific Isotope In Mnutes. 
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1 ? , 

Site* there ere 2.22 x 10 ' dpn/curle, a constant for all Isotopes, 
division of C Q 'ey this valu* vili give the concentration In curleo/llter. 


C A* 


f 2.22 X 10 U 


2.22 x 10 




where Is th= concentration of radioactivity In curies per liter. 

Multiplication of hy 10^ will convert curlcs/liter to mlcrocurie.?/ 
milliliter. To* entire conversion factor for etch Isotope la then 

X. 


K. 


2.22 x 1CT 


The conversion factor, X^, for ench of the six isotopes being studied 


ate; 


Kq 9 . 4.25 x 10' 15 


90 


^106 


K 13l 


2.12 x 10 


- 5-95 x 10 


-17 


■16 


2.72 x 10 


-14 


K L37 • 1.93 x 10 


K 140 " 1 ' 70 x 10 


■17 


■14 


By Eulsiplylng each inversion factor hy the respective concentration 
of each isotope ic atoms per liter, the following values were obtained. 

0(69) - 2.7 x 10*3 uc / al 
C(50) - 2.4 x 10* 5 w c/ml 
C(l06) • j.C x 10* 1 * uc/ml 
C(l3i) - 3-0 x 10* 2 itc/tal 






Case II 


* 

. 

C(137) - 1.9 * 10* 5 Hfi/ml 

C(lUo) - 2.2 x 10* 2 uc/nl 

In this case the soluble fraction of the fallout particles deposited over 
the entire watershed was assuracd to have been transported to and homogeneously 
dispersed in the streams and reservoirs. Other assumptions were the same as 
for Case I. To continue the Denver example. In Case II, the watershed area of 
8.29 x 10 10 ft 2 was multiplied by the IT^ value for each isotope to give the 
total number of atorto of each Isotope. The remaining calculations were carried 
out in the sane fashion aa for Case I. 

Duluth. Minn. 

The entire water supply for the city of Duluth Is obtained from Lake Superior, 
and hence, this situation la different from that at the other cities studied. The 
Duluth Intake line extends out llj6 feet Into the lake and terminates at a lake 
depth of 72 feet. The fallout deposited over this Intake was considered to be 
evenly dispersed throughout this depth of water, and from this assumption the 
number of atoms per unit volume and subsequently the activity In yc/ml was calcu¬ 
lated. 

j. Conclusions 

In all of the caocs thus far studied, the time factor has not been considered. 
Even at on assumed wind velocity of 15 miles per hour, a number of hours would 
have to elapse before the arrival of particulate mntter at the far reaches of a 
diotant watershed. Even if the material landed directly in the water, much more 
time would be involved in transportation In streams and through reservoirs to the 






-181- 

intake of the water wOi’Ks. There would be further delay la getting hnto streaas 
thooe particles deposited Initially on diy surfaces. Another delay would be In 
the occurance of rainfall of sufficient intensity to dissolve and/or wash this 
naterlnl into streams. All of this time will alio, radioactive decay to proceed 
thus ameliorating the hazardous effects of the fallout. 
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I. De contamination of Water Supplie r* 

1. Removal by Conventional Wp.ter Treatment Processes 

Conventional municipal water purification processes generally Include 

aeration, chemical coagulation with sedimentation, rapid sand filtration and 

(79} 

chlorination' . To a lesser extent, lime-cods, osh softening and :,lov rand 
filtration are used. Fxcept for c-eintlon, the processes are capable of removing 
radioactive contaninatlon to some decree, either singularly or in combination. 

The decontamination capability of each type of process is discussed separately 
below, with particular reference to the six elements, barlun, cesium. Iodine, 
lanthanum, ruthenium and strontium. 

a. Chemical Coagulation 

The most common coagulants in water treatment, aluminum and Iron salt*, 
form aluminum or ferric hydroxides which precipitate as floe. The chemical 
floe acts as an efficient scavenger by ndcorblng, entrapping or otherwise 
bringing together suspended matter, particularly that which Is colloidal in 
nature. The artificial increase of the alkalinity in water may also form the 
hydroxides of heavy metals, which co-precipitate with the aluminum or ferric 
hydroxide. 


ORNL 


( 80 ) 


Coagulation, followed by sedimentation, has been extensively studied by 
with rcp-'i-tcd removals of only 36 $ and 51$ for Cs and Sr, respectively. 


It appears that the method Is effective for the removal of suspended or colloidal 
material ar.d for most cations of valence 3, •* or 5, including the rarc-earths 

group. Removals In excess of 98 $ have been reported for r- 32 as the 

. . , < 8 l)( 82 ) 

orthophosphate' ' . 

(83) 

Removal of fa, I and La by coagulation was investigated by Lacy , 
with a maximum of 70.7$ found for both Ba and La, but only ^5$ for I. 


References begin on page 1J5 
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Matsuraura^®^ 

coagulation. 


states that Ru can be removed up to JA by chemical 


b. Rapid Sand Filtration 

The amount of radioactivity removed by filtration vlll vf.ry dep» r.ili'.i! ur. 
the nature of the material. ^ cad Zr, probably present in the colloidal ~ta*e, 
veie removed from up to 935 “ by sand filtration alone as reported in a com- 
preher.slve C.RTIL report^® 0 ^, whereae materials such as Sr ) and C~ (’r-Oi) 
present in true solution, were cot greatly reduced by passage through siund 
filters. They further reported that Pa and la could be removed up to sv.d 
iH, respectively. 

No data were reported on the removal of ! or Ri by rapid suad 
flit ration. 


c. Chlorination 

Hannah, et studied various methods for the removal of 1*131 

from water, and found that email dosages of chlorine in the presence of 100 ppm 
activated carton produced up to 8056 removal. The authors^®'’^ concluded that 
the only effective method found for removing 1-131 with materials normally 
available in water treatment plants involves chlorination followed by adsorption 
of liberated iodine on activated carbon. The o ptimum chlorite dosages were 
quite small, 1,c., 0.05 to 0.1 ppm. The removal of 1-131 decreases to less 
than ZOf> when the dosage of chLorir.e increases to ippra. Therefore, normal 
pre-chlorination employed by the water treatment plant could not be used for 
lodlr.e removal because the chlorine residuals would generally exceed i no re¬ 
quired dosage for activity removal. Static iodine, in dosages greater than 
0.01 ppm, Inhibited removal of 1-131 with chlorine ar.d Ajua liuchiar A 
(activated charcoal). Variation of pH, achieved by adding sulfuric acid and 
rodlun nydroxlde, was found to have little effect on the removal of iodine 
by chlorine and activated charcoal. 








d. Lime-Soda Ash Softening 

Softening plants remove excessive amounts of scale-forming, soap- 

consuming compounds, chiefly comprised of th* cations of calcium and magnesium, 

by the addition of lime and soda ash which precipitates calcium as a carbonate 

and magnesium as a hydrate, Because of the similarity In the chemical picperties 

between strontium and calcium, the former ion is co-preclpitated. For most 

h 5 

satisfactory removals of strontium, Poyl found that excess dosages of both 

, (g„ ) 

lime and soda ash were required. Removal efficiencies of 99*7? are 
possible under favorable conditions. Other Isotopes which can be effectively 
removed by the lime-soda ash are: 3a, Cd, f, Sc, Zr and Nb. The removal of Ba 

/p- \ /r a \ 

(99$) and La (90$) was reported by McCauley, et al 1 , 0FC.L reported 

that Cs could be removed by lime-soda softening, tut not efficiently (25$). 

e. Slew Sand Filtration 

Downing, et &1 reported that the slow ^and filter was very 

effective in the removal of radiostrontiutt in the first few days of operation, 
but that the rff iciency decreesed rapidly and essentially drops to zero at 
days. In another paper, presented by Eden, et al , on experiments carried 
out in the sane Water Pollution Research Laboratory In Britain, the slow sand 
filter appeared to be effective for only a few hours. The activity of the 
effluent rose steadily and reached 30 $ of the initial (unflltered) water in 
one day, 70$ in 2 days, and 95$ after 7 days. Eden, et nl ^ also reported 
that iodine was not efficiently ( 50 $) removed by the slow sand filter. 

2. Removal by Non-Convent lonal Treatment Methods 

Below are listed several methods of decontamination of water not 
commonly employed in municipal treatment plants. Economic considerations nay 
limit the use of many of these methods, except in emergency cases, 
a. Ion Exchange 

A few water treatment plants use loc exchange as a softening measure. 
Resins have been found to provide one of tbs cost effective methods of 





decenttnination. Over r /)-9f> pc-bcvsI of 3a, C&, Cc, Cs, I, La ead over 99^ 

reaeval of V and Y,r vuo i 4 uotad by OK r u ^ | ar«l Airphlett and S&mmoa ^ ^ 

U s 

reported vork using nix-d bed** A derontamiration factor Of 10 to 10' may be 

fe- ) 

achieved if Indulge ie «. j i situated . The natural greens an da, used In 
tome municipal coitculrtg Mietcllaticne, nre cot effective La the removal of 
anions. 

Although the Ion oxtfi&uge process offers one of the isost promising 

methods for rc-jcvlng tv 1 hi-contaminanta, the coat may preclude widespread 

application unless sulLetlo regeneration techniques ore developed. 

The po. sii.-mty of using hone vatcr softeners, of the ioa exchange type, 

A>? 5 

has teen discurru-d in th- literature v ’’ 
b» Phosphate Q.v ij’-il atloa 

Seca-use the conventional method employing alum or Iron salt as the 

(o'* ) 

coagulant does not affectively remove Sr and many other Isotopes, Lauderdale 

Investigated phosphate o?( l fi' 1 l«tton as a mean* of increasing removals. A more 

iy ) 

comprehensive study yoe ir*»'!v later by h’esbitt, et ai ^ . Based oa the 

f s' 

theoretical conslderatlonf*, ftfl veil 05 experimental observations, Nesbitt ?y 

concluded that the me.d.or 1 t*mn of the removal of rudloactivity arc ccprccipitatlon 

and adsorption. Ke.-.bitt C' p.'.’tcd a strontium removal. Udie three variables 

vhlch wort found t.c exert Mi® greatest influence on cocgHstion «re pH, 

calcium cor.centr-xtion, w.d j'ho'iphr.t.r concentration. Careful control of there 

vuriablto vac necessary t.c* achieve optimum strontium rcaovaln. 

(yi) 

Tan’ira and StruM n f b reported that low cost natural phosphates 

could be uroJ in plane of phosphate. The authors reported that over 

yCji removals vore (•chic*.'-.’' 1 Wth these minerals. 

Removal of (Jv.?'/, fVi :-n, r /}f> for I/i rr.d for Ru have been reported 
by Katcunrum \ Cow-rel- a - ^ ^ report only yrft removal of Cs using 


phouptiate coagulation. 



At leaut two modified coagulation processes have teen studied, o« of 
„H‘h Is the continuous addl.tion of city. One modification of the piccese la 
the addition of small amounts of silver rltiate to increase the rcmovil of 
lc-11dc. The removal of 95-5$ 1-131 has teen reported by Eden, et al y . 

Lacy and de lagu"' rcpoi’t rcocvnlc of 9?-7$ for Ba cr.d Ce, and 

. 97 ) 

yy,li$ for Su when coagulation is supplemented by cloy olditlon. Cowsor, .et al 
u.ijr.g a llme-c'oda ash softening process with the addition of ciuy, reported 
t .t'/Kenlum removals of 76$- The removal of cesium averaged 86$ with clay doses 
of "'00 ppm. 

4. Metallic Dusts 

lacy' 1 ' 1 ' conducted laboratory Jar-teet studies for the removal of 
radioactive contaminants by employing various concentrations of pcnnlered 
nl'/Tilr.um, copper, iron and sine. la general, highest removals were obtained 
vlth Iron powder. Removals of 37-2$ of Z ve.re reported In 19 &, a 

pi.t ciit van grented to Sllker ^for removing ?. As, Kn, the rare earth metals, 
and e.ctlr.ides from aqueous eolations by sorption on particles of aluminum metal. 

e. O lay Materials 

Straub, et al and others v l0t v .1C 1 *) the .^ c of various 

t.ypeu of obeys (kaolinite, montmoi'illor.ltc nrid shales) for the removal of 
rpu'-lfi'.’ radioisotopes from water and vacte solutions. The clays vere found 
to he especially suitable for the removal of Cs, 7. r and lib. With the addition 

/ Qa\ 

of ,COO ppm of clay, removal of over 9&$ of t’s vao report’d ' v . Amphlctt ar.d 
S.-jiiium ^ reported 93$ removal of Sr. This method r.ot be economically 
j u tll'icd because large volumes of clay must be handled both initially and as 
c. ■•.i;.huni:mtcd sludge. T':ii3 method has its greatest potential use where 
gr >*'iglc ar.d hydrologic conditions arc such that the ground itself may be 
v»i for disposal of the eludge. Studies have also beet made on the 

ujo of clay for high level vaetc dispoc:d by fixing specific radioisotopes 
t„' o clay by means of .lgn temperat-ir- trevtscent. 
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f• Oxidation-Deduction 

This method was especially designed to improve !.u removal. Runhenlum 
Is often present in the fora of a nitrosyl-ruthenium conplex which ia difficult 
to remove ty coagulation, procedures. An approach to the problem, described 

/ qjA 

by Amphlett and Samson ' , involves pre-treatment v ith an oxidising agent to 

destroy the ruthenium complexes, iolio.ci uj a calcium phoephate-ferrouu 


phosphate floe in the presence of to maintain she ruthenium in a 

reduced form. As high as 99^ of ruthenium was removed when 100 ppa KaHSO^ was 
used. 


g. Foaming 

A foam separation technique, based on surface activity phenomena, waa 

recently developed by Radiation Appli eat tons, Inc., Long Iolund City, Lew Yorh^ 0 ^. 

A reduction of the concentration of radioactive strontium-yO and cesium-137 
-5 -12 

from 10 to 10 molar in lav-level nuclear waste was claimed. The concen¬ 
tration of the foaming agent is one of the most critical factors in determining 
the separation factor (surface adsorption factor). A detailed discussion of 
the theory as we LI ao experimental observations of the foam separation technique 
was presented by Schoen and Mazzella ' . 

Laqy'i^ Investigated the flotition process in 1957 and reported 
removals of up to 78 $ for oa-lbO-La-lbO, 8[# for Cs-137> 7*+$ for Ru and 68 In for 
Sr-90. These studies were conducted using cetylpyridinlum chloride bg the gut- 
face-active agent and Teaming by aeration. 

3. Kncrgency Methods for the Dec or. Lamination of Radioactive Water 

According to the literature, very few emergency decontamination unlt 3 
are now available for municipal use, although the U.S. Arsy has developed 
several mobile units which could find vide application for public use. Several 
s mall er decontamination unit 6 are corrercially available, most of which employ 
ion exchange techniques. 







Supplementary treatment to the conventional water treatment processes 


would probably he the cost efficient menns o:' supplying potable water to a 
larje number of people, following radioactive contamination of the water. 

a. Municipal Gita teconts'B.tr.fttton Unite 

* A: o) 

Woodward and Pnbeck v ' report on ar. ion cjcehcrge process tnat mr.y be 
u£ti to supplement the normal water treatment procedure. An Ion exchange 
cclvnui Le the form of a cartridge is inserted into .10 system, and idl.:paced 
of after 'n-eriX'/hroo/h occurs. The type cf rerIn v-lli vary according to the 
time after detonation that the water will bo p-t into use. A mlx>-1-bod :-.-cnin 
’•til be necessary for irmrediate u.-o of tbo voter, whereas a notion rvslu will 
.'j'fl'.e for long-term usage, due to the decay cC short-lived iodine-131. 

Tee life of t-a:h lor. cartridge will deptri upon ',‘c total sciido present. 

Although large, cent rally lo-uhed, da font xninat ion ur.it: an: more 
c:jr.omical, smaller decentralized units would b? preferable os the problem 
cf water diatributioa .xmld be lessened. 

b. hie Id Decontamination Til ts 

fti 1) dlfl 

JThe 1.3. Ar*v has developed e:\v~rnl mobile decontamination 

utits prlnm lly for use by troops in Mv* field, but these could also be -.mr.i to 
s.pply water for small population group: in case of nuclear attack. One such 
unit is comprluod of a flocculator, filters, dual-bed Ion exchange column 
and a cblc-rlnnfcoi, 'ihe entire unit is mounted or. two 2-1 fi ton trucks. A 
1,500 gpfc output la obtainable, regeneration ol the ion exchange l-eclno la 
effected by hydrochloric acid and soda ash. 

c. Mlr c -llftntons lei or. bend nation 'Jr.lt3 

laaotts* Ohenlccl Company ha: Ic.'eiOpvi a mixed ion exchange jvcfn 
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vater. »e -jolt coeulrif of a ted of Ariterlit- realn !3-3 ClH-120 and J/'A-kiaJ 
and la capable of producing 10 gallons of vater of triple distilled quality. 


A -f.lt dt.-.eloped by Lauderdale and Etco ri 


consists of tv!(i 


columns of steel root,, burnt clay cxd activated carbon In one col—a a:.<1 a 
ul.-i-l bed ion exchange resin in the other. Approximate!;.' jO Ute.-o cl' outer 
can be dei:c>>A«.,<no*«^ V ith as exceeding 99.9$. 


N..-CSO 


< 5 ii. lA-y ~ r\tzz rite e rr/^o-J cf rrtHo^ciive vntci' >i<?con- 


taainatlon whereby the water is passed through a column of natural. rv.Wtalo, 

.-,ueh as clay, leaves, huff'^, gravel and sand. Pesovals of over 9^ were 
reported. O.e problem Wild be to obtain ur.cor.tuainate-d materials from the 
environment for use In ouch a col-nan. 

**• Diaeussian and ConeImlonB 

— e more important radioisotopes that appear initially it ccritOfiMated 
vc.ter are Ba, Cs, I, La, Ru ar.d Sr. T.ot>oratoiy data indicate that the 
combination of conventional water treatment nvithods, nt-t.’ly, chemical coagu¬ 
lation, lime-coda ash soft ning, and sued filtration, Is capable cf rcw>/lng 
59$ ctrontium. The rctayval of iodine con be accomplished by the addition 
of silver ions. 86$ cesium can be renu>vcJ by the llme-.-oda ash softening procedure 
when proper amour, to of clay are added. 

'fhe efficiency of gross activity removal lice generally tetve'.'i '/) and 
75$ in plant scale cporntlo:.***^. The lower plant cce.ie efflcienolcc, In 
comparison to the Irtborntoiy Jesuits, are bvll-.-vod to result from the rant 
t.lvit; (a) only one or two decontamination procc-ii;en, .-jsccg those Ber.t'ovid, 
tre employed by the ur.nnl vater trent sunt plart, aul t.hc operating ci/vU*. Irmti 
m.-\y not be optimised; aj-.l (b) in the laboratory, the radioactive rsitarial > arc 
generally present lr. the form of simple salts. In contrast, the rad lot." Uve 





*15''• 


materiad-i reaching the vcter treatment pit et are generally believed to be ir. 

ehesieel forms which are mors difficult to remove beeducc the portion which ifl 

readily remove! has already been eliminated by natural decor.taetnation 

(llA 

mechanisms before reaching the water plant. Setter and .boa 3 ell'' stated that 

the longer the activity has been subjected to natural purification, the more 

difficult it is to remove the remaining soluble activity by coagulation. 

~b rho'sld b-c pointed out here that the decor.lamination data reported 

ir the literature, both on laboratory and plant scale, ar« based on or.o of the 

follcviug three scarcer of radioactivity: 

(l) added radioactive salts 
( 2 J lev-level radioactivity waste 

(3) long-range (vorl’.-vide) fallout from atomic bomb teats 
Che chemical and physical characterlotlea of radian chive materials fw.t 
any of the above sources is different from that proceut ir. local fallout. The 
latter is characterized by Its lor water solubili ty of 1 $ to 23 C in comparison 
to the or higher solubility of radioactive materials cf the above three 
categories. Low water solubility generally will incresr.e the removal of local 
fallout by conventional water treatment mtbods which ara efficient in dealing 
with particulate or colloid materials. 

The atxUnuin decontamination factor for etch of six selected elements 
above ere listed in Table X. Thcce values arc baced on the best available data 
in the literature. It is emphanieed that the percentage of removal Is greatly 
dependent on the chealcnl form and physical state of the radioactive element as 
well 2.3 the concentration of the treatment additives and other related conditions 
such ar- p'i ml temperature of the water. As an example of the decontamination 


factors required to reduce activity to the Mc’J, based on hypothetical 
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Tnble X 


Mnxlmun Decontamination Kac‘ors 

Peported 

in the Site rv 

-r? for 

Selected Iso'-rc* 

Decontaai*. 
nut too. 

Method 

Decontamination Factors (D.F.) and Pe-f 
B'i-ltO Cs-137 1-131 fa-lltO 

* recces 

106 

Sr-89,90 

1. Che-raical 

3-4 

1-5 

1.8 

3.4 

10 c 

2.0 

Coagulation 

(83) 

( 80 ) 

(83) 

( 8 f) 

(84) 

(So' 

2. Llme-3od» 
Softening 

io a 

(87) 

(83) 

m 

10 l 

(87) 

ifR 

3-3 x i; L 
(So) 

3 . Rapid Sand 

PiItration 

2 x 10 1 
( 80 ) 

(30) 

m 

3-8 

( 80 ) 

SP. 

1.04 

(So) 

4. Ion 

10* 

10* 

10 3 

10 3 

ICR 

3.3 x 1: 3 

Exchange 

( 80 ) 

( 80 ) 

(9i) 

( 80 ) 

( 80 ) 

3• Phosphate 

3.2 

1.5 

NF 

5 x 10 1 

■c-7 

4.5 t i? L 

Coagulation 

(84) 

(9?) 


(84) 

(34) 

(94) 

6 . Cloy 

ifR 5 

X 10* 

(fio) 

tf) 

m 

ICR 

2 x 10 1 
(9U 

7- Metallic 

Dusts 

NR 

hr 

1.6 

( 80 ) 

NR 

ICR 

HR 

8 - Coagulation v/ 
Silver Ions 

jrR 

HR 

2.2 

(96) 

HR 

SR 

NR 

9* Oxidation- 
Reduction 

im 

HR 

NR 

HR 

10 2 

(91) 

HR 

10. Coagulation 
with Clay 

3.3 x 10 l 
(99) 

3-3 x iO 1 
(99) 

' HR 

nr 

3.3 

(99) 

HR 

It. Lime Soda 
with CLay 

nr 

7.1 

(97) 

IfR 

NR 

4.2 

(97) 

HR 

l?. Foaming 

4.5 

6.7 

HR 

•*•5 

3.9 

3.1 

(Flotation) 

( 109 ) 

(109) 

U09) 

( 109 ) 

( 109 ) 

13. Glow Sand 
Filtration 

NR 

HR 

(90) 

IfR 

SB 

HR 


Notes; 

(1) Percent Removal » 100 . (gi--)ioo 

(2) RtJ - not reported iu literature revieved 


. » 1 





■ 'stji % j >' 











of the Springfield, Haas. vatershcd (aa reported eloevhera la 
ir.it rrpiirO, Table XT vna prepared. It la apparent from thin tools that ion 
Grange u. the only £l_^le decontamination procesa that will adequately remove 
■<v* alnlsna. concent rati on reported for the six celeettd Isotopes* and that no 
pro-'eaa .111 adequately remove the taxlcroa radlocontaminant, namely I-I. 31 . 
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